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The  application  of  the  theory  of  collectively  compact  operator 
approximations  to  Fredholm  integral  equations  gives  theorems  on  the  uniform 
convergence  of,  and  error  bounds  for,  the  approximate  solutions  of  these 
equations  obtained  by  a  direct  method  of  numerical  integration.   It  is  the 
purpose  of  this  paper  to  investigate  the  possibility  of  extending  this  concept 
to  linear  singular  integral  equations  in  a  Holder  space. 

In  this  paper  it  is  shown  that  a  sequence  of  approximate  operators, 
obtained  from  the  singular  operator  by  a  sequence  of  interpolatory  quadrature 
rules,  converges  pointwise  only  on  a  certain  subset  of  the  Holder  space. 
Hence  the  convergence  of  the  approximate  solutions  of  the  singular  integral 
equation  obtained  by  a  direct  method  cannot  be  guaranteed  by  the  collectively 
compact  operator  theory.   However,  it  is  shown  that  a  direct  approximation  of 
the  reduced  equation  does  satisfy  the  hypotheses  of  the  abstract  operator 
approximation  theory  and  that  the  sequence  of  approximate  solutions  converges 
in  the  desired  space.   It  is  concluded  that  this  is  an  appropriate  method  of 
obtaining  approximate  numerical  solutions  to  a  linear  singular  integral 
equation. 
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1.      INTRODUCTION  AND    SUMMARY 

1.1     Statement   of   the    Problem 

This    paper    is    concerned  with   a  method   of   obtaining   approximate 
solutions    of    the    linear   singular    integral  equation  discussed  by 
Muskhelishvili  [33].      This   equation  has    the    form 


a(s)cp(s)   +  J*  k<£>st)  cp(t)dt  =  f  (s) 


where  L  is  a  smooth  line  and  a(s) ,  f(s),  k(s,t)  are  given  functions  satis- 
fying a  Holder  condition. 

There  are  a  number  of  references  to  the  approximate  solution  of 
singular  integral  equations  (for  example  [18,  p.  167],  [23,37,41])  which 
generally  involve  methods  (such  as  the  method  of  moments)  carried  over  from 
regular  equations  to  singular  equations.   However  these  methods  are  not  always 
the  most  convenient  for  numerical  computations.   This  paper  investigates  the 
problem  of  extending,  from  Fredholm  integral  equations  to  singular  integral 
equations,  the  method  of  approximating  the  integral  equation  directly  by 
numerical  integration. 

Consider  a  Fredholm  equation  of  the  second  kind 


1 
x(s)  -  '  k(s,t)x(t)dt  =  y(s)      ,  O^s^l 
0 


where  y(s) ,  k(s,t)  are  continuous  real  or  complex  functions  for  O^s ,  t^l  and 
a  continuous  solution,  x(s),  is  sought. 

In  a  classical  method  of  approximate  solution  based  on  numerical 
integration,  this  equation  is  replaced  by  an  algebraic  linear  system 


n 


x  (t  .)  "  2  w  .k(t  .,t  .)x  (t  . 
nv  niy    ._.  nj  v  ni'  nj'  nv  nj 


B  y(tni)     i=l,...,n. 


If  this  system  has  a  solution,  then  an  approximate  solution  x  (s) ,  O^s^l, 
of  the  original  integral  equation  can  be  obtained  by  interpolation. 

The  technique  of  replacing  an  integral  equation  by  a  finite  system 
goes  back  at  least  to  Fredholm,  one  of  the  founders  of  the  theory.  However 
there  arises  the  question  of  the  relationship  of  this  approximate  solution, 
x  (s) ,  to  the  exact  solution,  x(s),  and  the  nature  of  the  convergence  of 
x  (s)  to  x(s).   Hilbert  gave  convergence  proofs  for  approximate  solutions 
using  the  rectangular  quadrature  formula.   Other  convergence  theorems  and 
error  bounds  for  approximate  solutions  in  particular  cases  and  under  various 
assumptions  have  been  obtained  by  a  number  of  investigators  including 
Brakhage  [ 10]  and  Kantorovich  and  Krylov  [26].   Specific  references  to  other 
investigations  can  be  found  in  Anselone  [5]. 

In  a  series  of  papers  [1-8]  beginning  in  1964,  P.  M.  Anselone  and 
a  number  of  co-workers  were  able  to  abstract  the  essence  of  these  results, 
and  to  develop  a  general  operator  approximation  theory  for  linear  and  non- 
linear operators  on  an  arbitrary  Banach  space.   The  operator  equation  is 
regarded  as  being  approximated  by  the  same  equation  with  an  approximate 
operator  replacing  the  exact  operator.   It  is  then  shown  that,  under  suitable 
conditions,  the  convergence  of  the  sequence  of  operator  approximations  implies 
the  uniform  convergence  of  the  solutions  of  these  approximating  equations  to 
the  solution  of  the  original  operator  equation. 

This  collectively  compact  operator  approximation  theory  is  then 
applied  to  a  wide  variety  of  integral  equations  [5].   Anselone  [2,4]  showed 
that  the  abstract  theory  could  be  applied  to  an  integral  equation  with  a 
"mildly  discontinuous"  kernel  by  considering  the  analysis  in  the  space,  R, 


of  bounded,  Riemann  integrable  functions  with  the  supremum  norm.   Atkinson  [7] 
further  extended  the  class  of  integral  operators  to  which  the  abstract  theory 
is  applicable  by  treating  the  case  of  an  integral  equation  with  a  weakly 
singular  kernel  (that  is  with  an  integrable  singularity).   Anselone  (in  a 
lecture  and  also  in  a  private  communication)  suggested  that  the  collectively 
compact  operator  theory  could  also  be  superimposed  on  approximation  techniques 
for  singular  integral  operators.   Such  an  extension  to  the  case  when  the 
kernel  is  strongly  singular  involves  the  computation  of  singular  integrals 
which  are  interpreted  in  the  sense  of  the  Cauchy  principal  value.   Gakho v 
observes  [18,  p.  156]  that  methods  of  approximate  computation  of  singular 
integrals  have  not  been  worked  out  in  detail  so  far,  either  from  the  theoreti- 
cal or  from  the  practical  point  of  view.   Gakhov  further  notes  118,  p.  167 J 
that  the  methods  of  direct  solution  (in  general  approximate)  of  the  complete 
singular  integral  equation  have  so  far  been  little  elaborated. 

The  theory  of  singular  integral  equations  has  practical  applications 
in  a  number  of  important  areas.   Among  these  are  a  variety  of  problems  in  the 
theory  of  elasticity  (Muskhe  lishvili  L33,  Chapter  13]),  the  problem  of  flow 
past  an  arc  of  given  shape  (Mikhlin  l 3 1 ,  p.  325],  Noble  [36]),  the  theory  of 
tides  (Poincare  [39])  and  more  recently  in  the  fields  of  dispersion  relations 
and  scattering  theory. 

1-2   Possible  Algorithms  for  Direct  Solution 

Any  such  algorithm  for  direct  solution  will  clearly  be  based  on  an 
approximate  quadrature  formula  for  principal  value  integrals.  References  to 
computational  experience  with  such  integrals  are  not  widespread. 

Voight  [46]  discusses  the  numerical  solution  of  a  Fredholm  integral 
equation  in  which  the  known  functions  are  themselves  specified  by  principal 


value  integrals.   These  singular  integrals  are  evaluated  by  using  a  Fourier 
series  technique  given  by  Collatz  [ill.   An  appropriate  transformation  followed 
by  the  expansion  of  the  numerator  in  a  Fourier  cosine  series  transforms  the 
integral  into  a  sum  of  principal  value  integrals  whose  values  are  known 
explicitly.   Voight  gives  no  discussion  of  the  conditions  under  which  this 
process  is  justified.   He  notes  that,  in  the  example  given  from  annular 
airfoil  theory,  the  computational  results  for  the  evaluation  of  the  kernel 
of  the  free  term  were  rather  unsatisfactory.   Noble  [36,  section  14]  outlines 
essentially  the  same  method  for  the  numerical  evaluation  of  this  singular 
integral. 

Delves  [14]  proposed  a  Newton-Cotes  type  method  which  is  basically 
an  application  of  the  approximate  product  integration  technique  of  Young  [47], 
which  was  used  by  Atkinson  [7]  in  the  case  of  weakly  singular  integrands.   The 
error  bounds  given  by  Delves  depend  on  the  numerator  having  a  sufficient  num- 
ber of  derivatives  and  on  certain  constants  depending  on  the  fixed  value  of 
the  singularity.   However  it  is  not  clear  that  such  bounds  will  hold  uni- 
formly as  the  singularity  varies  in  the  range  of  integration. 

Another  general  method  of  attacking  this  problem  is  the  method  of 
subtracting  out  the  singularity.   In  this  method  the  singular  term  is  subtracted 
out  of  the  integrand  and  this  integral  is  evaluated  analytically  while  the 
remaining  integral  (usually  with  an  integrable  singularity)  is  evaluated 
numerically.   An  example  is  the  modified  Gaussian  rule  mentioned  by  Delves 
[ 14] .   A  variation  of  this  idea  is  the  "method  of  symmetric  pairing"  suggested 
by  Bareiss  and  Neuman  [9].   In  this  method  the  function  values  of  the  inte- 
grand are  "paired"  for  arguments  that  are  symmetric  with  respect  to  the 
singularity.   This  again  produces  an  integrand  with  an  integrable  singularity. 


Finally  we  mention  Stewart's  method  [42]  which  can  be  regarded  as 
a  "mixed"  method.   The  singularity  is  subtracted  out  to  produce  a  weakly 
singular  integrand  but  then  the  numerator  is  interpolated,  which  is  essentially 
an  application  of  approximate  product  integration. 

Hence  the  available  methods  can  be  divided  into  two  basic  classes: 
(a)  Approximate  product  integration  techniques.   These  will  include 
approximation  of  the  numerator  by  polynomials  (Delves),  Fourier 
series  (Collatz)  or  some  other  appropriate  functions  such  as 
Chebyshev  polynomials  (compare  Elliot's  approach  [16]  for 
Fredholm  equations), 
(b)  Reduction  of  the  integrand  to  a  function  with  an  integrable 

singularity  and  the  application  of  a  direct  quadrature  formula. 
In  case  (a)  it  is  necessary  to  consider  the  convergence  of  the  sequence  of 
interpolating  polynomials  (see  section  2.2)  and  also  the  convergence  of  a 
sequence  of  principal  value  integrals  (section  3.4.2).   m  case  (b)  the  con- 
vergence of  quadrature  formulas  for  unbounded  integrands  must  be  considered 
(section  2.5) . 

In  the  application  of  these  quadrature  methods  to  the  direct 
solution  of  singular  integral  equations,  the  first  class  seems  to  be  the  most 
natural.   An  example  is  Gabdulhaev's  investigation  [  17]  of  the  method  of 
constructing  solutions  by  polygonal  approximation.   Gabdulhaev's  paper  does 
not  consider  the  convergence  of  this  quadrature  process  and  does  not  give  a 
numerical  example  to  illustrate  the  method.   However,  there  is  a  theoretical 
justification  of  the  convergence  of  the  approximate  solutions  of  the  singular 
integral  equation  in  a  Holder  space,  obtained  by  this  method.   This  justi- 
fication is  based  generally  on  Kantorovich 's  general  theory  of  approximate 
methods  [25,  chapter  14].   The  paper  and  its  implicit  assumptions  are  dis- 
cussed in  section  3.5. 


Bare  is s    and  Neuman  apply   their   method   of    the    second   class    to    the 
solution  of    the   homogeneous    singular    integral  equation   arising   from  Milne's 
problem   in   radiative    transfer,    but    the    authors    do   not    attempt   a   theoretical 
justification  of   their  results. 


1.3   Convergence  Theory 

It  was  noted  above  that  the  integral  equation  is  replaced  by  an 
algebraical  linear  system  by  means  of  a  numerical  integration  formula.   The 
solution  of  this  linear  system,  or  approximate  equation,  can  be  regarded  as 
an  approximate  solution  to  the  original  equation.   However  it  is  necessary 
to  ensure  that  such  an  approximate  solution  does  converge  to  the  exact 
solution  in  some  appropriate  sense.   In  order  to  avoid  the  necessity  of  con- 
sidering individual  cases,  unified  theories  of  approximation  methods  have 
been  constructed  in  a  functional  analysis  framework.   We  are  primarily  con- 
cerned with  Anselone's  collectively  compact  operator  approximation  theory  but 
also  consider  Kantorovich's  general  theory  of  approximation  methods. 

1.3.1  Anselone  Theory  [5] 

Let  X  denote  an  arbitrary  real  or  complex  Banach  space.   The  closed 

unit  ball  in  X  will  be  denoted  by  13    =  [  xeX:||x||^l}  .   Denote  by  [x]  the  Banach 

space  of  bounded  linear  operators  T:X"*X,  with  the  usual  operator  norm, 

i|t||  =   Jj|Tx||  .   Convergence  in  norm  of  the  sequence  of  operators  is  written 

||T  -t||— 0  and  the  notation  T  -»T  will  be  used  for  pointwise  convergence,  i.e., 

||  T  x-Tx||-0  as  n-*°  for  each  xeX. 
n 

In   the    standard   approximation   theory    (for  example,    Taylor  [43], 

or  Kantorovich  [25])    the   operator   T    is   uniformly   approximated   in  norm  by   a 

sequence   of   approximating  operators   T    .      Then  based  on   the  well-known  Neumann 

series   representation  for    (I-T)"1   the   following  basic   proposition   is   obtained. 


Proposition  1.3.1  [5,  prop.  1.4] 

Let  T,T  e[x]  and  ||t  -tII-'O.   Then  there  exists  T~  e[x]  iff  for  some 

N  and  all  n^N  there  exist  T   e [x]  bounded  uniformly  in  n,  and  in  this  case 

n  J  ' 

IIt        -T     1 1  — >0   as   n-00. 
ii    n  ii 

However,    the    condition  of   norm  convergence    is   not   met    in  many   cases 
of   practical    interest,    and   in  particular    in   the    case   of   the   approximation  of 
integral  operators   by  quadrature   formulas.      But  Anselone    shows  how  the  weaker 
hypothesis   of   pointwise    convergence,    T  -,T,    combined  with   certain   compactness 
conditions   can  be   used   to   deduce    an   analogous    theorem  giving   pointwise    con- 
vergence  of    the    inverse   operators    together  with  useful   error  bounds.      The 
following   concepts   are   needed  where   X  has   been   assumed   complete. 
Definition:  [5] 

An  operator   Ke [x]    is   compact    (or   completely   continuous)    iff    the 
set  YB    is   precompact. 
Def inition:[ 5] 

A  set  of   operators  ^r[x]    is   collectively   compact    iff   the    set 
3/3  =  {Kx:Ke?  ,xeZ?}    is   precompact. 

With    these   definitions,    the   development   of    the    theory  may  be    sum- 
marized by   stating   the    following   convergence    theorem  whose   proof   may  be   found 
in   the    papers   of   Anselone,    particularly   [3,5]. 
Theorem  1.3.2    [5,   Theorem   l.ll] 

Let   K,    K  e[x]    and   suppose    that    the    approximate   operators   K      satisfy 

1.  Pointwise   convergence,    K  ~*K 

2.  [k    :n^l]     is    collectively   compact,    which   together    imply   that 

3.  K    is   compact 

Suppose   also    that    (I-K)        exists    and   define 


8 


An  -||<I-KrVl|<VK)KJ| 


Then  A  -0  and  thus  for  all  n  sufficiently  large  A  <1,  in  which  case  (I-K  ) 
n  n  n 


exists  and  is  bounded  uniformly  in  n, 

-ll 


and 


x  _  ||  <X-K)"L!|  .l|Kny-Ky|l    +  Ajld-K)"1^! 


x„-x||    =  ||  (I-Kn)      y-(I-K)      y||   *  j^E '   ° 


n 

n 


1.3.2     Kantorovich  Theory  [25,    Chapter   14] 

Kantorovich's    general   theory  for  equations   of   the   second  kind 
projects    the   exact   equation  Kx  =  x-Hx  =  y    in   the    normed   space  X      (not 
necessarily  assumed  to  be   complete)    into  an  approximate   equation 
Kx  H  x-Hx  =  y   in  the   complete   sub space,   X,of  X.      The    theory   involves   the 
following   three   conditions   connecting  the   normed  space  X   and   the   space  X  of 
the   approximate   equation,   where  X   is    the   projection  of   the   space  X  by   the 
linear  operator  P,    and  H   is   the    linear  operator   in  X  corresponding   to   the 
general  operator  H   in  X. 

I:       (Condition   that  H  and  H  be   neighboring  operators) 
For  every  xeX,   ||pHx-Hx||    £  7]||x|| 
II:       (Condition  for   the   elements   of   the   form  Hx  to  be   approximated 
closely  by  elements   of  X).      For   every  xeX,    there   exists  xeX 
such  that  ||hx-x||    <.  T]  ||x|| 
III:      (Condition  for   the   close   approximation  of   the   free    term  of   the 
equation)      There   exists   yeX   such   that  ||y-y||    £  T]Jly|| 


The  following  theorem  gives  conditions  for  the  convergence  of  the 
sequence  [x  j  of  approximate  solutions  to  the  exact  solution  x* . 
Theorem  1.3.3  [25,  Theorem  3,  p.  549] 

If  the  following  conditions  are  satisfied: 

1.  The  operator  K  has  a  linear  inverse, 

2.  The  existence  of  a  solution  of  Kx=y  for  every  yeX  implies 
its  uniqueness, 

3.  The  conditions  I,  II,  III  are  satisfied  for  all  n  and 
11-0,  TIJIpJ-O,  Tl2||pJho  as  ir~, 

then  the  approximate  equations  are  soluble  for  sufficiently  large  n  and  the 
sequence  of  approximate  solutions  converges  to  the  exact  solution: 


|x*-x*||    £   Qo7]  +  Qj^IIpJI    +  Q2T]2l|Pnll    -   0   as   n— , 


where  Q_,  Q1 ,  Q„  are  constants. 

Thus  Kantorovich  considers  the  solution  of  an  approximate  (or 
projected)  equation  in  a  convenient  "approximation  subspace"  of  the  space  of 
the  exact  solution. 

1.4  Summary  of  the  Paper 

It  has  been  noted  that  the  pointwise  convergence  of  the  approximate 
operator  basically  depends  on  the  convergence  properties  of  the  approximate 
quadrature  scheme  that  is  applied.   Thus  in  order  to  weaken  the  sufficient 
conditions  (i.e.,  the  hypotheses  of  theorem  1.3.2)  for  the  applicability  of 
the  theory  to  the  solution  of  integral  equations,  we  consider  the  convergence 
of  methods  of  numerical  integration  and  the  convergence  of  the  resulting 
approximate  integral  operators. 
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In  Chapter  2  the  convergence  of  approximate  integration  formulas  is 
considered,  using  the  convergence  properties  of  the  interpolation  polynomials 
and  the  modulus  of  continuity  of  the  integrand.   These  results  are  applied 
to  the  integral  operator  with  a  continuous  kernel.   It  is  shown  that,  for  an 
interpolatory  quadrature  rule,  the  sufficient  conditions  for  the  application 
of  the  uniform  approximation  theory  in  the  space  C  may  be  recast  as  conditions 
on  the  partial  moduli  of  continuity  of  the  kernel  of  the  integral  operator. 
To  extend  the  theory  to  the  case  of  a  mildly  discontinuous  kernel  Anselone  [2] 
carried  out  the  analysis  in  the  space  R  of  proper  Riemann  integrable  functions 
with  the  supremum  norm.   In  extending  the  above  approach  to  the  space  R  we 
obtain  the  more  general  sufficient  conditions  that  the  sets  of  functions 
obtained  by  regarding  the  kernel  as  a  function  of  one  variable  only,  are  both 
regular  sets.   In  the  case  of  a  weakly  singular  kernel,  the  difficulties  with 
the  numerical  quadrature  of  singular  integrands  (Davis  and  Rabinowitz  [ 13]  ) 
preclude  the  application  of  a  direct  quadrature  rule.   However  the  application 
of  a  rule  of  approximate  product  integration  (Atkinson  [ 7] )  ,  regarding  the 
singular  part  of  the  integrand  as  part  of  the  operation  of  integration, 
reduces  this  case  to  one  of  the  above  cases. 

This  approach  is  applied  to  the  case  of  a  singular  integral 
operator  in  Chapter  3.   It  is  shown  that  an  interpolating  polygon  of  a 
function  in  a  Holder  space  converges  to  its  generating  function  only  on  a 
certain  subset  of  the  Holder  space.   Although  the  singular  integral  operator 
is  a  bounded  operator  on  the  Holder  space,  it  is  demonstrated  that  a  sequence 
of  approximate  operators  obtained  by  approximate  product  integration  converges 
pointwise  only  on  an  equi-Holder  continuous  subset  of  the  Holder  space. 
Hence  it  is  concluded  that  the  uniform  approximation  theory  cannot  be  applied 
in  this  case  since  the  hypothesis  of  pointwise  convergence  is  not  satisfied. 
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Gabdulhaev's  direct  method  L 17]  for  solving  singular  integral 
equations  in  a  Holder  space  is  considered  in  section  3.5.   It  is  shown  that 
his  discussion  of  the  convergence  of  the  method  cannot  be  considered  as 
satisfactory  since  the  convergence  obtained  in  the  larger  Holder  space  does 
not  imply  convergence  in  the  solution  space.   Gabdulhaev's  approach  is  based 
on  Kantorovich 's  general  theory  of  approximate  methods.   This  theory  is  con- 
sidered in  section  3.6  for  the  case  of  an  integral  equation.   It  is  shown  tbat 
the  conditions  for  the  application  of  the  convergence  theorem  may  be  reduced 
to  the  same  as  those  obtained  for  Anselone's  theory  in  Chapter  2.   From  this 
analysis  it  is  noted  that  the  implicit  assumptions  made  by  Gabdulhaev  are 
not  necessarily  justified. 

In  Chapter  4  it  is  shown  that  a  uniform  approximation  may  be  ob- 
tained by  considering  the  reduced  equation.   It  is  noted  that  the  hypotheses 
on  the  kernel  need  to  be  modified  slightly  to  obtain  a  consistent  equation. 
The  reduced  equation  equivalent  to  the  singular  equation  may  be  considered  in 
the  space  C.   In  section  4.3  we  define  approximate  operators  by  approximate 
product  integration  of  the  component  singular  operators.   These  approximations 
involve  function  evaluations  only  of  the  given  functions.   It  is  shown  that 
these  approximations  are  collectively  compact  and  converge  pointwise  in  C. 
Hence  the  Anselone  theory  is  applicable  and  it  is  shown  that  the  approximate 
solutions  converge  uniformly  to  the  correct  solution  although  error  bounds 
are  not  easily  computed. 
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2.   THE  UNIFORM  APPROXIMATION  THEORY  FOR  INTEGRAL  EQUATIONS 

2.1  Introduction 

In  order  to  determine  whether  Anselone's  uniform  approximation  theory 
may  be  extended  to  the  case  of  the  singular  integral  equation,  the  basis  for 
the  sufficient  conditions  in  Theorem  1.3.2  is  examined.   It  is  shown  that 
these  conditions  may  be  recast  so  as  to  depend  on  the  properties  of  the  kernel 
of  the  integral  operator.   This  suggests  the  argument  necessary  to  deal  with 
the  operator  with  a  singular  kernel. 

The  abstract  theory  requires  certain  compactness  conditions  on  the 
operators.   This  is  to  compensate  for  the  inadequacy  of  the  weaker  condition 
of  pointwise  convergence  of  the  approximate  operators  obtained  by  quadrature. 
Anselone's  theory  obtains  uniform  convergence  of  the  operator  sequences  by 
restricting  the  class  of  admissible  functions  in  the  space,  C,  of  continuous 
functions  to  the  bounded,  equicontinuous  set  YI3 .      Since  the  operator  K  is 
bounded  by  hypothesis,  Anselone's  theory  in  the  space  C  effectively  requires 
only  that  the  operator  sequence  converges  uniformly  on  an  equicontinuous  set 
of  functions,  and  that  the  set  YJ3   be  equicontinuous.   In  order  to  generalize 
this  result  it  is  necessary  to  consider  what  are  the  equivalent  conditions 
on  other  spaces.   Thus  the  convergence  properties  of  a  sequence  of  approxi- 
mate quadrature  formulas  for  less  smooth  integrands  are  investigated. 

The  uniform  or  supremum  norm  is  also  the  norm  for  the  space,  R,  of 
bounded  Riemann  integrable  functions.   Suppose  that  the  Fredholm  integral 
equation  in  C,  with  a  continuous  kernel  and  a  continuous  given  function,  is 
regarded  as  an  equation  in  the  space  R.   Then  we  may  seek  a  solution  in  the 
space  R  and  such  solution  is  automatically  continuous.   This  follows  since 
for  an  equation  x(s) -Jk(s ,  t)x(t)dt  =  f(s),  if  feC  and  the  kernel  k  is  such 
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that  J|k(s,t)|dt  =  g(s)eC,  then  x£C.   In  particular  equations  whose  kernels 
are  mildly  discontinuous  [  2]  or  weakly  singular  [7]  satisfy  this  condition. 

If  g(s)eR  or  f(s)eR  then  the  equation  is  in  the  space  R  and  the 
solution,  if  it  exists,  is  also  in  R.   Thus  by  obtaining  the  conditions  for 
the  convergence  of  the  approximate  quadrature  formulas  in  the  space  R,  the 
application  of  the  general  theory  is  extended  to  equations  in  the  space  R 
(section  2.4).   It  is  clear  however,  that  we  cannot  infer  by  this  argument 
a  bounded  integrable  solution  if  k  is  a  singular  kernel  since  the  Cauchy 
principal  value  is  not  absolutely  integrable. 

The  results  of  this  investigation  provide  a  basis  for  the  considera- 
tion of  the  convergence  of  approximations  to  singular  operators  in  a  Holder 
space  in  the  following  chapter. 

2.2   Convergence  of  Formulas  for  Approximate  Integration 

It  was  observed  by  Krylov  [27]  that  the  convergence  of  an  inter- 
polator quadrature  rule  is  closely  related  to  the  convergence  of  the 
interpolation  process  itself.   In  particular  if  the  interpolation  process 
converges  uniformly  in  the  range  of  integration  then  the  quadrature  process 
will  also  converge.   Krylov  obtains  regions  of  holomorphy  of  the  function  in 
order  to  obtain  uniform  convergence  of  the  interpolations.   However  we  are 
more  concerned  here  with  obtaining  convergence  results  to  enable  the  identi- 
fication of  those  sets  of  functions  over  which  the  sequence  of  quadrature 
functionals  will  converge  uniformly.   Anselone's  theory  is  basically  con- 
cerned with  this  property  since  the  convergence  of  quadrature  functionals  is 
not  uniform  on  bounded  sets  in  C. 

Consider  a  sequence  of  Lagrange  interpolating  polynomials  on  an 
arbitrary  sequence  of  nodes  in  a  closed  interval.   The  basic  convergence 
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result  for  a  continuous  function  x(t)  is  given  by 
Lemma  2.2.1  (Natanson  [35,  p.  47]) 

If  x(t)ec[a,b]  and  E   is  its  best  approximation  by  polynomials  of 
degree  not  greater  than  n,  then  the  interpolation  polynomial  L  x  satisfies 


L  x-x  ^  L   .E   n 
n   "     n"   n-1 


where  L   is  the  norm  of  the  Lagrange  interpolation  operator. 

n 

The  class  of  functions  for  which  the  sequence  converges  uniformly 
includes  entire  functions  135,  p.  10]  .   But  it  is  substantially  more  restricted 
than  the  class  of  continuous  functions  since  the  Faber-Bernstein  theorem 
[35,  p.  27]  gives  the  inequality  ||l  ||  >  (log  n)/8/rr  in  general. 

However  by  making  the  restriction  to  piecewise  polynomial  inter- 
polation a  more  useful  result  may  be  obtained. 


Lemma  2.2.2 


If  L  denotes  the  Lagrange  interpolation  operator  of  degree  m,  and 


m 


P   denotes  the  piecewise  m-th  degree  polynomial  operator  over  n  sub intervals, 
then 


P     x-xj!    £  mill  ||u)(x,6) 
n  "  v    m"    v        ' 


u  . 


where   uq(x,6)    is   the  modulus   of   continuity  of  xec[a,b]    and  6    =  0(— )    is    the 
norm  of    the   net   of    interpolation   points. 
Proof : 

For  an  arbitrary  value   t   in  one   of   the   n  subintervals 


1      m  1  m 

|P     x(t)-x(t)|    =    |L  x(t)-x(t)|    =    I   S   x(t   )i    (t)-x(t)| 

k=0       k     k 
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where  (t-  tQ)  .  . .  (t-y^t-^)  .  .  .  (t-t,) 

k  ' "  <vv-(vvi)(vw-ivy 

m 
Thus    since     E  jg,     (t)    =   1,    the   right   side   becomes 
k=0  k 


m 
|      2   (x(t  )-x(t))£    (t)| 
i=0         k  k 

tn  , 

u;(x,m5)     2    |i    (t)|  where  6    =  O(-) 

k=0     k  n 


^  m     L       uo  (x,6 ) 

"    m" 


Although   the  Faber   tbeorem   implies    that   a  uniformly   convergent 
interpolation   process    for   every   continuous   function  cannot   be   obtained,    the 
situation    is   more    favorable  with   respect    to   convergence    in   the   mean: 
Lemma   2.2.3    (Natanson  [35,    p.    56]) 

For   an  arbitrary  orthogonal   polynomial    system  with  weight   function 
p(t) ,    the    inequality 


2  2    2  1 

I    p(t){Lnx(t)-x(t)}    dt^Cuj    (x.^TY)J    P(t)dt 


holds    for   every   continuous    function  x(t)  ,   where   C    is    an   absolute   constant. 

These   results   may   then  be   applied   to   obtain   convergence    results    for 

the   associated    interpolatory   quadrature    rules.      Approximating   the    integrand 

by  a   sequence   of    interpolations   of    increasing  order    leads    to   a   sequence   of 

Newton-Cotes    formulas. 

1 
Denote  <f>x   =  J    x(t)dt   and   the   approximate    formula  $   x   =  $L  x   = 

r1  ° 

J    L  x(t)dt.      Thus 
0 


j$x-$    xl    =    U  (x-L  x)  I    £  !|$||    !|x-L  3 
n    '  n      '  n 


16. 
It  is  well  known  (Krylov  [27])  that  this  sequence  does  not  converge 
for  xeC  since  the  coefficients  are  asymptotically  unbounded.   However  from 
Lemma  2.2.1  we  may  deduce  uniform  convergence  over  a  class  of  sufficiently 
smooth  functions.   For  example  (Nathanson  [35])  if  the  nodes  are  the  zeros  of 
the  Chebyshev  polynomial  and  the  function  x  belongs  to  the  Dini-Lipschitz 
class  {x:cju(x,6)  log6-0  as  6-0}  ,  then  the  interpolation  sequence  converges 
uniformly.   Thus  in  this  case,  using  the  Jackson  theorem,  we  have 

|§x-$   xl    £  II $11 1| x-L  x||    <■  C.||$||  .uo(x,-)log- 

1     n  '    Mil.,    n  !  i      ii.ii     n'   °n 

and  the  quadrature  converges  uniformly  over  the  Dini-Lipschitz  class. 

The  interpolation  quadrature  rule  associated  with  a  polygonal  approxi- 
mation is  simply  the  compound  trapezian  rule.   Similarly,  piecewise  quadratic 
interpolations  will  lead  to  a  compound  Simpson's  rule,  and  more  generally, 
piecewise  interpolations  by  polynomials  of  the  m-th  degree  leads  to  a  corres- 


ponding  compound  Newton-Cotes  rule. 
Theorem  2.2.4 

A  sequence  of  composite,  interpolatory  quadrature  functionals  con- 
verges uniformly  on  each  equicontinuous  set. 

Proof:   Consider  a  composite  quadrature  rule  based  on  piecewise  interpolation 
by  polynomials  of  the  m-th  degree.   Then 

|$x-$   xl    =    |$(x-P  mx)|    £  ||$||.||x-P  mxl| 

^  ||$||.m.||L  ||.U)(x,6)  by    lemma   2.2.2 

ii   ii   ,     ii  mil  J 

Since  a  subset,  E  cz  C[0,l]  is  equicontinuous  iff  sup  uj(x,6)-0  as  6-0,  the 

-.  xeE 

result  follows  since  6  =  0(— ). 

n 

The  theorem  applies  equally  to  a  composite  rule  involving  a  weight 
1 
function  p(t).   Let  Qx  =  J  p(t)x(t)dt.   Interpolating  the  operand  x(t)  as 

before  we  obtain 


m 


1 


Q  x  =  QP  x 
n      n 


m 


t) 


p(t)P  x(t)dt 
r    n 
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which  is  the  form  of  approximate  product  integration  used  by  Young  [47]  and 

Atkinson  [7],  where  p(t)  is  supposed  to  be  Lebewque  integrable. 

Corollary  2.2.5 

A  sequence  of  approximate  product  integration  functionals  converges 

uniformly  over  an  equicontinuous  set. 

Proof: 

1 

|Qx-Q  x|  =  |Q(x-P  mx)|  £  J|q||.||x-P  mx||  where  IIqII  -  |  |p(t)|dt. 
'     n  '    '      n1  n  "  n 

The  result  follows  by  an  application  of  lemma  2.2.2  as  in  the  theorem. 

A  particular  case  arises  when  the  nodes  of  the  Lagrange  interpolation 
operator,  L  ,  are  chosen  as  the  roots  of  an  orthogonal  polynomial  system  with  re- 
spect to  the  weight  p(t).   This  gives  the  Gaussian  quadrature  formulas  for  which 
pointwise  convergence  for  each  xeC  is  usually  shown  by  utilizing  the  Weierstrass 
approximation  theorem  (for  example  Natanson  [35],  Davis  and  Rabinowitz  [12]). 
Theorem  2.2.6 

Let  fO  }  denote  a  sequence  of  Gaussian  quadrature  rules.   Then  Q 
n  n 

converges  uniformly  to  the  functional  Q  on  each  equicontinuous  set. 

Proof: 

b 

|Qx-Q  x|  =  If  p(t)[x(t)-L  x(t)}dt| 

1     n  '    '  ■' o  n       ' 

*{J  |p(t)|dt]1/2{J-  |p(t)|  |x(t)-L  x(t)|2dtj1/2 
a  a 

by  the  Holder  inequality 

1     b 
^  C.(ju(x,— jr)  (*   |p(t)|dt    by  lemma  2.2.3. 

a 

'hus  |Qx-Q  x|  ^  C .  I  j  Ql  I  s  up  (A)(x, — -)  ,  and  the  result  follows  as  in  the  previous 

xeE 

:heorem. 

Davis  and  Rabinowitz  [ 12  J  also  show  that  the  convergence  of  the 
luadrature  formulas  is  also  valid  if  x(t)  is  a  bounded  Riemann  integrable 
unction  and  if  the  formula  has  positive  weights.   (Gaussian  rules  and  low 
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order  composite  Newton-Cotes   formulas).      That   is  $nx-§x  as   rr*°   for  each  x6R, 

where  R   is    the   space   of   properly  Riemarm   integrable   functions   on  a  finite 

interval  with  the   sup  norm,    and  $      is   a  positive    linear   functional.      Hence 

to  extend   the   above   results   to   the   space  R,  we   define   the   concept  of   a 

regular   set,   which   gives   a  kind  of  uniform  integrability  condition. 

Definition: 

A  subset  S  C  R   is   regular    iff   for  each  xeS    and  each  m=l,2,..., 

there   exist  x     and  x      in  C   such  that  x  ^x^x     and  $  (x   -x  )-»0  uniformly  for 
m  m  x  nr 

xes   as  nr*30 ,    and,    for  each  fixed  m,    the   sets  S      =  {x   :xes}  ,   S      =  {x   :x£Sj    are 
'  '  mm 

equicontinuous . 

Note   that   this   definition   is   a  weaker  form  of  Anselone's   definition 
of  a  regular   set   in  that   the   sets  S    ,   S      are   required  only   to  be  equicon- 
tinuous  and  not  necessarily  precompact.      However    it   still  follows   that  an 
equicontinuous   set   is   regular  but  not  necessarily  conversely.      (For  example, 
the   set  of   all  characteristic   functions   of   intervals   on  [0,l]).      A  bounded 
set  may  not  be   regular.      The   properties   of   regular   sets   given  by  Anselone  [5] 
are    unchanged  by  this   modified   definition. 
Theorem  2.2.7 

A  sequence  of   positive,    interpolatory  quadrature   functionals   con- 
verges  uniformly  on  each  regular   set  sc  R. 
Proof; 

Since  xeR,    for  each  m=l,2,...,    there   exist  x    ,   x  eC   such  that 

x  %^x     and  $  (x   -x  )-»0  as  nr*0 .      Since   $,   $      are   positive   functionals 
m  m  n 


(x    -x    )    +  $    x    -$x     ^  §    x-$x  ^  §    X    -$x      +  $  (x    -X   ) 
m  nmmn  n  K  m 


If  xes,  a  regular  set,  then  $  (x  -x  wO  uniformly  for  xeS.   Further  an 

m-'  J 
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application  of  theorem  2.2.4  or  theorem  2.2.6,  as  appropriate,  implies  that 

$  -4  uniformly  on  the  equicontinuous  sets  S  ,  S  .   Hence  the  inequality  above 

n  m  n     J 

implies  $  -»$  uniformly  for  xe  S . 

Although  the  above  theorems  are  implicitly  presented  in  terms  of 
real  functions,  the  results  extend  immediately  to  complex  functions  since 
many  assertions  reduce  to  the  corresponding  real  cases  (Anselone  [5]). 

In  practical  applications  it  will  generally  be  more  convenient  to 
use  a  simple  composite  quadrature  rule,  although  the  theorems  of  this  section 
indicate  a  more  general  application.   Thus  the  operator  P  will  refer  to  any 
of  the  polynomial  interpolation  operators  applicable  to  these  theorems.   For 
definiteness  it  will  often  be  considered  as  a  piecewise  polynomial  operator 
of  some  fixed  degree  (often  linear) . 


2.3   Integral  Operators 

These  results  are  now  applied  to  the  case  of  an  integral  operator, 
where  by  an  approximate  quadrature  rule  it  is  to  be  understood  a  rule  of  the 
type  discussed  in  the  preceding  section. 

Consider  initially  the  integral  operator  K,  with  continuous  kernel 
k(s,t),  as  an  operator  from  C  into  C.   Write 

1 
y  =  Kx  where  y(s)  =  J"  k(s  ,  t)x(t)dt . 

0 

We  may  write  k  (t)  =  k(s,t)  =  k  (s)  where  k  ,  k  eC.   We  will  use  the  con- 
vention that  the  subscript  s  will  refer  to  the  first  variable  of  the  kernel 
k  and  the  subscript  t  will  refer  to  the  second.   Thus  Kx(s)  =  $ (k  x)  and 
Knx(s)  =  $  (k  x)  where  $   is  the  approximate  quadrature  rule  as  before. 
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It  can  be  easily  verified  that  the  function  k(s,t),  continuous  as  a 
function  of  two  variables  on  the  closed  unit  square,  implies  that  the  sets 
{k  :0£s^l},  [k  :0^t^l}  are  bounded  and  equicontinuous .   The  converse  result 
is  not  true  in  general  since  a  function  of  two  variables  may  be  continuous 
in  each  variable  without  being  continuous  in  both. 
Definition:  (Timan  [44]) 

The  partial  modulus  of  continuity,  uq  (6)  of  the  function  k(s,t) 
regarded  as  a  function  of  s  is  defined  by 


(6)  =  sup   |k(s+h,t)-k(s,t)  I  =  sup  uo(k  ,6) 
|h[<6 


(JU 

S     (Ks,t<l  t 


The  function  u>  (6)  is  a  modulus  of  continuity  and  the  partial 

modulus  of  continuity,  oj  (6 )  ,  is  defined  similarly.   It  follows  that  if  the 

t 

partial  moduli  of  continuity  exist  as  simple  moduli  of  continuity  then  the 
sets  [k  },  {k  }  are  bounded  and  equicontinuous. 

In  the  previous  section  it  was  noted  that  the  sequence  of  approximate 
quadrature  rules  with  integrand  x  converges  pointwise  if  ||x-P  x||-*0,  and  uni- 
formly on  each  equicontinuous  subset. 
Lemma  2.3.1 

The  sequence  [k  },  of  approximate  operators  obtained  by  approximate 
quadrature,  converges  pointwise  to  the  integral  operator  K  if  the  partial 
modulus  of  continuity  of  the  kernel,  oo  (6 )  ,  exists  as  a  simple  modulus  of 
continuity. 
Proof: 

By  lemma  2.2.2  we  have 
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k  x-P  k  x     <;  Coi(k  x,6) 
s        n  s   "  v   s 


£  c[||x||co(k    ,6)   +  ||k  ||u)(x,6)] 


^  c[||x|'|    sup  uj(k    ,6)   +  sup||k  ||(ju(x,6)] 

S  S 


=  G||x||u)  (6)   +  C.M.ua(x,6)  (2.3.1) 


t 


where  C  is  a  constant  depending  on  the  fixed  degree  of  the  interpolating 
polynomial  and  M  is  the  absolute  bound  of  the  kernel  k(s,t).   Hence,  arguing 
as  in  theorem  2.2.4,  for  the  integral  operator  K, 

|Kx(s)-K  x(s)|    =    I  ($-$    )(k  x)|    =   |${k  x-P  k  x]  I 
£  ||*||  .||ksx-Pnksx|| 
^   C.||§||  .[||x||u)    (6)    +Muu(x,6)] 

by  the  inequality  (2.3.1). 

Thus  for  each  fixed  xeC,  the  right  side  converges  uniformly  with  respect  to  s 

if  {k  j  is  an  equicontinuous  set  (i.e.,  if  uo  (6 )  exists  as  a  simple  modulus 
s  t 

of  continuity).   It  is  concluded  that  ||kx(s)-K  x(s)||-'0  as  n-00  for  each  xeC 

II     v   /    n   v   /ll 

if  iu  (6)-0. 

From  the  proof  of  this  lemma  it  is  also  clear  that  this  operator 
convergence  will  be  uniform  on  each  bounded,  equicontinuous  set,  E,  of 
operands.   While  the  argument  of  theorem  2.2.4  again  implies  that  the  linear 

functionals  $   converge  uniformly  on  the  set  fk  x:0£s^l,  xeE}  ,  whether  or  not 

n       b  J  s 

,  this  set  is  bounded,  it  does  not  follow  in  the  above  proof  that  the  operators 
j  Kn  converge  uniformly  on  the  equicontinuous  set  E,  whether  or  not  this  is 
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bounded.   However  since  bounded  linear  operators  are  being  considered,  a  pre- 

- 

compact  range  will  be  obtained  in  any  case.   It  is  concluded  from  the  lemma 
that  the  equicontinuity  of  the  set  {k  :0^s^l}  is  sufficient  for  the  pointwise 
convergence  of  the  operators  K  to  K. 
Lemma  2.3.2 

The  integral  operator  K  is  compact  if  the  partial  modulus  of  con- 
tinuity of  the  kernel,  id  (6),  exists  as  a  simple  modulus  of  continuity. 

S 

Proof : 

The   operator  K  is   compact    if   {Kx:xe/?}    is   precompact,   where  IS  is   the 
unit   ball.      Thus   for    |s,-s2|<5 

1 
|Kx(s   )-Kx(s    )|    ^J    |k(s      t)-k(s      t)|.Jx(t)|dt 

0 

1 

^  ||x||.J'  u)(k   ,6)dt  £  ||x||.u)   (6). 
0  s 

Hence  {Kx:xe/?}  is  equicontinuous  and  since  the  operator  K  is  clearly  bounded, 
the  result  follows. 

Consider  a  collection  of  subsets  of  the  space  C. 
Definition:   The  collection  of  sets  of  functions  [E  }  is  collectively  equi- 
continuous iff  sup  sup  cu(x,8)-*0  as  6-*0. 
n  xeE 
Definition:   The  collection  of  sets  of  functions  [E  }  is  called  collectively 

precompact  if  [E  }  is  uniformly  bounded  and  collectively  equicontinuous. 
This  follows  from  the  Arzela-Ascoli  lemma  and  it  follows  analogously, 
that  the  set  of  operators  [K  :n^l}  is  collectively  compact  iff  the  collection 
[K /5:n^l}  is  collectively  precompact.   This  gives  the  same  definition  of  a 
collectively  compact  set  of  operators  in  [C]  as  the  more  general  definition 
given  in  Section  1.3.1. 


23 


Lemma  2.3.3 

Let  {k  }  be  the  sequence  of  operators  obtained  from  the  integral 
operator  K  by  a  sequence  of  approximate  quadrature  rules.   Then  [k  }  is 
collectively  compact  if  the  partial  modulus  of  continuity  of  the  kernel, 
U)  (6),  exists  as  a  simple  modulus  of  continuity. 
Proof: 

Since  the  sequence  of  approximate  quadrature  rules  considered  in 

section  2.2  is  pointwise  convergent,  the  Banach-Steinhaus  Theorem  [25,  p.  252] 

implies  that  the  sequence  of  operators  is  uniformly  bounded.   Thus 

11$   II    £  B*^00   for   some   B   and   all  n.      Thus      K  xll    =   sup|$    (k  x)  I    <■ 
ii    nii  ii    n  ii  sr|nvs/| 

||$   ||  .||x||  sup||k  ||    £   B.M.||x||  .      Hence   {K./S'inSl}     is   uniformly  bounded, 
s 
Note    that   a   definite    value    for    the   bound   B   may  be   obtained  directly 

if   required.      In   the    case   of   a   composite    rule    (cf.    corollary  2.2.5)    it    is 

easy    to   show   that   B  ^     Q    .    L  Ii    where   L      is    the    fixed    interpolation  operator   of 

J  iixii     ii     m  |  m  r  r 

degree  m.   In  the  case  of  a  Gaussian  rule,  applying  the  Holder  inequality  as 
in  Theorem  2.2.6  and  noting  the  orthogonality  of  the  basis  polynomials  1,  (t), 
shows  that  B  ^  ||q||  . 

For  two  arbitrary  points  s.,  s„  such  that  |s.-s„|<6, 


Knx(sl)"Knx(s2)l    =    l*n(ks  X"ks  x)l    *  ||$nl|  .||x||  sup|  k(sr  t) -k(s2  ,  t)  | 

12  s 

£   B.||x||  sup   0)(k    ,6)    <>  B.||x||uj   (6). 
t  t  s 


Hence   {k  x:||x||^1,    n^l]     is    collectively   equicontinuous    if   ou  (6)-'0   as  6-*0 . 

n      ll     ii  J         i  s 

Thus   {K/?:n^l}     is    collectively   precompact   and  {l<    :  n^  l}     is    a   collectively 
compact   set   of   operators. 
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From  the  previous  two  lemmas  the  same  condition  on  the  partial 
modulus  of  continuity  of  the  kernel  implies  both  compactness  of  the  operator 
K  and  collective  compactness  of  the  set  of  approximate  operators  {k  }  .   Thus 
in  summary  we  obtain  the  theorem 
Theorem  2.3.4 

If  K  is  an  integral  operator  on  C  with  a  continuous  kernel  k(s,t), 
and  {k  }  is  a  sequence  of  approximating  operators  obtained  by  applying  a 
sequence  of  approximate  quadrature  rules,  then  it  is  sufficient  for  the 
application  of  the  collectively  compact  operator  theory  that  the  two  partial 
moduli  of  continuity  of  the  kernel  exist  as  simply  moduli  of  continuity. 
Proof: 

1.  ao  (6)-»0  as  6-*0  implies  that  K  is  compact  (lemma  2.3.2)  and 
[k  :n^l}  is  collectively  compact  (lemma  2.3.3) 

2.  uo  (6)->0  as  6-*0  implies  that  K  ->K  pointwise  as  n-*00  (lemma  2.3.1) 
The  result  follows  immediately  by  theorem  1.3.2.   Since  k  is  (uniformly)  con- 
tinuous on  the  unit  square  its  modulus  of  continuity,  uo(k;6  ,6„)->0  as 
61'62~*°*   since  max{  uj  (6,),oj  (6„)}  ^  (jo(k;6  ..  ,62)  ,  the  theorem  follows  for  a 
uniformly  continuous  kernel. 

It  should  be  noted  that  this  theorem  for  bounded,  integral  operators 
requires  that  both  the  sets  {k  rO^s^l],  {k  -O^t^l]  be  equicontinuous .   It  was 

S  L 

observed  that  if  k(s,t)  is  uniformly  continuous  then  this  condition  is  satis- 
fied.   It  might  be  suspected  that  in  principle  it  would  be  sufficient  to 
require  k  to  be  continuous  in  each  variable  separately  without  being  con- 
tinuous in  both.   However  from  Timan  [44]  we  have  the  relation 


max{UOs(61),CJUt(62)}  £  oj(k;61,62)  *  oug  (6  L)  +  uut(62) 
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Thus    if  both   sets    are    equicontinuous    then  cu    (6    )—0   and  tx>   (69)-*0   and  hence 
uo(k;6 ,  ,6  9)~*0   as   6    ,69-*0.      That    is,    k    is   uniformly   continuous    as   a   function 
of    two   variables. 

2.4     Equations    in   the   Space   R 

The    last   theorem  does    not   reveal   a   new   result  but    it   does    shift 
the   emphasis    for   sufficient   conditions    onto   the    properties    of    the   kernel   k. 
In   this   context  we   may   investigate   sufficient   conditions   under  weaker   con- 
ditions  on   the   kernel.      Thus    an   extension   to   the    results    obtained   by 
Anselone   [2]    for   the    case   of   a  mildly   discontinuous   kernel    is    suggested.      The 
analysis,    in   this    case,    is    carried  out    in    the    space  R. 

It    is    shown  [2]     that    if    the   kernel  k(s,t)    is   uniformly   t-integrable 
then  KR  c   C,    K.    is    compact    (in   the    space   C)    and   { K.  }     is    collectively   compact. 
Thus   all    the    conditions    for    the   uniform  approximation   theory   to   be    applied 
are   satisfied. 

These    conditions    can  be   weakened    to    require   only   that   K  maps   R 
into  R.      A  parallel   approach   to   that   used    in   the    last   section   for   a   con- 
tinuous  kernel   suggests    that   conditions   might  be   obtained  which   are    similar 
for  both   variables    s,t    (compare    theorem  2.3.4). 

By   theorem  2.2.7   a   sequence   of    positive,    interpolatory   quadrature 
functionals    converges    uniformly   on  each    regular   subset. 
Lemma   2.4.1 

The    sequence   [  K  j    of   approximate   operators    obtained   by   applying 
positive,    approximate    quadrature    rules    converges    pointwise    to   the    integral 
operator  K    if   {k   :  O^s^l}    is    a   regular   set. 
Proof: 


If  {k  }    is   regular   then  fk  xrO^s^l]    is   a  regular  set   for  each   fixed 

o  S 


xeR.      Thus 


26     ; 


K  x-Kx||    =   sup|  ($    -§)(k  x)|~>0   as   n-°° 


by   theorem  2.2.7. 

If    the    analogy  with    the   continuous   case   were    to   be    continued  we 
would  require  {k  }    regular   to  obtain  a  condition  analogous   to    the   compactness 
of    the    operator  K.      This    symmetry    is    implicitly   required   to   obtain   the 
necessary   convergence    condition  ||  (K  -K)K||-*0   since 


2  r*1! 

||KnK-K  ||    =   supj     |kn2(s,t)-k2(s,t)|dt 

s   0 

1 

where   k„(s,t)    =  J  k(s ,u)k(u, t)du  =  i (k  k  ) 

n 

and  k   „(s,t)    =     E     w    .k(s,u    .)k(u    .,t)    =  §    (k  k  ) 
n2v    '    y         .    ,      ni   v        nr    v   nx     '  nv   s   t 

i=l 


Thus  ||  (K  -K)KJ|    =   supj    |  ($    -$)(k  k  )|dt.      The    integrals   converge   provided  the 
s    0 

integrands    are    uniformly   convergent.      By   theorem  2.2.7    this    follows    if 


{k  } ,{k  }    are  both  regular   since   the   pointwise   product   of   regular   sets    is 
regular.      Since   the   convergence    is   uniform  on   regular   sets  we   require   a  con- 
dition to  obtain  a  regular  operator. 
Definition: 

An  operator  K  on  R   is   defined   to  be   regular    if   the   set  {kx:||x||^1} 
is   bounded   and   regular.      That    is   K  maps   each  bounded   subset   of   R   into   a 
bounded  and  regular   set. 

It   is   shown  [2,   Theorem  3.1]     that  K  maps  each  bounded   subset   of  R 
into  a  regular   set    if   the   kernel   is    (properly)   Riemann   integrable  with  re- 
spect  to    (s,t)    and   also  with   respect   to   t   for   all   s.      By   the   above   definition 
such  an  operator   is   regular.      However   this  hypothesis   may  be  weakened 
slightly. 
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Lemma   2.4.2 

Let   k(s,t)    be    (properly)    Riemann   integrable  with   respect   to    t   for 
all   s    and   also  with   respect   to   s    for   all   t.      Then   if  K    is    the    integral 
operator   with  kernel   k,    KRc  R,    K    is   bounded   and    if    the    set   [kfc  (s) :  <*t*l]    is 
regular    then    the   operator   K    is    regular. 
Proof; 

The   hypotheses   on    the   kernel  k(s,t)    imply   the    existence   of    the 
repeated    integrals   [22,    p.    516]     and    thus   KRC   r.      These    conditions    also 
imply   that   k(s,t)    is    bounded    and    thus 

1 
I!k||    =   supj    |k(s,t)|dt   ^  M  where  M   =      sup        |k(s,t)| 
0  0^s,t^l 

Hence   K    is    a   bounded   operator. 

From  the   definition  of   a   regular   set,    {Kx:l|x|M     is   regular    if   for 
each   m=l,2,...    there   exist   continuous    functions    (Kx)m,    (Kx^   such    that 
(Kx)m  *   Kx  ^    (Kx)m  and  *[  (Kx)m-  (Kx^-O   as   nr~,    uniformly   fl   xe/?,    and    the 
sets   {(Kx)m:xe/?},    {  (Kx^rxcflJ    are   equicontinuous    for   each   m. 

Since   {kt}     is    regular    there   exist   continuous    functions    k  m(s) , 
kmt(s)    such    that   kmt(s)    *   kt(s)   *   k™(s),    •(Ik-.fcjH)   as   —  ,    uniformly 
with   respect    to    t,    and    the    sets   {kj,    {  k*}    are   equicontinuous    for   each  m. 
Since   x(t)    is    integrable,    then   so   are    the   non-negative    functions 
*    (t),    x    (t)    where   x(t)    =  x    (t)-x'(t).      Thus 

km(s,t)x+(t)    ±   k(s,t)x+(t)    <.   km(s,t)x+(t), 

kn(».Ox"(t)    ^   k(s,t)x-(t)    <.  km(s,t)x"(t)  and   so 

,  m  + 

(s,t)x    (t)    -    km(s,t)x    (t)   £    k(s,t)(x"t"(t)-x"(t))    =   k(s,t)x(t)    (2.4.1) 
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Hence  defining  the  function 

(Kx)m(s)  -J  kra(s,t)x'r(t)-km(s,t)x-(t)dt 


+  Ct^-k  fs.tU~(t)dt  (2.4.2) 


0 


and    (Kx)    (s)    similarly,  we   obtain  for  each   integrable  xe/?,    the   continuous 

ITi 

functions    (Kx)™,    (Kx)m  such  that 

(Kx)m(s)    S    (Kx)(s)    fe    (Kx)m(s)  (2-4.3) 

Since   the   sets  [k     }  ,   fk™}    are  equicontinuous   for  each  m,    the   sets 

{r1km(s,t)x+(t)dt:0^x+^l}    and  [J   k   (s  ,t)x"  (t)dt:0^x"^l}    are   also  equicon- 

0 
tLuous  by   the   argument  of    lemma  2.3.2.      Since   the   sum  of   equicontinuous   sets   . 

is   equicontinuous,    it    is    inferred   from   (2.4.2)    that   the   sets  {  <Kx)n:||x||£l} 

are  equicontinuous   for  each  m.      Similarly  {  <»c)    :||*M    are  equicontinuous 


' 


for  each  m.     Finally 

*[  (Kx)m-(IOc) J    =  J1{J(km(s,t)x+(t)-km(s,t)x-(t))-(kin(sJt)x+(t)-kin(s,t)x-(t))dts 
0     0 

11  .      , 

=  J   {J(km(s,t)-kin(s,t))|x(t)|dt}ds 

0      0 

By  the  definition  of  the  integral  operator  (2.4.2),  it  is  implicitly  assumed 
(without  loss  of  generality)  that  the  functions  km(s,t)  and  km(s,t)  are 
(properly)  Riemann  integrable  with  respect  to  t  for  all  s.   These  functions 
are  defined  to  be  continuous  (and  hence  Riemann  integrable)  functions  of  s 
for  all  t.   Since  xeR,  these  conditions  imply  the  existence  and  equality  of 
the  two  repeated  integrals  as  before,  and  the  above  integral  may  be  rewritten 
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1   1 

J  {J'(km(s,t)-k  (s,t))ds}|x(t)|dt  (2.4.4) 

0   0 

Since  {k  }  is  regular  the  inner  integral  converges  to  zero  as  m-*00 ,  uniformly 
with  respect  to  t.   Hence  $[  (Kx)  -(Kx)  1-0  as  m-*> ,  uniformly  for  xe/9,  and  thus 
K  is  a  regular  operator. 

It  is  noted  that,  corresponding  to  the  fact  that  continuity  with 
respect  to  s  and  t  separately  does  not  imply  continuity  with  respect  to 
(s,t),  the  hypotheses  of  integrability  with  respect  to  s  and  t  separately 
do  not  imply  the  existence  of  the  double  integral  (compare  for  example 
[22,  example  4,  p.  519])  . 

The  analogous  condition  to  a  collectively  equicontinuous  family  of 
subsets  in  the  space  C  would  be  a  collectively  regular  family  of  subsets  in 
the  space  R. 
Definition: 

The  family  of  sets  {e  }  in  the  space  R  is  collectively  regular  iff 

the  set  (U  E  ]  is  regular.   That  is,  for  each  x  eE   and  each  m=l,2,...  there 
n  n       °  n  n 

exist  x   ,x  ec  such  that  x  ^x  ^x  ,  $  (x  -x   )-0  as  nH30 ,  uniformly  for  x  eE 
nm  n  nm  n  n    v  n  nnr  n  n 

and  uniformly  in  n,  and  for  each  fixed  m,  the  sets  {e  }={x:xeE  ,n^l} 

and  {e   }  =  [x   :x  eE  .n^l]  are  collectively  equicontinuous. 
nmJ     nm  n  n'  j      ^ 

Definition: 

A  set  { K.  :nsl]  of  linear  operators  on  the  space  R  is  collectively 
regular  iff  the  family  of  sets,  {k  x:||x||^l,n^l}  ,  is  uniformly  bounded  and 
collectively  regular. 
Lemma  2.4.3 

Let  the  integral  operator  K  satisfy  the  hypotheses  of  lemma  2.4.2. 
Let  {k  }  be  the  sequence  of  approximate  operators  obtained  from  the 
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operator  K  by  a  sequence  of  positive  quadrature  rules.   Then  the  set 

{k  :n^l]  is  collectively  regular. 

Proof; 

Since  K  x(s)    =  $    (k  x)    =  $ (P  k  x) ,    the  hypotheses   on  the  kernel 
nv/  nv   s   '  v  n  s   '  J r 

imply  K  R  c  R  as    in   lemma  2.4.2.      Since  k   is  bounded,    each  K  s[r]  .      The   se- 
r  J      n  n 

quence   of   positive   quadrature   rules    is   pointwise   convergent  on  R  and   thus 
the   sequence  {||$   ||}    is   uniformly  bounded  by  the  Banach-Steinhaus   theorem  as    in 

lemma   2.3.3.      It    is   concluded     K        =  sup     $    (k  x)    ^B.M  and  {k  }     is   uniformly 

"    n"      1 1    1 1 ;.  -I      "     s   '  '  n  J 

||x||sl 

bounded.  O^s^l 

Applying   a  positive   quadrature   rule    to  obtain  the   approximate 

operator,   K   ,    define    the   continuous    functions    (compare    (2.4.2))    (K  x)    (s)    = 

$    (k  x   -k     x   ).      Since  $      is   positive,    (K  x)    (s)  ^    (K  x)  (s)   by   the    inequality 
nvs  ms      '  n  r  vn/v/         vn/v/y  n  J 

(2.4.1).      (K  x)    (s)    is   defined   similarly. 

For   two  arbitrary  points   s,,s_    such   that     sn-s„     <  6 

11  1     z 

I ,     ,,  m     +     m     +   I         II  i    ||    ||    -hi  .,  m  „ .  , 

$    (k     x   -k     x  )     <     $       .   x  I  sup  au(k,_,6)  .      Thus 
in\s  s        /i         1 1    ni  i    1 1      "tr      vt'y 

sup  sup   uo($  (k  x  ),6)  ^  B  sup  uo(k  ,6) 
n  QSxTsl    n  s  t         '   Z 


Since  (k  j  is  equicontinuous  for  each  m,  sup  cu(k  ,6)-*0  as  6->0  and  thus  the 

t 
set  {$  (k  x  ) : O^x  ^1]  is  collectively  equicontinuous  for  each  m.   Similarly 


nx   s 


{$    (k     x  ) : O^x  £lj    is   collectively  equicontinuous    for   each  m  since  (k  J    is 
nv  ms      '  .  mt 

equicontinuous    for   each  m.      Hence    the    sets  {  (K  x)    (s)  rxe/S^n^l}  , 

v  n  ymv  ' 

{  (K  x)  (s)  :x.eB,r&l}    are  collectively  equicontinuous  for  each  m. 
n  1 

Finally  $  [  (K^)1"-  (Kx)  J     «  J  {^[  (kV-k^x")  -  (k^-kV)]  }  ds 

=    r   {$    (km-k      )|x|}ds 
Jq       nv   s     ms/ '     ' 

1 

=  $n(ix|r  km(s,t)-km(s,t)ds). 
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1 

Thus    U[  (K  x)m-(K  x) J  I    ^   B.|jx||sup|[    k™(s)-k      (s)dsl    and   the    final   term  -0 
1       xn  v    n      in    '  IJ       tv/      mt  ' 

t      o 

as   nr*03 ,    uniformly    for  xe/9  and  uniformly  with    respect   to   n,    since   {k  }     is 
regular. 

Thus   [k  x :  xe/? ,  n^  l]     is    collectively   regular   and   uniformly  bounded, 
and  by   definition    the    set   of   operators   [k   :n^l]     is    collectively   regular. 

From  these    results   we   obtain   a   theorem   for   the    space   R   analogous 
to   theorem  2.3.4. 
Theorem  2.4.4 

Let   K   be   a   bounded    integral   operator   on  R.      Let   fl<  }    be   a   sequence 
of   approximating   operators   obtained   by   a   convergent   sequence   of   positive 
quadrature    rules.      Then    it    is    sufficient    for    the    application   of    the    con- 
vergence   theorem    (theorem   1.3.2)    that    the    sets   [k    ,0->s^l},    [k    -O^t^l]    be 
regular . 
Proof: 

1.  (k  ]    regular    implies    that   K  -*K   pointwise    (lemma   2.4.1) 

2.  {kjj    regular    implies    that   K    is   a   regular   operator    (lemma   2.4.2) 

and  {k    :n^l]     is    a   set   of   collectively   regular   operators 
n 

(lemma  2.4.3) . 
Since    integral   operators    on   R   converge    uniformly   on   regular   sets,    the    proof 
and   conclusions    of    theorem   1.3.2    follow. 

Thus   we   have   obtained    the   irore    "natural"    conditions    for   the    appli- 
cability  of    the    general    theory    in   the    space   R;    namely   that    the    sets   {k   },[kt) 
are   regular.      The    condition   of   {k   }    regular    implies    that   k(s,t)     is  a    (prop- 
erly)   Riemann    integrable    function   of    s    for   all    t   by    the   definition   of    the 
Riemann   integral.      Similarly   [k   }    regular    implies    that   k    is    a    (properly) 
Riemann    integrable    function   of    t    for   ail   s. 
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The  result  of  Theorem  2.4.4  implies  that,  in  general,  we  may  regard 
the  Fredholm  integral  equation  in  the  space  R,  and  seek  solutions  in  the 
space  R  (compare  section  2.1).   In  the  case  that  the  equation  is  actually  in 
the  space  C,  the  solution  obtained  will  automatically  be  continuous,  and  if, 
in  fact,  the  kernel  is  continuous,  then  the  conditions  [k  },  [k  }  regular 
reduce  to  the  conditions  [k  },  [k  }  equicontinuous  (theorem  2.3.4). 

These  conditions  are  more  general  than  those  obtained  by  Anselone  in 

[2].   The  condition  that  [k  }  be  regular  requires  the  existence  of  continuous 

functions  k  (t) ,  k   (t)  but  does  not  make  any  particular  requirement  on  the 

variable  s,  other  than  integrability .   The  condition  of  uniform  t-integrability 

requires  that  these  functions  also  be  continuous  in  s,  which  leads  to  the 

compactness  condition:   $(  k  -k   )-*0  as  s,-*s0. 

Sl   S2         l     2 
Since  the  regularity  conditions  do  not  require  that  the  kernel 

k(s,t)  be  continuous  in  s ,  or  that  the  discontinuities  of  k  (t)  lie  on  a 

curve  which  is  a  continuous  function  of  s,  the  operator  is  not  necessarily 

compact  Lin  CJ.   For  example  the  function 

,   0,   S<l/2 

k(s,t)  =  I  1,  s^l/2 

is  neither  uniformly  t-integrable  nor  the  kernel  of  a  compact  operator  (in 
Anselone's  terminology)  but  does  satisfy  the  conditions  [k  },  [k  }  regular. 

Thus  this  is  the  kernel  of  a  regular  operator  to  which  the  general  theory 
is  applicable. 

It  is  clear  from  the  definition  that  a  compact  integral  operator 

(in  C)  is  a  regular  operator.   Since  a  uniformly  t-integrable  kernel  implies 

that  the  operator  is  compact  [2,  theorem  3.2],  a  uniformly  t-integrable  kernel 

also  gives  a  regular  operator  (this  also  follows  directly  from  the  definition 
using  a  construction  like  that  in  (2.4.2)). 
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While  the  sufficient  conditions  of  theorem  2.4.4  are  neat  and 

symmetrical  they  can  in  fact  be  weakened.   Although  [k  }  equicontinuous 

implies  that  the  operator  K  is  compact  (lemma  2.3.2),  the  compactness  of  the 

1 
operator  also  follows  from  the  much  weaker  condition  j  |  k(s  ,  t)  -k(s  J  t)  |  dt-*0 

0 
as  s  '-»s ,  uniformly  for  0^s,s'^l  L  2 , 5]  .   Analogously,  although  [k  }  regular 

implies  that  K  is  a  regular  operator  (lemma  2.4.2),  this  property  also 

follows  from  weaker  conditions. 

Lemma  2.4.5 

Let  the  kernel  k(s,t)  be  (properly)  Riemann  integrable  with  respect 

to  s  for  all  t  and  with  respect  to  t  for  all  s.   Since  k  eR,  there  exists 

r  s 

k  ,  k  eC  such  that  k  ^k  ^k  and 
s   ms  ms   s   s 


§  (k  -k   )-*0  as  m-*0  for  each  s.  (2.4.5) 

v  s   ms 


Suppose  that  k  ,  k   satisfy  the  compactness  condition 

s    ms       J  r 

^1 

J  |k  (s1,t)-k  (s_  ,t)  |  dt-*0  as  s.-*S„,  OrSs-.S-^l,  for  each  m.       (2.4.6) 
0 


Then  the  integral  operator  K  is  regular. 

Proof: 

1        , 

As  in  (2.4.2)  define  (Kx)m(s)  =  j  k  (s,t)x  (t)-k  (s,t)x  (t)dt  and 

o 

(K  )    (s)    similarly.      The    inequality    (2.4.3)    follows   as    in    lemma    (2.4.2). 
x,   m 

From  [.4,    proposition  4.1]     the    compactness    condition    implies    that    the    con- 
vergence   in    (2.4.6)    is    uniform  and    the    operators   with   the   kernels   k    ,    k      are 
compact.      Thus    the    functions    (Kx)m(s) ,    (Kx)    (s)    are    (uniformly)    continuous 
and    in   fact    the    sets   {  (Kx)m:||  x||^l}  ,    { (Kx)    :l|x||^l]    are   equicontinuous    for 
' each   m. 
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1 

Since   the   function  J    k(s,t)dt   is   a  continuous   function  of   s   on 

0 
[0,l]    by    (2.4.6)    the   functions   of    (2.4.5)    are   boundedly  convergent.      Hence 

1   1 
$[  (Kx)in-(Kx)J^!|x||  .J  J    (km(s,t)-k   (s,t))dt  ds->0  as  m-°° ,    uniformly  for  xe#. 

0    0 
Thus  K   is   a  regular  operator. 

Definition:    L  10] 

A  sequence  {e  }    of   sets   of   functions    in  R   is   asymptotically 

equicontinuous    if   sup   |  x(s)  -x(t)  |-»0  as   s-*t  and  n-*00 ,   uniformly  for   te[0,l]  . 
xeE 
n 

This    is   a  rather  weaker  condition  than  collectively  equicontinuous 
defined  above.      Similarly  we   may  define   a  weaker  form  than  collectively 
regular. 

Definition: 



A  sequence  {e  }  of  sets  of  functions  in  R  is  asymptotically  regular 

if  for  each  fixed  m,  the  sets  {e  } ,  {e   }  (see  the  definition  of  collectively 

n      nmJ  v  \ 

regular)  are  asymptotically  equicontinuous  and  $  (x  -x  )-»0  as  m-^°  and  n-*° . 
J    r        J  v  n  nnr 


Lemma  2.4.6 


Let  {E  }  be  asymptotically  regular  and  each  E   a  regular  set  in  R. 


n 


Then  [e  }  is  collectively  regular. 
Proof; 

Since  each  E   is  regular  the  sets  E  ,  E   are  equicontinuous  for 

n      b  n'   nm      n 

each  m,n.   Further  $  (x  -x   )-*0  as  tcH* ,   uniformly  with  respect  to  x  eE  for 

v  n     nm  J  n     n 

each  n.      Since   [e  }    is   asymptotically  regular   the   sets   [e   } ,   {e     }    are 
asymptotically  equicontinuous.      Anselone's   result  using   a  construction  based 
on   the   usual  proof  of    the  Arzela-Ascoli   theorem  L5,    proposition  2.4]    implies 
that  U   E      is    equicontinuous.      That   is,    in   the    present    terminology,    {e   }, 


[E      }    are   collectively   equicontinuous    for   each  m. 


nm 
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Further   $  (x    -x      )~»0   as   m-*30   and  n-*° .      Thus,    given  e>0,    there   exists 
v   m     nm  '    b  ' 

M,  ,N   such   that  $  (x    -x      )<e    for  n>N ,    m>M,  .      Since    this    convergence    also  holds 
1  v   m      nm  1  to 

for   each   n,   we    also  have    for   n^N ,   $ (x   -x      )<e    for  m>M(n)  .      Thus   $  (x   -x      )<e 

for  m>M   and   for   all   n  where  M   =  maxCmax  M(n) ,M, ) ,    i.e.,   $  (x   -x      )-*0   as   m  ■>CD , 

i  m      nm 

n^N 

uniformly  with   respect   to   n.      Hence   [E   }     is    collectively    regular   by 
definition . 


Lemma   2.4.7 

Let  the  kernel  k(s,t)  satisfy  the  hypotheses  of  lemma  2.4.5  and 

also  suppose  that  {k  }  is  regular.   Then  (k  :  n^  l]  is  collectively  regular 

where  IK  }  is  the  sequence  of  approximate  operators  obtained  from  K  by  a 

sequence  of  positive  quadrature  rules. 

Proof; 

Defining   the    functions    (K  x)    (s) ,    (K  x)    (s)    analogously   to   those    in 
o  v    n  n     m  °  J 

lemma   2.4.5  we   have 

|  (K  x)m(s)-(K  x)m(s,)|    <.  $    (|km-km,|)||x+||    +  *    (Ik      -k      ,|)l|x"||. 
lvnv/vn/v      'x  nV|s      s|/M      "  nV|    ms      ms    '  " 


Since   {k   }    is    regular    the    sets   {k   }  ,    {k      }    are   equicontinuous ,    and   hence 
s  s  ms 

$    (Ik   -k    ,|)  -<fr(|k   -k    ,|)    as   n~*° .    uniformly    for   0^s,s'^l   by    lemma   2.4.1. 
nlssllssl 

By  (2.4.6),  $  (|  k  -k  ,  |  )~'0  as  s  '-s  ,  uniformly  for  O^s^l  (compare  lemma  2.4.5). 

o  S 

Thus    (K  x)m(s)-(K  x)m(s')-0   as   rr1™    and   s'-s,    uniformly   for  ||x||^l 
and   O^s^l.      This    implies    that    the    sequence   of    sets   {  (K  x)    :||x||^l]    is 

asymptotically   equicontinuous    for   each    fixed   m. 

1 

Further  *[  (K  x)m(s)-(K  x)    (s)]    =   [  ($    (km-k      )|x|)ds.      As    rr*° ,    the 
nnm  «Jxns      ms1 

0 

right   side   converges   uniformly   with    respect    to      s    to 

11 

J  J  (km(s,t)-k  (s,t))|x|dt  ds  by  the  same  argument  as  above.   This  expression 

0  0         m 
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converges  to  zero  as  m-"30  using  (2.4.5)  (compare  lemma  2.4.5). 

Thus    the   sequence   of   sets  {k  x:||x||^l}    is   asymptotically   regular. 

Since   each  of   the   sets   {k  x;    x  £l}    has   a   finite  dimensional   range   it    is 

n  ^j 

regular.      Thus  by   lemma  2.4.6   the   sets    {k  x:||x||^l,n^l}    are   collectively 

n 

regular.        Uniform  boundedness   follows   as    in   lemma  2.4.3   and  hence  {k  :n^l} 
is    collectively  regular. 

In  connection  with   this   proof,    note    that,    as   defined,    the   sets 

m 
{  (K  x)    (s):xe/6>]    are   not  equicontinuous   for  each  m,n.      This  may  be   compared 

to  Anselone's   treatment  of  discontinuous  kernels   [5]    where  KR  C  C  but 

a 
KR<^   C    in  general.      However   since   k   (s)eR,    there   exist  k   (s)  ,    k.    (s)eC   such 

that  k„    <=k  ^k{   and  $(kX-k,    )-0  as  A-«   for  each   t.      Thus   {$    (k^x+)  :  0^x+£l]    is 

it      t   t        t  it  nv  s   ' 

equicontinuous  for  each  i   and  each  n  since  it  involves  a  finite  number  of 

c  i  , 

continuous  functions.   Hence  {  (K  x)  (s):xe/S>j  (defined  in  terms  of  the 

i 
functions  k  ,k„)  is  equicontinuous  for  each  i    and  each  n. 

i' 

Hence  in  order  to  see  that  the  lemma  2.4.6  is  in  fact  applicable, 
we  may  take  for  the  majorizing  and  minorizing  sets  of  functions  (K  x)  , 
(K  x)   where  p=m  for  n>N  and  p=i    for  n^N.   For  p=m  the  collection  is 
asymptotically  equicontinuous.   Thus  given  e>0,  there  exist  6(e),  N(e) 
such  that 


sup     sup       I (K  x)F(s)  -  (K  x)p(s')|  <  e 
n>N(e)  xe/?        '  n  n       ■ 

|  s-s  '  |<6  (e) 


For  n^N(e)  we  have  a  finite  collection  of  equicontinuous  sets  with  p=i 
and  thus 


SUP     sup       | (K  x)P(s)  -  (K  x)P(s')  <  e 
n^N(e)  xe/?         n  n 

| s-s  *  |<6 (t) 


I 
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Combining    these    two  collections   of   sets   of   functions    it    is    clear    that 

{  (K  x)    :xe/?}     is   collectively   equicontinuous    for   each   p. 

11  l  l  I 

For  n^N(e),  J($  (k  -k  )|x|)ds  =  $  [J  k  -k  ds]  |x|  which  converges 

to  zero  as  l-*™   for  each  n  since  the  inner  integral  converges  for  each  t.   For 

n>N  the  expression  converges  as  in  the  lemma  2.4.7  and  hence 

$[  (K  x)  -(K  x)  ]-»0  as  p-"00  uniformly  with  respect  to  n,  as  required  in  the 

proof  of  lemma  2.4.6. 

Collecting  these  results  we  obtain  the  more  general  convergence 
theorem 
Theorem  2.4.8 

Let  the  kernel  k(s,t)  satisfy  k  eR  and  [k  }  regular.   Suppose 

further  that  the  majorizing  and  minor iz ing  functions  k  ,  k  eC  satisfy  the 

condition  J  1 1 (s , t) -I (s ' , t) | dt^O  as  s '-s  ,  0^s,s'^l,  for  each  m. 
0 
Then  K  ~*K,  K  is  regular  and  {k  }  is  collectively  regular.   Thus 

the  convergence  theorem  (theorem  1.3.2)  follows. 

Proof: 

Since  {k  }  is  regular,  pointwise  convergence  of  the  approximate 

operators  follows  by  lemma  2.4.1.   K  regular  and  {k  }  collectively  regular 

follow  from  lemmas  2.4.5  and  2.4.7  respectively.   Since  integral  operators 

on  R  converge  uniformly  on  regular  sets,  the  proof  and  conclusions  of  theorem 

1.3.2  follow  as  in  theorem  2.4.4. 


2.5  Quadrature  for  Unbounded  Functions 

Since  we  are  ultimately  concerned  with  kernels  with  an  unbounded 
singularity,  a  next  step  in  progressively  relaxing  the  conditions  on  the 
kernel  is  to  consider  the  case  where  the  kernel  is  absolutely  integrable. 
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This  leads  to  an  investigation  of  numerical  quadrature  for  integrable 
singular  functions.   As  in  section  1.2,  this  problem  may  be  considered  either 
by  approximate  product  integration  or  by  a  direct  method. 

The  direct  method  is  the  process  of  "ignoring  the  singularity"  which 
was  studied  by  Davis  and  Rabinowitz  Ll3,40]  .   That  is,  if  x(t)  becomes 
singular  at  a  point  §,  then  one  defines  x(§)  =  0  (or  any  other  finite  value) 
and  then  approximates  the  integral  by  a  usual  numerical  quadrature  rule.   It 
is  shown  that  the  process  may  or  may  not  be  valid  depending  on  the  nature 
and  position  of  the  singularity.   Thus  if  the  singular  point,  § ,  is  a 
"rational"  point  in  the  range  of  integration  then  compound  quadrature  rules 
do  approximate  the  integral  when  the  integrand  is  monotone  near  £ •   On  the 

r- 

other  hand  if  x(t)  is  monotonic  but  £  is  an  irrational  point  then  L 13 , 
theorem  5]  the  process  may  not  be  theoretically  valid  and  in  general  only 
"lim  inf  convergence"  can  be  asserted.   These  results  were  generalized  by 
Miller  [32]  who  replaced  the  assumption  of  monotonicity  by  the  more  general 
condition  that  x(t)  can  be  dominated  by  a  monotone,  integrable  function. 
However  the  restriction  to  rational  singularities  is  retained. 

Consider  the  weakly  singular  operator.   The  integrand  in  this  case 
t-s |  ,  0<Q<1.   Regarding  s  as  temporarily  fixed, 
the  above  remarks  imply  that  convergence  for  an  irrational  value  of  s  cannot 
be  asserted  for  a  sequence  of  approximations  obtained  by  a  direct  method. 
Thus  if  Kn  denotes  the  approximate  operator  obtained  by  a  quadrature  of  this 
type,  we  certainly  cannot  assert  K  x(s)-»Kx(s)  uniformly  in  s.   That  is,  we 
cannot  assert  pointwise  convergence,  and  the  first  hypothesis  of  the  con- 
vergence theorem  1.3.2  is  not  satisfied. 

However,  the  method  of  approximate  product  integration  neatly 
solves  this  problem.   The  singular  part  of  the  integrand  is  absorbed  into 
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the  integration  operator.   Let  the  kernel  k(s,t)  =  p(s  ,t)i  (s  ,t),  where  p(s,t) 
(p(s,t)  =  J  t-s I   ,  0<Q'<1)  is  absolutely  integrable  and  £(s,t)  is  continuous. 
Lemma  2.5.1 

The  sequence  [K  }  of  approximate  operators  obtained  by  approximate 
product  integration  converges  pointwise  to  the  weakly  singular  operator,  K, 
if  the  partial  modulus  of  continuity,  m  (6 ) ,  of  the  continuous  part,  £(s,t), 
of  the  kernel  exists  as  a  simple  modulus  of  continuity. 
Proof: 

Using  corollary  2.2.5 

||k  x-Kx||    =   sup|$[p  fp  I    x-£   x}]  |    <■   C    sup{|$(p   )|u)(i    x,6)] 

n  e  S       II    a  s  S  S  S 


^   C.||Q||[||x||u)t(6)    +  L   iu(x,6)]  ,    by    (2.3.1), 


where   C    is    a   constant   depending   on   the    fixed   degree    of    the    interpolating 

polynomial,   Q    is    the   weakly   singular    integral   operator  with  kernel   p,    and  L 

is    the    absolute   bound   for  4(s,t).      If   lz   }    is    equicontinuous  we   have   con- 

s 

vergence  for  each  fixed  xeC. 

Note  that  in  practical  cases  part  or  all  of  the  function  i.(s,t) 
may  be  absorbed  into  the  integral  operator  in  order  to  obtain  better  con- 
vergence properties  of  the  polynomial  approximations.   Also  the  Gauss  type 
formulas  may  be  applied  using  Theorem  2.2.6. 
Lemma  2.5.2 

The  integral  operator  K  is  compact  if  {&    }    is  equicontinuous  and 

1 
J'  |p(s,t)-p(s',t)|dt-0  as  s'-s,    0^s,s!£l  (2.5.1) 

0 


i 
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Proof: 

The  convergence  in  (2.5.1)  is  uniform  and  max  $(|p  |)  exists 

s      s 

[4,  Theorem  4,13  .   Thus  K  is  bounded  and  1 1 k| | £L . max  $(|p  |)  =  L.||q||.   For 


|srs2|<6, 

|kx(Si)-Kx(s2)|  <■   J  |p(s1,t)A(s1,t)-p(s2,t)A(s2,t)|.|x(t)|dt 

1 
£  ||x||if(|p(s1,t)|u)at,6)+|p(s1,t)-p(s2,t)|.U(s2,t)|)dt 

1 
*  |W|((0  (6).||q||  +  L.J  |p(S]L,t)-p(s2,t)|dt) 


Thus  K9  is  bounded  and  equicontinuous  using  (2.5.1)  and  [jl   }    equicontinuous . 

The  lemma  was  expressed  in  this  form  to  show  the  analogy  with  the 
previous  lemmas,  but  it  is  clear  that  if  [l   }  is  not  equicontinuous,  the 
operator  will  still  be  compact  if  k  satisfies  (2.5.1). 
Lemma  2.5.3 

The  sequence  {k  }  of  approximate  operators  is  collectively  compact 
if  (2.5.1)  holds  and  £(s,t)  is  continuous. 
Proof: 

By  corollary  2.2.5,  the  sequence  {Q  }  of  approximate  quadrature 
rules  is  pointwise  convergent  and  the  Banach-Steinhaus  theorem  implies  that 
the  sequence  is  uniformly  bounded  as  in  lemma  2.3.3.   Thus 


Ik  xll    =   suplo    (l   x)  I    £  ||o  II  .||x||  .suplU   II    &  B.L.||x| 
1    n  "  o        nv   s    /j        l|Xn"    "    "      „r|1    s"  "    ' 


Hence  lK/9:n^l}  is  uniformly  bounded. 
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For   two   arbitrary   points    s.,s      such   that    | s    -s    |<6 , 

1 
(Kx)(Sl)-(Knx)(s2)|    =    |J,p(s1,t)Pn(^(s1,t)x(t))-p(s2,t)Pna(s2,t)x(t))dt| 

i 
^  ||x||[f    |p(s1,t)-p(s2,t)|.Pn(U(s1,t)|)dt 
0 

1 
+  J|p(s2,t)|Pn(|i(Sl,t)-^(s2,t)|)dt] 


Since   [i   }    is   equicontinuous    the    integrals    converge   uniformly  with   respect 
to   s    (by    lemma   2.5.1)    to 


1  1 

[{  |p(Sl,t)-P(s2,t)|  .|jj(Sl,t)|dt  +  f  |P(s2,t)|  .|j&(Sl,t)-ja(s2,t)|dt] 


This    term  converges    to   zero   as   s.-*s_    by    lemma   2.5.2    and   hence 

(K  x)(s1)-(K  x)(s9)-*0   as   n—      and   s..-*s9,    uniformly    for   ||x||^    1   and  O^s   ^1. 

That   is   { K  13}     is    asymptotically   equicontinuous.      Since   each  K     has    a   finite 

n  J    r  J      n  n 

dimensional   range    it    is   equicontinuous,    and  hence   tR/S^nSl}    is    collectively 
equicontinuous . 

It   follows    that  {k    :n2:l}    is    collectively   compact. 

Note    that    if    the    approximate    product    integretion   rule    is    obtained 
by   piecewise    polynomial    interpolation,    the    collective   equicontinuity   is 

easily  obtained   by 

1 

|Knx(Sl)-Knx(s2)|    =    |,fj)(s1,t)P^(s1,t)x(t)-p(s2,t)P^(s2,t)x(t)dt| 

^    f1|p(s1,t)-p(s      t)||p*Vs      t)x(t)|dt 
0  L  1 

+   |'    |p(s2,t)|  .|p™(A(s1>t)-£(s2,t))x(t)|dt 

1° 

<  ||x||[l.||lJ|J  |p(s  ,t)-P(s2,t)|dt  +  ||q||.||lJ|o)s(6)] 

0 


42 


which   gives    the    result   as    in   lemma   2.5.2. 

Hence  we   obtain  Atkinson's   result  [7]    for  weakly   singular  kernels 

in   terms    of    the   current   development. 

Theorem  2.5.4 

Let  K  be  an  integral  operator  such  that^the  kernel  k(s,t)  - 

p(.,t)X(.,t)  satisfies  p.,^(0,l)  **   <**«  «*  J  [«*VWV»\««>  » 

s   -s  ,  OSS  ,s  *l,      If  fKn]  is  a  sequence  of  approximating  operators  obtained 
hy  allying  ^sequence  of  approximate  product  Integration  rules,  then  It  is 
sufficient  for  the  application  of  the  collectively  compact  operator  theory 
that  the  two  partial  moduli  of  continuity  of  *(.,t)  exist  as  simple  moduli 

of  continuity. 

Proof: 

The  theorem  follows  from  lemmas  2.5.1,  2.5.2,  2.5.3  as 

Theorem  2.3.4. 

The  remarks  in  this  section  suggest  that  the  approximation  of  a 
strongly  singular  Integral  operator  is  probably  best  achieved  by  some 
method  of  approximate  product  integration. 
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3.   APPLICATION  TO  SINGULAR  INTEGRAL  EQUATIONS 

3.1  Preliminary  Remarks 

It  has  been  shown  in  the  preceding  chapter  that  the  conditions  for 
the  application  of  the  abstract  convergence  theorem  are  the  pointwise  con- 
vergence of  the  approximate  operators  and  the  identification  of  those  sets  over 
which  this  convergence  is  uniform.   These  conditions  were  ultimately  reduced 
to  the  requirement  that  the  kernel  of  the  integral  operator  have  suitable 
properties . 

The  application  is  relatively  straightforward  for  continuous  kernels. 
In  the  case  of  a  weakly  singular  kernel  the  convergence  properties  of  approxi- 
mate quadrature  formulas  for  improper  integrals  cause  problems.   These  are 
avoided  by  approximation  with  an  approximate  product  integration  rule  (section 
2.5).  However,  this  technique  cannot  be  immediately  extended  to  singular 
equations  with  a  Cauchy  kernel  since  this  kernel  is  not  absolutely  integrable. 
"urther  the  integral  must  be  interpreted  as  a  Cauchy  principal  value  and  hence 
conditions  for  the  existence  of  this  integral  as  well  as  the  convergence  of  a 
sequence  of  principal  value  integrals  must  be  considered.   Stewart  [42]  has 
)bserved  that  even  if  a  sequence  of  functions  is  uniformly  convergent  in  C,  and 
:he  Cauchy  integrals  of  all  functions  in  the  sequence  exist,  the  sequence  of 
lauchy  integrals  may  not  converge. 

Hence  the  properties  of  the  singular  integral  operator  are  considered 
n  order  to  determine  conditions  analogous  to  those  in  Chapter  2,  for  the 
lointwise  and  uniform  convergence  of  a  sequence  of  approximations  to  such 
perators . 

It  is  assumed  that  the  definition  and  properties  of  a  function  satis- 
'  i-ng  a  Holder  condition  with  index  u,  (denoted  cpeH(p,))  are  known  (for  example 
133,  chapter  1]  )  .   Functions  which  satisfy  the  Holder  condition  with  larger 
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exponents  are  "more  benign"  than  those  which  satisfy  the  condition  with  a 
smaller  exponent  (that  is  cv<(3  implies  H(P)  C  R(a)  ,  [34,  chapter  4]).   The 

Oi 

statements  cpeH(|i)  and  uo(9,S)^M6   are  equivalent  for  an  arbitrary  constant  M. 
The  class  of  functions  satisfying  the  H()J,)  condition  form  a  Banach  space  with 
the  appropriate  norm  (Lorentz  [29,  p.  50]  )  . 


The  Riemann  integral  (with  respect  to  t)  of  the  function  cp  (t)/  (t-s) 
may  not  exist  even  as  an  improper  integral  if  cpeC,  but  (Mikhlin  L  31 ,  p.  115]) 
it  does  exist  in  the  sense  of  the  Cauchy  principal  value  if  cpeH(|J,)>  CKu.^1. 
The  idea  of  a  principal  value  is  easily  extended  to  contour  integrals 
(Muskhelishvili  [33]  ,  Pogorzelski  [38])  and  we  note  the  important  Privalov 
theorem: 
Theorem  3.1.1  (The  Privalov  Theorem)[38,  p.  439]  ,  [33,  p.  46] 

If  a  complex  function  cp  (t)  ,  defined  on  a  smooth  arc  L,  satisfies 
the  H((a,)  condition,  then  the  values  of  the  Cauchy  singular  integral, 
$(s)  =  [T  y^dt,  at  points  s  on  the  line  of  integration,  except  in  arbitrary 
small  neighborhoods  of  those  ends  where  cp(t)?fe0,  satisfy  the  H(p,)  condition 
for  |ji<l,  and  the  H(l-s)  condition  for  \i=l,  where  e  is  an  arbitrary  small 
positive  number. 

We  note  also  the  Poincare-Bertrand  transformation  formula  for 
iterated  singular  integrals  [33,  p.  57]  :   If  L  is  a  smooth  arc  or  contour, 
cp(t,u)  a  function  of  the  two  points  t,u  on  this  line,  satisfying  a  Holder 
condition  with  respect  to  t  and  u,  and  let  s  be  a  fixed  point  on  L  not 
coinciding  with  one  of  its  ends,  then 

r  -£  r  ^xi^du  =  -TT2cp(s,s)  +  r   r     ?<fc>u>    dt  du . 

J  L   t-s    JL  u-t  Y  v    '    '        J  L  -JL    (t-s)(u-t) 
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3.2   Properties  of  the  Singular  Integral  Operator 

The  Hilbert  transformation  of  a  function  f(t)  is  defined  by 


g(s)  =  Hf(s)  =  Hs[f(t)]  =^pf  ^-   dt. 


In  the  terminology  of  Fourier  series  this  is  often  called  the  conjugate 
function  to  f (t) .   Since  the  Cauchy  principal  value  integral  is  conditionally 
but  not  absolutely  convergent,  the  existence  of  g(s)  is  not  trivial  even  if 
f(t)  is  continuous.   Zygmund  [49]  notes  (by  means  of  the  conjugate  Fourier 

00         g  ~\r\     TIX  °^         COST1X 

series:      E.  — : ,     Z~  — : )    that   the    function  con  "jugate   to   a  continuous 

n=2   nlogn        n=2   nlogn7 

function  need  not  be    continuous    or  even   bounded.      Thus    for    the  Hilbert 
operator,    H,    it   cannot   even  be    asserted   that   HC  c   R, in   contrast   to  Anselone's 
requirement    that   KR  c   C   for   a  mildly  discontinuous   kernel. 

However   the   Privalov    theorem  shows    that    the    linear   operator  H 
transforms   H(pJ)    into    itself    for   CX(JL<1 .      Mikhlin  [3l]  ,   Mus  khelishvili  [33] 
and   Pogorzelski  [38]     assume    the   known   functions    are   Holder   continuous   and 
derive    solutions    of    the    singular    integral   equation    in   the   space   of   Holder 
continuous    functions.      Thus    the   apparent   natural   analogy   for    the   space   C    in 
the   case   of    the   Fredholm   integral   equations,    is    the   Banach   space  H(|~i)    in   the 
case   of   singular    integral   equations. 

The   function   feH(M-)    will  be    said    to   be   of   exact   order  p,    if    it   does 
not  belong   to   any  smaller   space,    i.e.,    f(t)^H(v),   \>>\x  .      Then  we   have    the 
minor   extension   to   the   Privalov   theorem: 
Lemma  3.2.1 

Let   f(t)eH(M-)    be   of  exact   order  |i  ,    (K+J.<  1 .      Then   the   Hilbert   trans- 
formation of   f    is   also  of   exact  order  \x  . 
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Proof : 

Let  g(s)  =  H  [f(t)].   Suppose  g(s)eH(v),  v>|Jb  .   (This  does  not 
contradict  the  Privalov  theorem  since  H(v)  c:  H(jj,)).   Then 


Hf 


ig(s)]  =  .f(t)  +-ipjLf(u)dU  p.r(s.t^u.s) 

But   as    is   easily  verified 

PJL(s-t)(u-s)    =  u-t  PJJs^t    "  7^ds   ~  ° 

Thus  H  Cg(s)]  =  -f(t),  and  the  Privalov  theorem  implies  that  since  geH(v) 
then  feH(v)  ,  contradicting  the  hypothesis  that  f  is  of  exact  order  p, . 

These  remarks  assume  the  existence  of  functions  belonging  exactly 
to  a  given  Holder  class.   Such  functions  have  been  identified,  for  example 
by  Loud  [30]  who  constructs  (based  on  Knopp ' s  construction  of  a  continuous 
non-dif ferentiable  function)  a  class  of  functions  satisfying  at  every  point  a 
Holder  condition  of  prescribed  order.   Also  by  Hardy  [20]  who  effectively 
obtains  a  similar  result  [20,  p.  303]  using  Weirstrass1  non-dif ferentiable 
function. 

3.3  Convergence  in  a  Holder  Space 

Let  H(j3  )  be  the  class  of  functions  defined  on  an  interval  of  length 

and  satisfying  a  Holder  condition  of  order  $.      Let  the  norm  be  defined 

||xl|g   =  M^x]   +M2[s;0],   xeH(B)   where 

jx(t1)-x(t2)| 
M    [x]    =  max|x(t)  |  ,    M 9[x;g]    =   sup  = . 

2  o<|trt2K     |trt/ 

Then  the  space  H(|3)  is  a  Banach  space  [29,  p.  50]. 
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By  an  obvious   extension  of    the    idea  of   equicontinuity ,   we    define    a 

set,   E,    of   functions    to  be   equi-Holder   continuous. 

Definition: 

The    set,   E,    of    functions    is   equi-Holder   continuous    of   order  [3 

[equi-H(g)]     if   given  e>0,    there   exists   6>0   such   that,    if    |  t   -t    |<6  ,    then 

x(t1)-x(t2)| 
-z <   e    for   all  xeE    and   tn,t„    in   the    interval.      Equivalently ,    in 

lVt/ 

terms  of  the  modulus  of  continuity,  if  sup  — W-2 — ^-  -*0  as  6-»0. 

xcE  6P 

Since  the  abstract  approximation  theory  implies  uniform  convergence 
on  compact  sets,  it  is  useful  to  identify  the  compact  sets  in  H((3).   We  need 
also  to  define  the  class  h(|3)  (apparently  introduced  by  Hardy  [2l]  as  the 
class  Lip  S)  • 
Definition :  L  2 1] 

The  function  x(t)  belongs  to  the  class  h(g)  if  x(t+h)-x(t)  = 

Q 

0(|h|    )    as    h-'O,    uniformly    in    t. 

Clearly   h(|3)    is   a   closed   subspace   of  H  (P )  .      Johnson  [24]     character- 
izes   the    relatively   compact    subsets   of   sets   of    functions   defined   on   a  compact 
interval   by 
Lemma  3.3.1  [24,    p.    154] 

Let  E   be   a   subset   of   h(l).      Then  E    is    relatively   compact    iff  E    is 

I  x ( s) -x ( t) I 
bounded   and      lim       J — ^-r2 r-^-1  =   0  uniformly   on  E. 

ls-tl-0    l3^' 

Corollary  3.3.2 

, X 

If  (Ko<P£l  then  the  unit  ball  of  H(3)  is  compact  in  h(a)  C  H(a). 

Thus  the  relatively  compact  subsets  of  h(cO  are  identified  as  bound- 
ind  equi-H(o')  sets.   These  sets  are  also  compact  in  the  larger  space  H(o0 
(with  the  same  norm).   Since  the  unit  ball  of  H(3)  is  compact  in  the  topology 
af  H(cy),  it  is  a  fortiori  compact  in  any  coarser  topology,  and  in  particular 
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in  the   sup-norm  topology  of  h(0)    (the   space   C) .      Thus   the   Arzela-Ascoli 
theorem  is    reclaimed  as   a  special   case. 

Q 

It  was   noted   that  xeH(B)    iff   the   expression  u)(x,6)/6      is  bounded 
for   all   sufficiently   small  values   of  6 .      We   assert   that  x(t)   has   exact  order 
|3    iff   there   exists   at    least  one   point,    t    ,    such   that  x  has  exact  order  3    at 
that   point.      Thus   modifying  Loud's   definition  [  30]    slightly  we   obtain: 
Definition: 

A  function  x(t)    satisfies   a  Holder   condition  of  exact  order  B    if 
there   exists   positive   constants  K, ,   K„    such   that 

(a)  For   any  t   and  any  sufficiently  small  6>0, | x(t+6 ) -x(t) |<K.6 

(b)  3   t     and  an   infinite   sequence,    {6   }  ,    such   that  6  -»0  as  nr°°   and 
v    '  o  m  m 

|x(t  +6    )-x(t   )|>K06P. 
1     v   o     nr       v   o/  '       2   m 

In  section  2.2  it  was  seen  that  convergence  of  the  approximate 
quadrature  rules  on  C  depends  on  the  convergence  of  the  piecewise  polynomial 
approximations  of  the  integrand  in  C.   Thus  we  are  also  concerned  with  the 
convergence  in  H(B)  of  the  polynomial  approximations  to  a  function  in  a 
Holder  space. 
Lemma  3.3.3 

Let  L  denote  the  Lagrange  interpolation  operator  of  degree  m  on 


m 


,m 


equally  spaced   intervals,    and   let  P     denote   the    piecewise,   m-th  degree   poly- 
nomial operator  over  n  sub  intervals   of    length  mh.      Then 


[x-P  xjlg    £  K(m)    sup        v   i    ' 
n  (K6^h       6P 


where  K(m)  is  a  constant  depending  only  on  m. 
Proof: 

By  lemma  2.2.2 
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M,[x-P™x]    ^  m||  L  ||u)(x,h)  (3.3.1) 


1  n  m1 


In  order    to   bound   the   Holder   constant,    let    px   =  x-P  x,    and   for   two 
arbitrary   points   s..,    s„    consider    the    ratio    |  px(s1 ) -px(s_)  | /|  s1 -s_|     .      If 
h  ^    ] s 1 - s „ J     then,    using    lemma   2.2.2 

|px(s    )-px(s    )|  2|px(t)|  «J(x,h) 

—       — *  supj-— —  *  2nj|Lj|— g-  (3.3.2) 


sx-*2\  IVS2' 


If    |s..-s„|<h  and    these    two    points    lie    in    the   same    subinterval  of 
length  mh,    then 

|px(s    )-px(s2)|  R 

i 5 £    lsrs2l       n^(s1)-x(s2)l+lLmx(sl)"LmX(S2)^         (3-3.3) 

|srs2| 

m 

Now    |Lmx(Sl)"LmX(S2)l     =    I      S    x(tk)(^k(si)^k(s2))l 

k=0 


=    |      S    (x(tk)-x(t))(j£k(Sl)-jfck(s2))| 

m 
for   an   arbitrary   fixed   point    t    in   the    interval    since     £  ^(t)    =    *• 

k=0  k 

m 
£  w(x,mh)     S    U    (s    )-4    (s    )| 

k=0     k     i       k     z 


Further  j£    (t)    is   a   polynomial   of   degree   m  and    this    is   dif f erentiable .      Thus 
Uk(s1)-^k(s    )|    <■  max|4  '  (t)|  .  |  s,-s2|  .      A  bound   for    this   derivative   can  be 
crudely   estimated   from   the   expression   for  I    (t)    in    lemma   2.2.2.      It    is    the 
sum  of  m  terms   whose   numerator    is    bounded   by   ml  (l+Hyf.  .  .  I ;_ ,  '^i"1"-  •  ,+m^n 
Che  denominator    is   k ! (m-k) !h    .      Thus 
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m 


s  \z:(t) 
k=o    k 


1  s 


m. 


,m 


k=0 


Hence  from  (3.3.3)  we  conclude 


px(s1)-px(s2)|         u)(x}|  s1-s2|) 


VS2 


1-3 


srs2! 


srs2| 


o —  +  uo(x,mh).2    (2+logm)  . - 


uXx,|  S1-S2I  ) 
|srs2|P 


.2m.(2+logm)^^-  (3.3.4) 


If    |s1-s9|<h  but   the    points    lie    in  different    intervals    they  must  of 


course    lie    in  adjacent    intervals.      Let   L     be    the   operator    in  one    interval 


m 


with  basis  polynomials  i,    and  nodes  t,  ,  and  L  be  the  operator  in  the  adjace 


nt 


interval  with  basis    polynomials  I.     and   nodes    t '  .      Suppose    t       <s  <t      = 

r      J  k  k  r  m-1      1     m 

t'<s  2<t[-      Then 


m 


k=0 


m 


S    (x(tk)-x(tm))Xk(s1)-(x(ti;)-x(t;))ik(s2)| 
k=0 


m-1 


m 


£   u)(x,mh)(  S    K  (s.)|    +S    U,  (s_)|) 


k=0 


kv    1 


k=l 


kv   2' 


m-1 


=   o)(x,mh)  Z    (Uk(s1)|    +   |^+1(s2)|) 
k=0 


(3.3.5) 


NowUk(s1)|+|Jek(s2)|    = 


(s1-t0)...(s1-tk_1)(srtk+1)...(s1-tm) 


(Vt0)...(tk-tk_1)(tk-tk+1)...(Vtin) 


+ 


(S2-t')...(s2-tk)(s2-tk+2)...(s2-ti;) 


ct'     -t'^.-.ft1     -t'^rt'    -t'    ">      ft1    -t'} 

k  k+l      o  k+l     kM   k+l     k+2J         ^k+1   V 


m!    ,m-l.  (s0-t')   m!   ,m-l 

— rh        (t   -Sl)  2      oVtt  h 


m-k 


"m     1 


-  + 


k+l 


using   the   properties 


iuni 


k:(m-k):hul  (k+l)!(m-(k+l))!hul       of   s1,s2 
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<Vsi>    iv^P 


k:(m-k-l)!hL.  .  ,2   "      2  J 
v     '    (m-k)      (k+1) 


But  m-k  ^  1  and  k+1  s  1  for  the  range  of  k  in  (3.3.5)  and  thus 


m 


VS1      m-1    S2'S1 


k!(m-k-l)!     h     mCk      h 


m~  *  - 1     - 1 
Thus  since  EC,    =2    ,  and  substituting  this  result  in  (3.3.5)  we  obtain 

k=0  k 


px(s1)-px(s2)|  |x(s1)-x(s2)|  |Linx(s1)-Linx(s2) 


srs2 


10 


srs2i 


3 


srs2 


iP 


0)(x,|  s1-s2|  ) 


srs2 


U)(x,|  si_-s2l  ) 


srs2l 


+  m  a)(x,h)  .ml 


2_m-l   uo(x,h) 
+  m   2  — *—£ — <- 


-1  |srs2 


1-3 


(3.3.6) 


Thus    collecting   the    results   of    (3.3.2),    (3.3.4),    (3.3.6)    we   obtain 
the   estimate 


px(s    )-px(s») 

MJ  x-P  x;PJ     =      sup  

2  n  rv-l  I  I  I 

0<ls1-s2l  Is].'    2' 


3 


*  Kl(m)    sup     ^fl 


(3.3.7) 


0<6^h       6 


m-1. 


where   K.(m)    =  max[2m||L  ||  ,  (  l+max{  2  (2+logm)  ,m}  m2    "    )] 


Although   this   estimate   may   be    able    to   be    improved    (for   example   more 

a 

careful   analysis   of    the    case   m=l   yields    a   coefficient   of    1   for   ou(x,h)/h      in 


(3.3.4),  and  hence  K, (1)=2),  it  is  sufficient  for  our  purposes 
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Finally  by    (3.3.1)    and    (3.3.7),  and   assuming  h^l 


'x-pmx||3    =Ml[x-P™x]    +M2Lx-prax;3] 

<r    /  Ik    II   _l  «.   /  w  oo (x, 6)  _  uu(x,6) 

^    (m||L J|    +  K,(m))    sup        -s  R?    '  =  K(m)    sup        ■>  q    < 

(K6^h   6P        CK6^h   6P 


Corollary  3.3.4 

1 

If  P  x(t)  is  the  piecewise  polynomial  approximation  to  xeH(a)  then 

l|x-Pmx|L  ^  K(m).M0Cx;a]ha"P,  (Kp^c^l. 
n  p         l 

Proof; 

If  xeH(a)  then  uo(x,6)  ^  M2[x;a]6Q'.   Thus  we  have  uo(x,6)/6   ^ 
M~[x;c/]6    and  the  result  follows  by  the  theorem. 

From  the  theorem  it  is  concluded  that  the  piecewise  polynomial 
P  x(t)  converges  to  its  generating  element  on  the  subset  h((3)  c  H((3).   It 
follows  from  the  corollary  that  this  subset  includes  functions  belonging  to 
any  space  H(Q')  ,  0<3<Q'^1.   Further  it  follows  easily  from  the  theorem  that  the 
convergence  is  uniform  over  any  equi-H(P)  set.   Thus  the  convergence  is  uni- 
form over  a  bounded  set  in  H(Cf)  . 

However  the  result  is  inconclusive  with  regard  to  a  function  be- 
longing to  the  space  H($).   This  contrasts  to  lemma  2.2.2  for  the  space  C 
where  P  x  converges  to  x  for  any  xeC.   That  the  result  for  the  Holder  space 
cannot  be  improved  in  general  may  be  seen  by  considering  a  function  belonging 
to  H(3)  but  not  to  any  smaller  space.   It  is  sufficient  to  consider  the  case 
m=l  of  polygonal  approximations. 
Lemma  3.3.5 

If  xeH((3)  is  of  exact  order  (3  then  its  polygonal  approximation, 
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P  x,    does   not   converge    to   x   in   the   norm  of   the    space  H(|3),    i.e., 

lip  x-xL  -4->0   as   n-*0 

ii    n       up 

Proof: 

Let   p   x   =  x-P  x.      We   use    the   notation  of    the   definition  above    that 
rn  n 

x(t)    satisfies   a  Holder   condition  of   exact  order  3 .      Clearly  K9<K,    and   since 

K9>0,    there   exists  e>0   such    that  (K«-e)/K1<l.      {6    }     is    the    infinite    sequence 

occurring    in   the   definition   that   x   has   exact   order  3 . 

By   choosing   a   subsequence    if   necessary,   we    can   assume    that    the 

points    t    ,    t  +6      both   lie    in   the    same    interval  [  t,  ,  t,  ,  J    of    length  h,    and 
o        o     m  k     k+1 

further    that  6      satisfies    the    condition  6  <((K„-e)/K, )    //n    QN  .h  (3.3.8) 

m  m   N     2     '      1     /  (1-p) 


Thus   p   x(t  +6    )-p    x(t   )    =  A  (x(t  +6    )-x(t.  ))(t.  ,,-t  +6    )    + 
"n    v   o     nr    rn    v   o7        h    v    v   o      m  k      v   k+1      o     m 


+    (t  +6    -t.  )(x(t  +6    )-x(t.  ,.)) 
o     m     kvo     m  k+1 


(x(to)-x(tk))(tk+1-to)-(Vtk)(x(to)-x(tk+1))} 


?-{(x(t  +6    )-x(t   ))(t._Ll-t   )-(t    -t.)(x(t  +6    )-x(t   )) 
h    v    v    o     m        v    o77V   k+1      o  o      k    v    v   o     m7       v   o77 


-   6    (x(t  +6    )-x(t.))   +  6    (x(t  +6    )-x(t,       ))} 
nr    v   o      nr       N   k  mom7  k+1 


Hence    |p   x(t  +6    )-p   x(t   )  1   ^  -r   |x(t  +5    )-x(t   )|  .h-6     I  x(t.  , . )  -x(t,  )  I 
1  rn    v   o     m7    rn    v   o7  *        h    '     v   o     m7       v   o7  '  m1         k+1.r  * 

The   first   term  >K06P  ,    and   the   second   term   is  <-^.K.hP   =  K,  {—)         .6      .      The 
2   m  hi  1   h  m 

inequality    (3.3.8)    implies    that    the    second    term   is    less    than    the    first   and 
hence 


p   x(t  +6J-p   x(t   )|    >  e& 


P 


no     m        n 


(3.3.9) 


54 

; 

Since  6  -*0  as  nr100 ,  we  may  always  choose  6   sufficiently  small  to  satisfy 

mm 

(3.3.8)    no  matter  how  large   n.      Thus   there   exists  e>0  such  that  for  every  N 

there   exists  ri>N   and  6      such   that    (3.3.9)   holds.      Hence  M„[  p  x;|3]    -f»0   as   n-*00, 

m  I      n 

and  we  have  the  required  result: 


|x-P  x||   +>  0  as  n-*°  for  x(t)  of  exact  order  |3  . 


In  the  following  section  we  wish  to  apply  these  results  on  a  con- 
tour L.   From  [33]  we  assume  L  a  simple,  smooth  contour.   This  implies  that 
the  curve  has  a  continuously  dif ferentiable  parameterization.   The  direction 
corresponding  to  an  increase  of  the  parameter,  s,  is  considered  the  positive 
direction  of  the  curve.   The  smooth  contour  is  rectifiable  [33,  p.  10]  .   For 
any  two  points  t,,t  on  L  we  have  CKk  ^|  t„-t  | /|  s  -s..  |^1,  where  k   is  a 
constant  and  s.  is  the  arc  coordinate  corresponding  to  the  point  t.,  i=l,2 
[33,  p.  47]  .   For  the  (closed)  contour  we  use  t  =t   .   If  the  interpolation 
nodes  are  determined  on  the  arc  coordinate,  then  it  is  clear  from  the  above 
inequality  that  all  the  results  of  this  section  hold,  with  the  possible 
exception  of  a  constant  factor,  on  the  contour  L. 

3.4  Operator  Convergence  in  Holder  Spaces 

With  these  results  on  the  convergence  of  polynomial  approximations 
in  a  Holder  space,  the  convergence  of  approximations  to  integral  operators 
on  a  Holder  space  are  considered,  similar  to  the  approach  in  section  2.3. 

3.4.1  Operator  with  a  Holder  Continuous  Kernel 

The  integral  operator  Kx(s)  =  I  k(s,t)x(t)  dt  maps  a  bounded  function 
(and  in  particular  a  function  in  H(g))  into  the  space  H(|3)  if 
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1 
J    |k(s,t)|dt  =  g(s)eH(P).      This  will  hold   in  particular   if  k   (s)eH(3).      P 
0  t  n 

will  be   piecewise   polynomial  operator. 
Lemma  3.4.1 

Let  K  be   an    integral   operator   on   C   such   that  k   (s)eH(3).      Let   the 
approximate   operator   K     be   obtained  by   approximate   product    integration.      Then 
the   approximate   operators    converge    pointwise    in  H((3) 
Proof : 

If K.  x-Kx|L    =  II K(P  x-x)|L 

ii     n  iig  n      \    n  /up 

1  1 

=  M1[fk(s)t)(Px(t)-x(t))dt]+M[Jk(s,t)(Px(t)-x(t))dt;P] 
1   0  0 

1  ffk(s      t)-k(s      t)|dt 

<■  {maxj  |k(s,t)|dt  +      sup  -Q n }C.U)(x,h) 

S    0  0<|srs2|<l  |srs2|P 

=  C.    K Lw(x,h),    where   C    is   constant   depending  on    the   order 

P 

of      interpolation  and  ||k|L    exists    since   k  eH($).      Since   xeC    the    result   follows 

P  ^ 

It    is   clear    that   convergence    is    uniform  on   an  equicontinuous   set. 
Since   KxeH(3)    implies    that   {Kx:||x||^l}     is    an  equicontinuous    set,    and 
IK  x:||x||^lj    is   collectively   equicontinuous,    by    lemma   3.4.1  we   have 


||  (K   -K)Kx|L    £C.||k|L    ua(Kx,h)    ^C.||k|L    sup       uo(Kx,h) 

P  3  P    ||x|l^l 

Thus   ||(K   -K)K|L    -'   0   as   n-«>    and   similarly  ||  (K   "K)K  |L    -  0. 

These   results   follow  exactly   the   same   for  xeH(3),    and   thus    these 
approximations    satisfy  Anselone's   convergence    theorem    in   the    space   H((3). 


Although   the    above    approximation   scheme   demonstrates    the   desired 
convergence    properties,    it    is   not   a   very  convenient    process   because   of    the 


L 
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necessity  of  evaluating  all  the    integrals   for   the    coefficients   of   the 
approximate   operator.      A  more   convenient  scheme  would  be   the   direct  appli- 
cation  of   the   quadrature   formula   to   the  whole    integrand  as    in  section  2.3. 
Lemma  3.4.2 

Suppose    the   kernel  k(s,t)   belongs   to   the   class  h((3)  .      If   the   approxi- 
mate  operators,   K   ,    are   obtained  by  composite,    interpolatory  quadrature 
rules,      then  K  x-'Kx  uniformly   in  H(8),    for  each  xsC,    as  n-*30 
Proof; 

MJK  x-Kx]    =  ||$(P  k  x-k  x)||    £   C.||$||[||x||u>    rh)-HMi)(x,h)]       by    (2.3.1) 
In  "    v   n  s        s   /"  '    tv   '         \    »    /  j    \  / 

Similarly  M9Lk  x-Kx;3]    £     sup  { C  .||$||[  sup|  k     -k      | /|  s   -s  J     .au  (x,h) 

n  0<|s1-s2|^l  C        Sl     S2  1     L 

+  ||x||u)(k     -k      ,h)/|s   -s    |3] 

^  1  S  ry  J-  J- 

Now   |k(s.,  ,t+h)-k(s„,t+h)-k(s1  ,t)+k(s0  ,  t)  |^do(k      ,h)-Ki)(k      ,h)^2supU)(k    ,h) 

I  Z  1  Z  S-  S_  s         '    s 

But   the   expression   is   also  bounded  by  uo(k     ,  ,  |  s.. -s    |  )-Kju  (k    ,|s--s9|)   ^ 

2supOJ(k    ,  |s    -s    |).      Thus 
t  fc        L     l 

sup  <u(k     -k      ,h)/|s   -s    |P   <:  2   sup  ~ — 

h^|  S]L-s2|^l     sl      S2  l     A  s  hP 


and 


i  iB  U)(kt,6) 

SUP  U)(k     -k      ,h)/|s1-s9|p  £  2   sup  sup — 

0<|s1-s2|<h     bl     s2  z  0<|s1-s2|<h  t         §P 


Since   keh(B),    it   is   clear   that  k  eh(B)    uniformly  with  respect   to   s   and 

s 

k.eh(B)    uniformly  with   respect   to   t.      Thus 
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a                     w(k    ,h)          o)(k    ,h) 
sup  uj(k      -k      ,h)/    s    -s    |p<2max(sup 5 ,sup 5 )^(h)-0  as   H 

(Ks.-sJsl        sl      s2  X      2  '    s        hP  t        hP 


(3.4.1) 


1      21 


Then  M„[  K  x-Kx;(3]  ^c||$||  (mo[  k;(3]  u)(x,h)+||x||  .e  (h))    and 
2n  2 

Ik  x-Kx|L    =  mJk  x-Kx]+M0[K  x-Kx;3] 
11    n  "3  In  2      n 

^  C.||*||{||k|Lu)(x,h)+||x||  .e(h)}-*0  as  h  =  -  -  0 


That    is      K  x-Kx|L-*0   as   n-00    for   each   fixed  xeC    if    the    sets   fk  }  ,    {k  }    are 
11    n  !l$  s    '        t 

both  equi-H(3).       (Compare    lemma   2.3.1) 


As    in   lemma   2.3.1,    it    is    clear   that   convergence   of    the   operators 
will   be    uniform  on   each   bounded,    equicontinuous    set    of    operands.       It    is 
trivially   verified   that   {Kx:||x||^l}     is    precompact   and   { K  x:  ||  x||^l  ,n^l]     is    col- 
lectively  precompact   for   keh(0).      Thus  ||  (K   -K) kJ|r~*0   and  ||  (K   -K)K  |L-0   as   n-*° 
as   before.      Hence   although    it   can  be    shown    that   K    is    compact    in  H(3) 
(Muskhelishvili  [33,    p.    133]),    these    remarks    and    the    above    lemma    imply    that 
these    approximations    satisfy   Anselone's    theory  and   the   convergence    theorem 
1.3.2. 

It    is    clear    that   computable   error  bounds    are   not    particularly 
practical   in   this    case.      However    the   Fredholm  equation    in  H((3)    may 
also  be    regarded   as   an  equation    in    the    larger   space    C.      Further    the    integral 
operator    transforms   xeC    into   KxeH((3).      Hence   we   may  apply   the    theory   to    the 

equation   in   C   and   obtain   computable   error  bounds    for     x-x    I .      The    following 

r  "    n 

lemma  shows  that  the  sequence  {x  }  does  in  fact  converge  in  the  space  H((3)  as 
required  but  the  rate  of  convergence  in  the  H(3)  norm  is  again  difficult  to 
estimate . 
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Lemma  3.4.3 

If  the  Fredholm  equation,  x-Kx=y,  is  in  the  space  H(|3)  and  x   is 

the  solution  in  the  space  C  of  the  approximate  equation,  then  x   converges 

to  x  in  the  H(P)  norm  if  {k  },{k  }  are  both  equi-H($). 

Proof: 

We  have  x  -K  x  =y;  x-Kx=y.   Thus 
n  n  n 


x  -x  =  Kx-K  x   =  (K-K  )x  +  K  (x-x  )  (3.4.2) 

n         n  n   v   n'     nv   n'  v     ' 


By  lemma  3. 4. 2, ||  (K-K  )x|L-0  for  xeC  if  [k  },{k  }  are  equi-H(P) 

up  st 


Kn(x-xn)||3    =  M^x-x^+M^x-x^p] 

£   sup  $ P  k   (x-x  )|    +     sup  IfP   (k     -k     ) (x-x   )l/|s,-s„ 

sri      n  sv        n' '       „j  L,i      n     s.      s   '  ny  '    '    1     2 

s  (K  s,-sJ^l  1        2 


3 


<■  ||$||||l  |||k||||x-x  |H|$|[||l  ||M0[k;P]l|x-x  II,  since  [k  }  is   equi-H(3) 

ii    ii ii   mn ii    mm        n"    "    ""    m"    2                    n1"                   t  n          v   ' 

^  ||$||||l     ||k||Q||x-x  ||-»0  as   rr*00   since  ||x-x  \\-*0   in   the  space  C  by 

M     nil     mMii     iipu          nii                                              II          nii  r                    j 

theorem  1.3.2. 

Hence  taking  norms  in  (3.4.2) 


x  -x  L  ^   (K-K  )x  L  +  K  (x-x  )  L-0  as  tr*0 . 
i  n  up   ii  \   n/  up   ii  nv   n/M3 


Computable  error  bounds  for  the  second  term  may  be  obtained  using  theorem  1.3.2, 
but  computable  error  bounds  for  the  first  term  are  again  a  problem  as  in 
lemma  3.4.2. 

3.4.2  The  Hilbert  Operator 

Various  schemes  have  been  proposed  for  the  numerical  evaluation  of 
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principal   value    integrals    (e.g.,    Longman  [28]  ,   Delves   [  14]  )  .      If   the   singu- 
larity  varies   we   have    the    case   of   principal  value    integrals   with  a   Cauchy 
kernel,    in  which  we    are    interested.      The   methods   of   Bareiss   and  Neuman  [9]     and 
Stewart  [42]    both   consider   the    finite   Hilbert   transform   (over   a    finite    inter- 
val  on   the    real    line).      Appropriate   modifications   are   made    for   the    end-points. 

Since   Stewart's   method    is   essentially   an   application  of    the   method 
of   approximate    product    integration,    the    primary   question    is    that   of    the    con- 
vergence  of   a   sequence   of   Cauchy    integrals.      It    is   easily   seen   that    if   { f   } 
is   a   sequence   of    functions    in  a  bounded   set    in   a  Holder   space    then   the 
(pointwise)    limit   of   the    sequence   also  belongs    to    this    same   bounded   set. 
This,    in   conjunction  with   the   Privalov    theorem    (theorem   3.1.1),    asserts    that 
the    limit   of   a  bounded   sequence   of  Holder    functions    is    Cauchy    integrable. 
However,    it    is    not   true    in   general    that   the    limit   of    a   sequence   of   Cauchy 
integrals    is    the   Cauchy    integral   of    the    (uniform)    limit   of    the    sequence,    or 
even   that    the   uniform   limit    itself    is   Cauchy    integrable.      An  example    is    given 
by   Stewart   [42]  . 
Definition  [  42]  : 

The    sequence   {  f   ]    of    functions   on  [a,b]     is    strongly    secant    con- 
vergent  on   the    subset   [[a,b]]    =  {(s,t),    a^s,t£b,    s^t}     iff    the    corresponding 
sequence   {f   ]    of   difference   quotients   of   {f   ]     is   uniformly   convergent   on 
[[a,b]]    where 


Fn(s,t)    =    (fn(s)-fn(t))/(s-t) 


Note  that  the  limiting  case  of  the  difference  quotient  is  the 
derivative  of  the  function  and  thus  in  a  formal  sense  this  condition  requires 
the  uniform  convergence  of  the  sequence  of  derivatives  of  f  .   However  it  is 
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immediately  evident   that   strong  secant  convergence    is  merely  convergence    in 
the  H(l)   norm   (compare   section  3.3). 

Stewart's    principal   result    states    that    if   { Pn3     is   a   sequence   of 
polygonal  approximations   to   a  continuously  different table   function,    f,    then 
the   sequence   { pn}    is   strongly   secant  convergent  and  hence   the   sequence 
(Hp   (s)}    converges   uniformly   to  Hf(s)   where  H   is    the  Hilbert  operator.      In 
the   present  notation  this   may  be   restated  as    the   fact   that   if   ffc'    then   the 
polygonal  approximations,   { Pnf}  ,    converge    to   f   in   the  H(l)   norm.      This   follows 
at  once   from  section  3.3   since   a   simple   calculation  shows   that 


M 


[f-P  f;l]     =  M9[p    f  ;l]    ^   max 


dkf<fc) 


£  u)(fjh), 


and  thus 

Hf-pflL  =  Up  £!U  =mi[p  £]4m9[p  f;i3  ^  uj(f,h)-ho(f;h) 

ii    n   1   "  n   1     J-  n    &.      u 


where  u)(fjh)  is  the  modulus  of  continuity  of  the  derivative  of  f. 

Stewart  concludes  that  such  a  sequence  of  approximations  to  a 
Cauchy  integral  will  converge  uniformly  if  f  is  "sufficiently  smooth"--a 
result  effectively  obtained  in  Hardy's  second  paper  [  19]  ,  although  not 
exactly  in  this  form.   However,  the  transformed  functions  have  greater  smoott 
ness  properties  than  merely  being  continuous  and  so  uniform  convergence  does 
not  ensure  convergence  in  the  norm  of  the  range  space  of  the  operator,  as 
required  for  the  abstract  convergence  theorem.   Thus  considering  an  approach 
for  the  Hilbert  operator,  H,  in  a  Holder  space  H(&)  similar  to  that  in 
section  2.5,  we  show  first  that  the  operator  is  bounded. 

Lemma  3.4.4 

Let  H  denote    the  Hilbert  operator  on   the   space  H(0),    0<P<1-      Then, 
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||H|L  -  sup  1 1 H  1 1 R  ^  C,  where  C  is  a  constant  depending  only  on  the  contour 

MlpSi 

of    integration   and   the    index  |3 . 
Proof: 

Let    the    integration  be   over   a  closed   contour,    L,    as    in   the    Privalov 
theorem.     Write,    as   in  the   proof  of   this   theorem 


Hx(s)   =  J.  |^1  dt  =  rri  x(s)   +   r  x(t)-x(s)dt,     Then 

L     L"S  ^  Lt  t""S 


M.Chx]    =   sup|Hx(s)|    ^  tt  M,[x]    +  C1M0[x;p]    where   C,    =   supP    It-sT"    Idtl 
1  c    '  ■  1  1   l  I  JL  ' 

b  s 

Similarly   from  the    proof   of    the    Privalov   theorem  we   obtain 

M„[Hx;$]    £   CpM„[x;p]  ,   where   C?    is    a   positive   constant    independent   of  x  but 

dependent  on   the    shape   of    the   contour  L .      Thus 


!|Hx|L    =  Mj[Hx]+M2[Hx;p]     <,  C(Mj[x]+M2[x;3]  )    =  C||x|L    where   C=m3x(n  ^-K^)  . 

Theorem  3.4.5 

Let  H  be  3  Hilbert  operstor  on  H((3).   Let  the  spproximste  operstor, 
H  ,  be  obtained  by  approximate  product  integration  (polynomial  interpolation 
of  the  operand).   Then  H  -H  pointwise  in  H(|3)  for  each  xeh((3),  and  uniformly 
on  each  equi-H(3)  set  of  operands. 
Proof: 


|Hx-Hnx||e    =  ||H(x-P™x)||p    ^  ||H||p||x-P™x||p    *   K(m)||H||pU)(x  ,h)  /h3 


using   lemma   3.3.3.      Since   ||h|L     is   bounded  by    the    previous    lemma,    the    approxi- 
mate  quadrstures   converge    in   this   norm   for  esch  xeh(P).      Trivislly   the    con- 
vergence   is   uniform  over   sn  equi-H(3)    set. 
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In  particular  we  have  convergence  in  this  norm  if  xeH(a)  C  H($), 

R       rv  -R 

Ck->3  since  then  uo(x,h)/h  £  Bh    ,  B  constant.   Alternatively,  for  xeH($),  we 
may  regard  this  as  indicating  convergence  in  any  coarser  topology. 
Corollary  3.4.6 


H  -»H  pointwise  in  H(y)  for  each  xeh(3),  y<$  . 


Proof: 


||hx-H  xll      ^  M,[Hx-H  x]+M„[hx-H  x;y]^M,[hx-H  x]-M0[hx-H  x;3]=||Hx-H  x|L    for  y>$ 
11  n  "y  1  n  2  n      T        1  n  2  n    '        !l  n  MP  ' 

Stewart's   result    is   a   particular  case   of   this   corollary  where 
H(Y)    =  h(0)    =  C. 


The   reciprocal   relation  of    the  Hilbert   transform   implies  H 


-1 


-H. 


-ll 


Thus   ||H      L    =     H  L      <  C  by   lemma  3.4.4. 

P      P  • 

Theorem  3.4.7 

The  approximate  operators,  H  ,  do  not  converge  pointwise  in  the 
H((3)  norm  if  xeH((3)  has  exact  order  3. 
Proof; 

lip  x-x|L  =  ||h~  (H(P  x-x))|L  ^  ||h"  ||J|h  x-Hx|L.  From  lemma  3.4.4 

II      n  IIR  M  \      \     n  y/|IR  II  URN      n  MR 

it    is    clear    that  ||h|L    ^  0,    since   ||Hx|L    =  tt   for  x(t)    =    1.      Thus 

||h  x-Hx  L  ^  t;.||p  x-x|L   +*  0  as  n— <*> ,    for  xeH(3)    of  exact  order  3   by   lemma  3.3.5. 
i    n         up        c  "    n        'p  v  *  J 


To  conclude    this   section  we   consider   the   most   direct  method   for 
the   numerical  evaluation  of   the  Hilbert   transform.      Bareiss   and  Neumann  L9] 
propose   the    "method  of   symmetric   pairing"   in  which   the   principal  value    integral 
is    reduced   to   an  unconditionally   convergent    integral   by    "pairing"   the    function 
values   of    the    integrand   for   the    arguments    that   are    symmetric  with  respect   to 
the    singularity.      Thus 
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*<.,.  tf  |^t  -  »»{(;*  :+l)^ 


b    , .  v  ,  .      rs  s-e 

a  s+e 


=    li°{r      X<S+U>-X^S'u)du  +  af  X<S^U)du, 
e-»0   lJ  u  -J         u  J 

e  m 


,+1    if  s-a<b-s 

where   m=min(s-a,b-s)  ,M=max(s-a,b-s)  ,  O  —  \         ... 

v         '        ''  v  L-l    if   s-a>b-s 

The    trapezoidal   and  mid-point   rules    are    applied   to   the    resulting 
unconditionally   convergent    integrals    (effectively   assuming   the    integrand 
constant    in   the    first    interval   for    the    improper    integral).      Thus   expressions 
are   obtained   for   Hx(a+kh)    and  Hx(a+(k+~-)h)    respectively,  where   h   =    (b-a)/n. 
Thus    this    is   essentially   an   application   of   the    process   of    ignoring   the    singu- 
larity  considered    in   section   2.5.      That    is,    for    the    trapezoidal   rule,    the 
function  Hx(s)    is   evaluated   at    the    rational   points   a+kh  by   a   compound   trape- 
zoidal  rule  with    the    singularity   at    zero    ignored.       (For    the    purposes   of    this 
formal   discussion  we   may   assume    (as    is    done    in  L  9]  )    that   x   vanishes   at   the 
end-points    in   order   to   avoid    the   end-point    problem) . 

However    in  [9]     the    authors    do   not   consider  whether    these    approxi- 
mations   converge    to    the    integral    for   each  k,    nor  whether    the    interpolation 
of    these    points   converges    to  Hx(s)    as   n-°° .      By  applying  Miller's    results   [32] 
we   show  that    the    convergence    is    uniform   in   s. 
Lemma   3.4.8 

Let   the   approximate   operator,   H    ,   be   defined   from  the    interpolation 

n 

of    the    points   obtained  by   applying   the   method  of   symmetrical   pairing  with   a 
trapezoidal   rule    to   the   Hilbert   operator  H.      Then  ||h  x-Hxlj-0   as   n-*00   for 
xeH(3)  ,    O<0£1. 
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Proof: 

Let  H     denote    the   operator  obtained  by   applying  the  method  of 
symmetric   pairing  at   the   rational   points   t,  ,    k=0,l,...,n.      By  Miller's 
result  [32,    theorem  2]  ,    the  compound   interpolatory  quadrature   rule  converges 
to   the   singular    integral   at    the   rational   points,    t    ,    if   the    integrand  can  be 
dominated  by  a  monotone    integrable   function   in   the    (one-sided)    neighborhood 
of   the   singularity.      Since   the    integrand,    (x(s+u)-x(s-u))/u,    is  bounded  by 
M~[x;|3]  .u  for  xeH((3),    and   this   bound   is    integrable,    it   follows   that  this 

quadrature   converges   for  each  fixed   t,  .      That   is  H  x(t,  )-»Hx(t,  )for  each  k, 
as   n-*00,    since   the   convergence  of   the   remaining   integral   follows   as    in 
section   2.2. 

We  need  to   show  that   the   convergence   of   the    improper    integral 
portion   is   uniform  with  respect   to  k.      Using  the   notation  of  Miller's   result 
for   the   error  [32,    theorem  3]  ,    given  e>0,   we   can  find  H  such  that 

H  H 

|J*   Pnf (t)-f (t)dt|    ^    (1+K)J  F(t)dt  <  e 
0 

where  K  is  a  positive  constant  depending  on  the  integration  rule  and  F(t)  is 

the  majorizing  monotone  integrable  function  for  the  integrand  f.  Further, 

we  can  then  choose  h  sufficiently  small  so  that  for  an  arbitrary  point 
t,  =  kh,  using  theorem  2.2.4, 

kh 
|f   Pnf(t)-f(t)dt|  *  OJ(f,-)(kh-H)  £  U)(f,-)((b-a)/2-H)  <  S 

where  the  modulus  of  continuity  of  f  is  calculated  over  the  interval 

[H,(b-a)/2]  ,  where  the  integrand  is  continuous.   That  is,  for  n  sufficiently 

kh 
large,  |J   p  f(t)-f(t)dt|  <  2e ,  independent  of  k. 
0   n 
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For   the    integral  above, F(t)    =  M„[x;(3]  t  anduu(f,-)    = 

uu((x(s+u)-x(s-u))/u,l/n)    ^    (M2[x;0]tT    +  Mj[x]  .h/H)  .2/H.      Thus,    in   this    case, 

maxjii  x(t   )-Hx(t   )|    <   2e    for   n   sufficiently    large, 
k^n 

The    corresponding   approximate    functions    are   obtained  by    interpolat- 
ing  the    points  H  x(t,),    k=0,l,...,n.      This    gives    the    piecewise    linear   function 

cp        H  x   =  H  x(s)  ,   where  cp         is    the    linear  operator   from   the   space,    m    ,    of 
onnv/  To  r  r        '      n 

bounded   sequences   of    length  n,    to    the    space,    C,    of   piecewise    linear   functions. 
Hence    in   the   usual   sup  norm  we   obtain 


|h  x-Hx||    ^  llcp"   H  x-P  Hxll    +  jlp  Hx-Hx||  (3.4.3) 

I      n  II  ,IT0         n  n         ,1  i|      n  II  V  / 

<.   max|H  x(t  )-Hx(t  )|  +   (Hx,-),  since  HxeC. 
k^n 


Hence  ||H  x-Hx||—  0  as  vr^    and  we  have  pointwise  convergence  of  these  operators 
in  the  sup  norm  for  each  xeH($). 


However  uniform  convergence  of  approximations  to  the  transformed 
function  was  obtained  previously  (corollary  3.4.6),  and  this  does  not  imply 
convergence  in  the  H(3)  norm.   It  was  necessary  to  assume  xeH($)  to  obtain 
this  result  and  the  error  bounds  above  involve  the  H(|3)  norm  of  x.   Compare 
theorem  3.4.7. 
Lemma  3.4.9 

The  approximations  obtained  by  the  method  of  symmetric  pairing  do 
not  converge  pointwise  in  the  H($)  norm  if  xeH(|3)  has  exact  order  3<1. 
Proof; 

Considering  the  inequality  (3.4.3)  in  the  H((3)  norm,  we  have  in 
an  alternate  form 


i 
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P  Hx-Hx  L  £   max  H  x(t.)-Hx(t.)   +  H  x-HxL 
n     "B   n  .  '  n  v  k    v  W  '    "  n    Mp 
k^n 


By  lemma  3.2.1,  if  x  has  exact  order  B  then  Hx  is  also  of  exact  order  B.   Thus 

■ 
by  lemma  3.3.5  ||p  Hx-Hx|L  +♦  0  as  rr*° .   Hence  the  above  inequality  implies 

||h  x-hx|L  +>  o. 

11  n    "B 

That  is  a  direct  approximation  method  leads  to  the  same  type  of 

difficulties  as  found  for  approximate  product  integration  in  theorem  3.4.7 

and,  in  general,  I  H  x-Hx|L  +»  0  as  rr^°   for  xeH(B). 
o  '  "  n    "B 

3.4.3  The  General  Operator  with  Cauchy  Kernel 

Consider  the  case  of  the  general  singular  operator,  K,  with  Cauchy 
kernel,  k(s , t)/ (t-s) .   We  consider  two  approximation  schemes.   Firstly  where 
the  kernel  is  integrated  exactly  and  approximate  product  integration  is 
applied  in  a  form  similar  to  lemma  3.4.1.   Secondly  where  only  the  "singular 
part"  of  the  kernel  is  integrated  exactly  and  the  continuous  part  of  the 
integrand  is  interpolated  for  the  approximate  product  integration  scheme. 

The  kernel  is  not  absolutely  integrable  and  it  is  first  necessary 
to  establish  that  K  is  a  bounded  operator  on  the  Holder  space.   This  depends 
on  the  generalized  Privalov  theorem  where  the  density  function  depends  on  a 
parameter  [33,  p.  49]  ,  [38,  p.  443]  .   Pogorzelski's  result  shows  that  if 
cp(u,  t)eH(a  ,B) ,  (Kc^l,  (KB<1,  then,  on  the  smooth  contour  L,  the  singular 
generalized  Cauchy  integral  $(u,s)  =  [   ^  I  »,  ,/dt  satisfies  $  (u,s)eH(a  ]B)  where 

Jj 

a'  is  an  arbitrary  positive  constant  less  than  ex.      Thus  in  the  case  that  u=s 
we  have  $  (s)  =  $(s,s)eH(a')  for  o^B  and  $(s)eH(B)  for  a>$  .   Hence  in  order 
that  the  operator  K  maps  the  space  H(B)  into  the  space  H(B)  it  is  necessary 
to  require  that  k(s  ,t)eH(a  ,y)  ,  where  a£>B  and  y^P  •   Note  that  this  condition 
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implies    that   k  sH(y)    uniformly  with   respect   to   s,    and  k  sH(Qf)    uniformly  with 
respect    to   t. 
Lemma  3 .4. 10 

Under   the    assumed   conditions,    the    operator  K   is   a  bounded,    linear 
operator   on  H(P).      That    is   ||k]|q    =      sup     ||kx|L    ^   C   for   some   constant   C. 
Proof: 

Following   the    proof  of    lemma   3.4.4,    but   using   the   generalized 
Privalov  theorem  in  this   case 


M,[Kx]    £  tt  Mj[kx]    +  Cl  M2C  (kx)    ;P]  ,    and 


M^KxjP]    *  hm2[  (kx)t;P]    +  C^L  (kx)g  ;P]    +  C3(m2[  (kx)g  ;p]    +  M2[  (kx)  fc ;a]  } 


where   C-,    is   also   a   constant    independent   of   k,x   but   depending  on   the   contour 
L.      Thus  with    the   bounds 


MjTkx]    =S  MjEk]  .Mj[x]  ; 

M2[(kx)t;p]    S  C4  M2[(kx)t;a]    *  C4  M2[kt^]  .M^x]    since  a>p  ; 

M2[(kx)s;p]    S  M^k]  ,M2[x;3]    +  M^xJ  .M2[k    ;0]  ,    we   obtain 

Kx||      =  Mj[Kx]    +  M2[l<x;p] 

S  {n  ML[k]    +   (C3-mc4)M2[kt;cv]    +   (Cj+C^+C^M^  kg  ;P]  Jm^x] 

+    (C1+C2+C3)M1[k]  .M2[x;p] 

£   C,||k||      R   M.[x]    +  cJ|k||M9[x;P] 
5"      a,p      1  6"    "    2 
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where  )|k||    denotes   the  upper  bound  of   the   function  k(s,t), 


k|L  0    denotes  ||k||-H40[k    :a]  +M „[k   ;3]  ,    since   k  eH(v)    implies   k  eH(B) 

dp  M    "      2      t  2      s  s  s 


Y^3, 

C,.   =  max(n,C3-HTC,  ,C1+C2+C3)    and  C&   =  €,+0+02. 

Thus      ||kx|L   ^  cllxIL    where   C  =  max(cj|k||      D,    cjlkll). 

II  Jig  II       llg  \      5II       IIq,   fi  gl|       11/ 

Lemma  3.4,11 

Let  K  be   a  bounded   integral  operator   on  H($)   with  a  Cauchy  kernel, 
k(s,t)/(t-s)   where   keH(a,$),  a>|3  .      Let   the   approximate  operator,   Kn,    be 
obtained  by  approximate   product   integration  where    the  operand   is   approximated. 
Then  the   sequence   of   approximate   operators   converges   pointwise  on  H(|3)    for 
each  xeh(3)  c  H($),    and  uniformly  on  each  equi-H(3)    set  of  operands. 
Proof: 


||Knx-Kx||p    =  ||K(Pnx-x)||      £  Nlg||PnX-x||      =  |t^|      KW    sup2i^  by 

lemma  3.3.3.      Since  ||K|L    is   bounded  by  the   previous    lemma,    it   is   clear  that 
the   approximate  operators   converge   pointwise    in  H(3)    for  each  xeh(j3),    and  in 
fact  converge   uniformly  on  each  equi-H(|3)    set  of  operands. 

From  this    lemma   it  cannot  be   asserted   that   the   approximate   operator! 
converge   pointwise   for  xeH((3).      But  we   can  make   the   statement: 
Lemma  3 .4. 12 

Suppose    that  k(s,s)    does   not  vanish  anywhere  on  the   contour  L.     The: 
the   approximate  operators,   K   ,    of   the   previous    lemma  do  not   converge   point- 
wise    in   the  H(3)    norm   if  xeH((3)   has   exact  order  |3  . 
Proof: 

Write 
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P  x(t)-x(t)      ,  k(s,t)-k(s,s) 

k(s,s)Jn  t_s dt  =  J-i^(Pnx(t)-x(t))dt-J — (Pnx(t)-x(t))dt 

Li  ij  L 


Thus  writing  b(s)  =  k(s,s)  we  have 

||b.H(Pnx-x)||p   ^  ]|Knx-Kx|lp   +  ||L(x-Pnx)||3  (3.4.4) 

where  H  is  the  Hilbert  operator  and  L  is  the  integral  operator  with  the 
weakly  singular  kernel  i,(s,t)  =  (k(s  ,  t) -k(s  ,s))  /  (t-s)  . 

The  term  with  the  operator  L  converges  by  the  argument  of  lemma 
3.4.1.   Although  the  kernel  does  not  belong  to  H(B),  the  same  argument  gives 

\k  (s      t)  -SL  (s      t)  I 

||L(x-P  x)||      S  {tnaxf|i(s,t)||dt|    +     sup  J ■  *'    q    ir    2' |  dt|  }u>(x,h) . 

P  s    L  0<|s    -s    |L        |Srs2'P 

Since  i,(s,t)    is   weakly   singular   the    first    term   in    the   braces    is   bounded.      By 

Pogorzelski's    lemma  [38,    lemma   1,    p.    492]    the   kernel,  4(s,t)   may  be  expressed 

*  i         i  1-e 

in   the    form    (k(s , t) -k(s ,s) ) / ( t-s)    =  k    (s,t)/|t-s|         ,    where   e    is    an   arbitrary 

positive   number    less    than  B    and   k    (s ,  t)eH(o/-e  ,B  -e )  .      Thus    the    second   term   in 

the   braces    is    bounded  by 


SVP      i  Is  -s  IbU^V^I 

<   s,-sj     |S1   S2,P    L 


I  1-e    i  i 1-e 

t-s.,  t-s, 


dt     + 


I  ~2'  i      °1'  '2 


r»i    *  *  I         dt 

+  r  k  (s.,t)-k  (s,,t)' 


1'    '         v-2'^i  ,  ,1-G 

L  t-sj 


i 


1-B-e   „    r  dtl 


s     sup  <    s,-s„  .D.       k ' ; 1 

'    l"s2'  L     I'"8!'         It_s2l 


I              \a-&  -B        r ,  *          ir          dt 
+    |srs2|  P.M2[kt;a-e]J— -1 L_ 


r 

_  ej 

L  | t-s2| 
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where  D  is  the  absolute  bound  for  k  .   These  integrals  are  clearly  bounded 
(compare  Muskhelishvili  [33,  p.  136]  for  a  discussion  of  the  first  integral), 
and  thus  by  choosing  e  ^  min(l-|3  ,a-3)  ,  the  entire  term  is  bounded.   It  is 
thus  concluded  that  L(x-P  x)L  ^  E  (D(x,h)  for  a  constant  E,  and  hence  this 

ll   v     n  /lip  \     / 

term  converges  to  zero  as  n-*33  for  xeH(|3)  . 

On  the  other  hand,  by  theorem  3.4.7,  ||h  x-Hx  L  ■**  0  for  x  of  exact 

J  "    n  "  p 

order  |3  .      Since  we   know    (corollary   3.4.6)    that  ||H  x-Hx||   -*  0,    this    implies 

M?[h  x-Hx;P]    I*  0  as   n-03 .      Thus   there   exists  e>0  such   that   for   all  N,    there 

exists   n>N  and   two   points   s..,s~    such   that     a    (s.)-a    (s„)/s,-s_|      >  e,   where 

r       1'  2  '  nv  V      n  2/ '  '  1   21 

o    (s)  =  H  x(s)-Hx(s).   Using  this  result  we  obtain 

lb(Sl)qn(Sl)-b(s2)qn(s2)|  ^C^^CBz)    ,    b(s1)-b(s2) 

j jp =     b<8l>— j |P ~I iT"'CTn<82> 

I sl"s2'  lsl-s2l  |S1"S2' 

K,      JVS1>~VS2>1         lb(sl)-b(s2>l,       ,     ft 

b<Sl>l j [1 j [3 lCTn<S2> 

|srs2|  |srs2| 

The   right   side    is   greater   than    (|b(s,)|e    -  M«[k;"3]T|)   where    \o    (s_)|<T],    an 
arbitrarily   small  constant,    for  n   sufficiently   large   since  \\o     -*0  as   ir*30. 

J  J  O  II      nii 

Since   b(s1)    f  0,   we   deduce   that,    for   sufficiently   large  N,  M?[b(H  x-Hx)  ;(3]    ^ 
|b(s1)CTn(s1)-b(s2)o-n(s2)|/|s]L-s2|P>€/2,    and  hence  ||b.  (Hnx-Hx)||p   +*  0  for 
xeH(p)    of  exact  order  |3  . 

In   view  of   these   two   results   and    inequality    (3.4.4)    it    is   concluded 
that  ||k  x-Kx|L  *   0  for  xeH(P)    of  exact  order  3. 

The   procedure    in   lemma  3.4.11   involves   evaluations   of   the    integral 

k(s,t)/(t-s)dt    in  order   to   obtain   the   coefficients    in   the   approximate  form 
L 
A   simpler   approximate    form  could  be   obtained  by   interpolating   the  whole  nume 

ator  of   the   general   singular  operator   and   integrating  only   the    "singular  par  • 
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Lemma  3 . 4. 13 

Let   K  be   a   singular    integral   operator   on  H(|3)    with   a  Cauchy  kernel 
k(s,t)/(t-s)  .      Suppose   k(s  ,t)eH(o',Qf)    =  H(cv)  ,  ar>(3  ,    and    let    the    approximate 
operators,    K    ,    be   obtained   by   approximate    product    integration  where    the 
singular   part,    l/(t-s),    is    integrated   exactly.      Then   the    sequence    of   approxi- 
mate  operators   converges    pointwise    in   the   H((3)    norm  for   each  xeh((3),    and 
uniformly   on  each  bounded   equi-H(P)    set   of   operands. 
Proof: 

From   the    remarks    at    the   beginning  of    this    section    it   follows    that 
K   is   a  bounded,    linear  operator   from  H(|3)    to   H($).      Using   the   generalized 
Privalov    theorem   as    in    the    proof    of    lemma   3.4.10, 

M,[K  x-Kx]     ^  rr  m,[  P  k  x-k  x]    +  C,    M„[  (P  k  x-k  x)    ;p]  ,    and 
In  Ins        s  1      2    v   n   s        s      s 


M 


2L  (Knx-Kx)  ;0]    *    (TTC4+C3)M2[(PnksX-ksx)t;a']    +    (C2+C3)M2[  (P^x-^x)  g  ;0] 


where  we  may  choose  an  arbitrary  a'    such  that  (3<cf  '<a .   We  have  the  following 
bounds . 

M.[P  k  x-k  x]    <■  c{  ||x||  supua(k    ,h)+||k||uj(x,h)}    by    (2.3.1) 
l     n   s        s  s  s 

uj(k  x,6) 

M?[(P  k  x-k  x)    :0]    s  K1(m)sup     5 by   (3.3.7) 

S        S      S  L     (K6^h  6P 

uj(k    ,6)  uu(x,6) 

^  K1(xn){||x||sup     g +  l|k||sup     = — }      by    (2.3.1) 
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M0[  (P  k  x-k  x)    -cv']    £  C   sup 
2ns        s't  r^.\ 


00  ((k     -k     )x,h) 
Sl     S2 


ii  \<y' 

(Klsrs2'     l  si_s2 ' 


oj(k     -k      ,h)  |k     -k     | 


S ..         S_  S-Srt 

^  C   sup  {||x|| — ■ — r-  +  sup— ~,oj(x,h)} 

0<|s   -s2|  |s   -s J  t      |s   -sj 

by   (2.3.1) 


^  C'{2||x||  .M9[k;a]ha_a     +  M„[k;a]u)(x,h)}  by    (3.4.1) 


Thus,    collecting  all  expressions 


K  x-Kx|L    ^  KlK  x-Kx]    +  M0[K  x-Kx;3] 
n  "p  In  2     n 


^cJ|x||M2Ck;c]h^'  +C6||k||a   sup  -ifeil 


where  C,-,Cr  are  constants  depending  only  on  the  order,  m,  of  the   i 
_>   b 

interpolating  polynomial  and  the  contour  L. 
Hence  the  sequence  of  approximate  operators  converges  pointwise  in  H((3)  for 
each  xeh(3),  and  uniformly  over  each  bounded,  equi-H(P)  set  of  operands. 

Thus  again  we  cannot  assert  pointwise  convergence  for  xeH(P).  How- 
ever, as  in  corollary  3.4.6,  we  trivially  have  convergence  in  any  coarser 
topology.   In  particular 
Corollary  3.4.14 

The  sequence  {k  x(s)-Kx(s)}  converges  in  the  sup  norm  of  space  C, 
uniformly  over  a  bounded  set  in  H(|3). 
Proof: 

Considering  the  bound  for  the  M,  term  in  the  above  lemma  we  are  at 
liberty  to  choose  any  index  ^  (3  in  the  second  term. 


73 


M,[k  x-Kx]    £  n  M.[P  k  x-k  x]    +  C'    M0[  (P  k  x-k  x)     ;p/2] 
In  lnss  1      2    N   n   s        s    's 

£   C2.K(m).M  [k  x;p]tT/2   by   corollary   3.3.4 

<;   C2   K(m){||x||M2[ks;3]    +  ||  k||M2[x;p]  }  h3  /2 

n    ii      ii    n      3/2  ii    ii 

£   C0  .    k      .    x  Lb.  -   0   as   n-°° ,    uniformly   for     x  L    ^   B 

3  ii    na  n      (3  J  "    "p 


Although  singular  integral  equations  are  usually  considered  in 
the  space  H(3),  the  above  remarks  concerning  pointwise  convergence  imply 
that  requiring  the  functions  to  belong  to  h(p)  might  be  a  more  desirable 
hypothesis.  However,  it  is  not  clear  that  the  Privalov  theorem  could  be 
extended  to  the  necessary  result  that  the  Hilbert  transform  of  a  function  in 
h((3)  also  belongs  to  h(3)  (Compare  the  remark  in  section  3.2  that  Hx  does 
not  necessarily  belong   to   C    for   xe C   =  h(0)). 

Following   the    analysis    in   Chapter   2,    the    remaining   condition   to   be 
considered    is    the    compactness   of    the    operator   K.      By   lemma   3.3.1   a   subset 
of  h((3)    is   compact    iff    it    is   bounded   and   equi-H($).      Further   a   compact   set 
in  h(p)    is    clearly   compact    in   the    larger   space   H(|3).      It   follows    from   the 
above    lemmas    that  ||  (K   -K)k|L— 0    if   K    is    a   compact   operator    (in  H(3]).      Since 
the   operator   K    is    bounded    in  H(3),    it    is    specifically   required   that 
lKx:||x|L^l}     is    an  equi-H(P)    set.      Since    the   general  operator   K  may  be   ex- 
pressed  as    the    sum  of   a  Hilbert   operator   and   a  weakly   singular  Fredholm 
operator,    it   suffices    to   consider  whether    the   Hilbert   operator   satisfies 
this   condition.      But   by    lemma   3.2.1,    if   x  has   exact   order  3    then  Hx  has 
exact   order   (3.      Thus,    from   the    definition    in   section   2.3,  ou(Hx,h)/h     *•   0 

as  h-0.      Hence    the    set   {Hx:||x||    ^l]     (and   by    inference   { Kx:||x|L^l}  )    is   not 

P  p 

equi-II(p)  . 
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Hence  we  cannot  assert  either  K  -*K  in  H(3),  or  ||  (K  -K)KJjg"*0  and  * 
thus  the  hypotheses  for  the  abstract  convergence  theorem  1.3.2  are  not  satis- 
fied.  On  the  other  hand  we  do  have  K  -+K  in  the  space  h(|3)  (by  lemmas  3.4.11 
and  3.4.14),  but  it  is  not  clear  that  we  can  assert  Hxeh(f3)  even  if  xeh(|3). 

The  situation  is  similar  to  that  of  the  Fredholm  equation  in  the 
space  R  (where  CC  RC  m,  the  space  of  bounded  functions)  discussed  in 
section  2.4.   The  concept  of  a  regular  set  was  introduced  in  order  to  extend 
the  region  of  uniform  convergence.   But  this  concept  depended  upon  the  fact 
that  the  integral  operator  is  a  positive  linear  functional  and  since  the 
Hilbert  operator  is  not  a  positive  operator,  it  does  not  seem  that  this  idea 
would  provide  a  fruitful  extension.   This  difference  in  behaviour  from  the 
ordinary  integral  operator  is  apparently  a  consequence  of  the  property  of 
the  Cauchy  principal  value  integral  being  only  conditionally  and  not 
absolutely  convergent. 

3.5  Remarks  on  Investigations  of  a  Direct  Method 

Gabdulhaev  [  17]  presents  some  investigations  of  the  conditions  for 
solvability  and  convergence  of  the  method  of  constructing  solutions  to 
singular  integral  equations  by  means  of  polygonal  approximation.   Although 
this  justification  of  the  method  is  generally  based  on  the  general  approxi- 
mation theory  of  Kantovovich  [25]  ,  the  form  of  the  approximation  is  much  the 
same  as  those  considered  in  the  preceding  sections.   The  convergence  theorem 
will  be  reviewed  in  the  light  of  the  difficulties  that  were  revealed  in  the 
previous  section. 

Corresponding  to  Kantovovich1 s  notation  the  polygon  P  x(t)  is 


denoted  by  x(t)  and  the  space  H(3)  is  the  set  of  polygons  corresponding  to 
xeH(P). 
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Gabdulhaev   considers    the    singular    integral   equation  with   a  Hilbert 
kernel  of    the   normal  L  38]     type: 

1      2TT 
Kx(s)    =   a(s)x(s)    +  Ux(s)    =    a(s)x(s)    +  —  J   b(s  , t)cot (-^)x(t)dt    =   y(s) 

0  (3.5.1) 

where  a(s),  b(s,t)  y(s)sH(a),  (Kcv^l.   Note  that  the  singular  part  of  the 
kernel  in  the  above  operator  has  the  form  cot ( (t-s) /2)  ,  which  is  called  the 
Hilbert  kernel.   However  it  is  easily  shown  (compare  Pogorzelski  [38]  , 
Muskhelishvili  [33])  that  tbis  operator  is  equivalent  to  an  operator  with  a 
Cauchy  kernel,  l/(t-s),  over  a  closed  contour  L,  which  we  have  been  investi- 
gating. 

An  approximate  solution  to  equation  (3.5.1)  is  sought  in  the  form 

n-1 

of  a  polygon  of  the  type  x(s)  =  T,    c.cp,  (s),  where  cp  (s)  =  (s-s,  .)/(s  -s   ..), 

k=0 

s,  ,^s^s.  ;  (s,  ,,-s)/(s,  ,,-s.),  s.  ^s^s.  ,,;  0  elsewhere.   The  coefficients 
k-1    k    k+1      k+1   k    k    k+1 

are  determined  from  the  system  of  linear  algebraic  equations 


Kx(sL)  =  y(s.)       (1=0,1,..., n-1) 


The    latter   system  can  be    rewritten    in   the    form 


Vi  +  Xbikck-yi      •   i=0-1 "-1  <3-5-2) 

k=0 


2n                      (t-s.) 
where   a.    =   a(S;L)  ,    y.    =  y(S;L)  ,    b ik   =   2tT  i^    Ms^^cot—^ cpR(t)dt. 

We  note  that  this  approximation  scheme  is  effectively  just  an 

application  of  the  approximate  product  integration  method  considered  in 

lemma  3.4.11. 
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1  s  1   Solvability  and  Convergence  of  the  Method 

Gabdulhaev  proposes  that,  since  K  is  a  bounded,  linear  operator 
fromH(0)  into  H(p),  (0<P«»  (compare  also  lemma  3.4.10),  the  equation  (3.5.1) 
can  be  considered  as  a  linear  functional  equation  of  the  for*  Kx  -  ax+Ux  -  y 
in  the  space  H<0).   It  is  supposed  that  under  the  given  conditions,  equation 
(3.5.1)  has  a  unique  solution  x*(s)  from  some  Holder  class.   Then  it  is     II 
stated  that  it  is  known  from  [33]  that  such  a  solution  satisfies  the  conditio, 
x*(s)eH(5)  «X6<fc).   The  numbers  or  ,0  and  6  are  fi^ed  so  that  (*$<6«> .   The 
system  (3.5.2)  is  rewritten  in  the  space  of  polygons  ~H(p)  in  the  form  of  an 
equivalent  linear  functional  equation 

KX  =  P„(ax)   +  Pn(*)    »  y  (3-5.3) 

n  n 

where 

K:H(B)-H<P)    and  ||k||      =     sup  II  P^ll^ll  PjIpNIp 

x  L^l,xeH(p) 
"P 


Note    that  ||pJL*3 


Convergence   of   the   method   is   established  by  the   theorem 

Theorem  3.5.1  L 17 ,   Theorem  1,    p.    219] 

Suppose   that  an  equation   (3.5.1)    is  uniquely   solvable   for   an 
arbitrary   right   side   and   that   the    linear  operator  F  V*,)    exists.      Then 
the   approximate   solutions,   x(s)  ,    converge   to   the   exact   solution  x»    in 
the    sense   of    the   metric  H(p).      In   this   connection 


' 


IX   -x 


VM\    s\™  (J.^J*  (3.5-4) 


*||     <    (l+3||Kl|p||K"l|p).3M2(x;6).h 


To  begin  the   analysis   of   this   section  we   note   that   the   stated 
hypotheses   of   the   equation   (3.5.1)    and   the   assumed  solution   space   do  not 
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seem  to  be   consistent.      If    it    is    supposed   that   the    solution   satisfies    the 
condition  x    (t)eH(5),    0<6<b/   with   exact   order   6,    then  x"(t)    does   not   belong 
to   any   space    smaller   than  H(6).      Thus   b(s,t)x    (t)eH(o,,6)    and  by    the    generali- 
zation of   Privalov's    theorem  [38,    p.    443]     (compare   also   section   3.4.3) 
Ux    (s)sH(6)    and   does   not  belong   to   any   smaller   space  H^)  .      Further   ax     be- 

-k 

longs  to  H(6)  exactly  and  thus  Kx  eH(6)  and  not  to  the  smaller  space  H^)  . 
On  the  other  hand  the  hypotheses  give  the  right  side,  yeH(a)  C  H(6),  which 
gives    a  contradiction. 

However    if    the   hypotheses   on   the    given   functions    are   modified   to 

a(s),y(s)   e   H(a)       ;    b(s,t)   e    U(y  ,a)    ,    0<2y<Y^l  (3.5.5) 

then  Ux  e    H(cv)    for   xeH(O')    and   the   equation   Kx=y   may  be    regarded   as   an 
equation    in   the    space   11(0"). 

Further,    under    these    conditions,    the    operator   K   satisfies 
K:H((3)-H(|3)  ,    CKP^o/    (compare    section   3.4.3),    and  hence    for   yeH(3)    the   equation 
Kx=y    is    an  equation    in   the    space   H(3).      However    the    condition  yeH("y)  c:   H(f3), 
0<$<y  ,    leads    to   a   contradiction   as    in    the    argument   above,  if    it    is    assumed 
that   a   solution  exists    in   H(j3)    and    in   no   smaller   space. 

Thus    in   order    to   determine    the    appropriate    solution   space    of    the 
equation    it    is    necessary    to   determine    the    smallest   Holder    index   satisfied   by 
the    given   functions.      Hence    if   yeH(6)    then  we   may   seek    the    solution  x  £H(6) 
in  Gabdulhaev's   equation   and    the    above    convergence    theorem  holds.      However 
we   obtain   convergence   only    in    the    coarser  H (P )    topology   and    it    is    clear    that 
theorem  3.5.1   does   not   assert   convergence    in   the    solution   space   H(6). 

The    convergence    in   theorem  3.5.1   also   depends   upon   the   existence, 
and  uniform  boundednes:;    for    sufficiently    large   n,    of   the    inverse   operator 
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K   .   Kantorovich  [25,  p.  546]  notes  that  this  condition  will  certainly  hold 
if  the  operator  is  completely  continuous  (or  compact),  but  that  the  analysis 
will  still  hold  if  this  result  can  be  established  directly.   It  has  been 
observed  previously  that  we  were  unable  to  assert  that  the  singular  operator 
is  compact  in  the  Holder  spaces.   In  succeeding  theorems  Gabdulhaev  obtains 
this  uniform  boundedness  under  quite  stringent  conditions.   However  the  con- 
vergence theorem  is  still  unsatisfactory  under  these  conditions  since  the 
convergence  is  not  in  the  metric  of  the  solution  space,  H^)  . 
Lemma  3.5.2 

The  Gabdulhaev  approximations  do  not  converge  in  the  solution 
space  H(cf)  . 
Proof; 

Since  the  functional  equation  (3.5.3)  is  equivalent  to  the  linear 
system  (3.5.2)  we  have  x  =  cp  x  where  x   is  the  solution  of  (3.5.2)  and 

cp    is  the  linear  operator  from  the  space  m  to  H(a).   By  theorem  3.5.1  we 

M  ~*  *n 
have  || x  -x  ||q~*0.   This  implies  convergence  in  any  coarser  topology,  and  in 

particular  uniform  convergence  in  the  space  C,  ||x  -x  ||~*0.   By  further  special- 

izing  this  result  we  obtain  max|x  (s.)-x  (s.)|-*0,  which  may  be  equivalently 

I-*      *       iSn     ~*    -1-, 
written  max|x  (s.)-x  (s.)|-»0,  since  x  =  cp  x*  is  the  interpolation  of  the 

i^n 

points  x*(s.),    1=0,1, ...  ,n-l.      Thus  we  may  write    (where  cp.(s)    is   defined   in 

the   previous   section) , 


*      ,  n-1     *  * 

P  x  -x*||      =1       2    (x''(s.)    -  JT(s.))cp.(s)|| 
x=0 

n' l 
£  max|x'  (s\)-x   (s^]  E  9i(s)||a 

i^n  i=0 

.    *  -*  .  n_1 

=  max|x    (s  .)-x    (s  .)  |  ,    since     Scp.(s)    =  P   (1)    =   1 

i^n  ii  i=Q   i  n 
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The  right  side  converges  to  zero  by  the  result  derived  immediately  above. 
Thus  we  may  write 

P  x    -x  II    £    P  X    -X        +Ux    -x        %x  x    (s.)-x    (s.)    +11  x    -x  I 

ii     n  na    ii    n  i\(y  na  i  -,_/  v     L/  i     i.  ,lQ, 

i>n 

Since    the    first    term  on    the    right    converges    to    zero   as    rr00    and    an    application 

of    lemma   3.3.5    implies    that   |lp  x    -x  |       t»   0   for   xeHto)    exactly,   we    conclude 

1  n     la 

that  |!x"-x"||   t*  0. 

The  argument  for  non-convergence  depends,  as  usual,  on  the  fact 
that  the  polygonal  approximations  converge  to  the  generating  function  only 
on  a  subset  of  the  Holder  space.   In  terms  of  the  Anselone  theory  it  has 
already  been  noted  that  this  method  of  approximation  fails  because  of  the 
failure  of  pointwise  convergence  of  the  sequence  of  approximate  operators 
(lemma  3.4.12  and  theorem  3.4.7)  for  the  same  reason.   In  the  following 
section  discussing  Kantorovich ' s  analysis,  it  will  be  seen  that  while  point- 
wise  convergence  of  the  approximate  operators  is  obtained  in  the  approximation 
space  H(ck')  (condition  I),  the  remaining  Kantorovich  conditions  (section  1.3.2) 
fail  to  be  satisfied  due  to  the  failure  of  the  polygonal  approximations  to 
converge  in  H  (a)  . 

From  lemma  3.4.11  we  cannot  assert  pointwise  convergence  of  the 
approximate  operators  in  H(o/).   Also  the  previous  discussion  on  compactness 
implies  that  the  operator  LI  is  compact  in  11(0/)  if  it  maps  a  bounded  set  in 
H(a)  into  II  (y),  v-W  (an  equi-H(a')  set).   But  the  discussion  of  the  generali- 
zed Privalov  theorem  at  the  beginning  of  section  3.4.3  does  not  imply  a 
relation  of  this  nature. 
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Thus   Gabdulhaev's   approach   is   essentially  to  obtain   "compactness" 

and  pointwise   convergence    in  the   coarser   topology  of  H(|3),   3<a .      Thus   for  x 

restricted   to    lie    in  H(of),   we   obtain     U  x-UxL^  mU  P  x-x  L-*0   as   n-^°   by 

v  "  n    Mp  "  "p"  n   "p  J 

lemma  3.4.11  and  corollary  3.3.4.   Further  with  x  restricted  to  H(a) ,  we 
have  UxeH(<*)  .   Thus  using  lemma  3.4.10  and  corollary  3.3.2,  [ux:||x  ^l}  is 
equi-H(3),  and  hence  U  is  a  "compact"  operator  in  the  space  H(3).   Thus, 
formally,  the  abstract  convergence  theorem  1.3.2  applies  in  the  space  H(3). 
But  convergence  in  the  H(|3)  metric,  as  in  Gabdulhaev's  theorem  3.5.1,  does 
not  imply  convergence  in  the  solution  space,  H(cy)  ,  of  the  equation.   That 
is  the  H(|3)  limit  of  the  sequence  of  approximate  solutions  implies  that  the 
limit  function  belongs  to  H(|3),  but  it  is  not  necessary  that  this  function 
belongs  to  H(cv),  ci>$  .      Hence  the  limit  of  the  sequence  is  not  necessarily  a 
solution  to  the  equation  3.5.1. 

This  situation  may  be  compared  to  the  well-known  fact  that  uniform 
convergence  (in  C)  of  a  sequence  of  functions  does  not  necessarily  imply 
convergence  of  the  sequence  of  derivatives  (i.e.,  convergence  in  the  smaller 
space  C').  More  specifically  compare  Loud's  construction  [30]  where  a  uni- 
formly convergent  (in  C)  sequence  of  functions  in  H(l)  converges  to  a  function 
that  is  in  H(P)  exactly.   It  is  clear  that  the  sequence  also  converges  in 
H(3)  in  this  case  but  not  in  any  smaller  space  including  H(o')  and  H(l).   The 
basic  problem  is  that  the  subset  H(Q')  C  H(|3)  is  not  closed  or  complete  in 
the  H(3)  norm  when  the  Holder  constants  of  the  subset  H(o')  are  not  fixed. 

3.6  General  Theory  of  Approximate  Methods 

In  this  section  Kantorovich's  general  theory  for  equations  of  the 
second  kind  (section  1.3.2)  is  considered  for  the  case  of  a  linear  integral 
equation.   This  approach  is  compared  to  that  of  Anselone  and  the  similarity 
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of  the  Kantorovich  conditions  to  those  of  Anselone  is  demonstrated.   Although 
Gabdulhaev's  convergence  theorem  (theorem  3.5.1)  does  not  directly  depend 
on  the  Kantorovich  convergence  theorem  (theorem  1.3.3),  his  approach  is  based 
on  the  Kantorovich  theory.   It  is  noted  that  an  appeal  to  the  Kantorovich 
theory  does  not  avoid  the  basic  problems  of  the  Gabdulhaev  approximations. 

Kantorovich  projects  the  exact  equation  Nx  -  x-Kx  =  y  in  the  normed 
space  X  into  an  approximate  equation  Nx  -  x-Kx  =  y  in  the  complete  subspace 
X  c  X .   Anselone 's  basic  approach  is  to  obtain  an  approximate  operator,  K  , 
to  the  exact  operator,  K,  both  acting  on  the  complete  space  X.   By  regarding 
the  equation  with  the  approximate  operator,  N  x  =  x-K  x  =  y,  as  the  approxi- 
mate equation,  similar  results  are  ultimately  obtained. 

Consider  the  Fredholm  equation  x-Kx  =  y,  with  continuous  kernel  k, 
in  the  space  X  =  C.  The  usual  algebraic  system  (obtained  by  piecewise  poly- 
nomial approximation  of  the  integrand) 

mn 
x(t  )  -  S   A  k(t  ,t  )x(t  )  =  y(t  )  ,  j=0,l,...,mn        (3.6.1) 
J    i=0      J  J 

will   be    regarded   as   an   approximate    functional   equation    in    the    space  X    of 
bounded   sequences   of    length    (mn+1) .      This   equation   may   be  written    in   symbolic 
form  x-Kx   =  cpy  where  cp   denotes    the    linear   operation  mapping  X    onto  X.      We 
take    as  X   the    space   of   piecewise    polynomial   functions   of   degree   m,    on   the 


system  of   equi-spaced   base    points.      Let  cp      be    the    linear   operation   providing 
a  one-to-one   mapping  of    the    complete    subspace  X   onto    the   complete    space   X, 
and  cp    coincides   with  cp      on  X.      In    this   case   condition    I    (section    1.3.2) 
becomes    condition   la  [25,    p.    558] 

la:    For   every   xr.X  ,||cpKx-Kcp   x||    ^  tj||x||  ,    where   T]  =   Tj||cp      || 
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We  now  show  that  the  conditions  on  the  kernel  which  imply  that  the 
Anselone  hypotheses  are  satisfied  (theorem  2.3.4),  also  satisfy  the  Kantoro- 
vich  conditions  (section  1.3.2). 

Condition  I  basically  expresses  the  convergence  of  the  approximate 
operators  in  the  approximation  space  X,  and  in  fact  shows  uniform  convergence 
of  the  operators  over  a  bounded  set  in  each  approximation  subspace  (compare 
lemma  2.3.1). 
Lemma  3.6.1 


If  the  approximate  operators  are  obtained  by  a  sequence  of  composite, 
interpolatory  quadrature  rules,  then  condition  I  is  satisfied  if  cu  (h)-*0  as 
h-0. 


Proof: 

Considering  the  equivalent  condition  la: 

1  mn 

||cpKx-Kq>  x||    =  max    |J*   k(t    ,t)x(t)dt   -  S   A  k(t    ,t   )x(t   )| 
j^mn   0  J  i=0  J 

1 

£  max  Ik     x   -   pV    xldt  (3.6.2) 

.^       ,J     '    t .  n  t  .    ■  v  ' 

J^nin  o  J  J 

For   an   arbitrary   value    t    in  one   of   the   n   sub intervals    (as    in    lemma  2.2.2) 


Ik  X(t)-Pmk  x(t)|    ^    |  (k    (tO-P^V    (t))x(t)|    +    \^\    (t).x(t)-pmk  x(t)| 
'    s    v  nsv/l         ,vsv/      nsv//v/|         'nsv/       v/      n   s  ' 

m 
Since  x(t)  =  pmx(t)  and  thus  x(t.)  =  x(t.),  and  further  E  I  .(t)  =  1,  then 

n  J       J  j=0  J 

using  the  notation  of  lemma  2.2.2 

m  m  m 

|p  k    (t).x(t)-pmk  x(t)|    =   |   E  k   (t.)A.(t)  S  x(t.)4.(t)    -  2  k   (t.)x(t.)4.(t) 
'nsv/       v/      nsv/|  ±=0   s      L      x        j=0        J      J  i=0   s      i  2/    l 

m 

M(t)i 
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)(k    (t   )-k    (t   ))(x(t   )-x(t   ))i    (t)H    (t) 

Z-jSISJ  J  11  J 


mfl 
where    the    sum   is    taken  over    the  C„    terms    consisting  of   all   pairs   of   values 

i,j    such   that    i^j.      A   crude   bound   sufficient   for   our   purposes    is   given  by 

(compare    proof   of    lemma   3.3.3) 


uo(k    ,mh).2||x||X]|j&.(t)j&  .(t)|    -   2m  u)  (k    ,h)  .||x]|2mC  .™C  . 
v    s        '      1  j  s  "    "         1      j 


2m   .  m 

=  m(JU(k    ,h).||x||(  S        C.    -  S    (mC.)    ) 
S  i=0        L      i=0        L 


=  uo(k    ,h)  .||x||  .m(22n1-(2m)  :/(m!)2) 


Hence    I  k  x  (t)  -p\  x(t)  I    <  uu  (k    ,h)  .||x||  f  mlk  II    +  m(2   m-(2m)!/(m!)    )}    using 
lemma   2.2.2.      Substituting    into    (3.6.2) 


!|9Kx-Kcp   x||    £  ||$||    max  B(m)u)(k      ,h)  .||x||    £  ||  $||  .B  (m)u)    (h)  .||x|| 
j<imn  j 

o  o 

where   B  (m)    =  m||L  ||    +  m(2      -(2m)!/(m!)    ).      Thus    condition   la    is    satisfied  with 

fl  =  ||$||  .B(m).U)    (h). 

The  condition  II  is  related  to  the  compactness  of  the  operator  K. 


(Compare  lemma  2.3.2) 
Lemma  3.6.2 


Condition   II    is    satisfied  with  T]      =  u>    (h)  (m||  L  ||  .i|$||  )  . 
Proof: 

We   may    take    the   element   x    to   be    P  Kx .      Thus 
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|kx-x||    =  Hkx-P^KxII    ^  m||L  ||u)(Kx,h)    by   lemma   2.2.2 

<■  mjJL  I!    ||$l|u)    (h)||x||    by    lemma  2.3.2 

|  —Jl  II  II  Q     \  /    I    I  J 


Finally  by  applying  a  similar  argument  to  the  function  y.  we  find 


> 


P  y  =  yeX   such   that  ||y-y||    ^  m||L  ||uo(y,h),    and  condition  III    is   satisfied  with 

\  =  n||Lj|u>.(y,ii)/l|y||. 

Thus  we  obtain  the  analogue  of  theorem  2.3.4,  showing  the  same 
conditions  on  the  kernel  are  sufficient  for  the  Kantorovich  convergence 
theorem  (theorem  1.3.3). 
Theorem  3.6.3 

If  the  integral  operator,  K,  on  C  is  approximated  by  a  composite, 
interpolatory  quadrature  rule,  then  it  is  sufficient  for  the  application  of 
the  Kantorovich  convergence  theorem  in  the  space  C  that  the  two  partial 
moduli  of  continuity  of  the  kernel  exist  as  simple  moduli  of  continuity. 
Proof; 

By  lemmas  3.6.1  and  3.6.2  conditions  I,  II,  III  are  satisfied  for 
all  n.   Since  the  equation  in  C  implies  yeC,  the  hypotheses  on  the  partial 
moduli  of  continuity  imply  Tj ,  T)  ,  T)  ->  0  as  rr*00 .   Since  K  has  a  finite- 
dimensional  range  it  is  a  compact  operator  and  hence  the  second  hypothesis 
of  theorem  1.3.3  is  satisfied.   Hence  theorem  1.3.3  applies  if  the  equation 
has  a  unique  solution. 

As  in  section  2.3,  the  continuity  of  the  kernel  k  implies  the 
required  conditions  on  the  partial  moduli  of  continuity.   Hence,  by 
theorem  1.3.3,  the  assumption  of  the  existence  of  (I-K)    implies  the 

algebraic  system  (3.6.1)  is  soluble  for  sufficiently  large  n  (as  in  the 

- 1-* 
Anselone  theory),  and  the  approximate  solutions  cp  x   (the  piecewise 
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polynomial    functions   obtained   from   the    solution  of    (3.6.1))    converge    to   the 
exact   solution. 

The    condition   III    is   not   generally   required   in  Anselone's    theory 
since    the    approximate   equation  has    the    same    right    side    as    the    given   equation, 
However   the   condition    is    in   fact   used  by  Anselone    and  Gonzalez-Fernandez 

[6,    p.    256]     in   the    case    of    the    iterated  operator   equation  where    the   new 

p-1 

right   side   would  become   u   =    (I+K+...+K        )y.      If   u    is   obtained  exactly   there 

is   no   problem  but   ordinarily  u  would   also   be    calculated  by   approximate    inte- 
gration,   yielding  an   approximate    value   u    . 

In   the    case    of   Anselone's    theory   the   approximate   equation    is   also 
considered    in   the    space  X,    so   that  X   =  X   and   conditions   II   and   III   are 
trivially   satisfied.      Condition   I    becomes    the    condition   that  K  and  K   are 
"close    in   the    space  X":      K-kI!"^  .      If   K     were    an    integral   operator  with   a 

r  ii  n 

kernel  k    (s,t)    such    that   k    (s ,  t)—k(s  ,  t)    uniformly   as   n-*° ,    then     K   -K  rO.    and 
nv  n  n 

the    standard   theory,    based   on  operator   norm  convergence,    for   the    convergence 
of    the    approximate    solutions   would   follow    (Proposition    1.3.1). 

However   operator   norm  convergence   does   not  hold  where    the   operator 
approximation    is   based   on  a  numerical   quadrature    formula   since    in    this   case 
[l,3]    ||  K   -k||    2;  ||k||,    n£l.      In   the    space  X   condition   la  becomes    (compare 
(3.6.2)) 


|cpKx-K£px||    ^   max    |     I  k      x    -    P  k      x   dt 
IT         ^    "  .-J't.  nt.' 

J      0  j  j 


Lc(u)t(h)-hu(x,h).||k||/||x||)]||x||         by    (2.3.1) 


The    term   in  brackets   may   be    taken   as   T] ,    but    in   this    case  Tj    clearly   depends 
on  x.      Thus,    as  with   Anselone's    theory,   we   have   only   pointwise    convergence. 
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of   the    approximate   operators   on  X.      Arise  lone    compensates    for   this    inadequacy 
by   showing   that   compactness    of    the    operator  gives   uniform  convergence    for   a 
certain   approximation  of    the    iterated  operator. 

The    primary  difference   between   the    two    theories    is    in   the    inter- 
polation  of    the   discrete    algebraic    system   (3.6.1)    to   form  an  approximate 
solution.      Kantorovich    interpolates    the    solution   of    (3.6.1)    into   the 
approximate   space  X,   while   Anselone   uses    the    approximate   equation   to    inter- 
polate   it    into    the   space  X.      However,    we   have    shown   that   for    integral  equations 
in   the    space   C,    the    same    conditions   are    sufficient    to   apply   either    theory 
(Theorem  2.3.4  and  Theorem  3.6.3). 

Finally,    considering  Gabdulhaev's  method  again,    from  the   point  of 

view  of    the   Kantorovich   theory,    it    is   noted   that   condition  I    (equivalently 

la)    is   obviously  satisfied  with  T]=0  since   the   operator  K   is   taken  to  be   PK. 

Thus    the   approximate   operators   converge    in   the   approximate   space  H(|3),   but 

as   noted   previously,    this    does   not    imply   convergence    in  H(|3).      We    still 

require  TL||p||,  T]J|p|Ho  for  the   convergence  of   the   approximate   solutions    in 

Theorem   1.3.3.      The  Kantorovich   theory   states    that    these   conditions   are 

satisfied    in  this   case    if   P  -*I  on  the   complete   space  H($)    and   the   operator  U 

is    compact  [25,    Theorem  6,    p.    555]  .      But    it  was    shown    (lemma  3.3.5)    that   it 

is  not    true    that   P  —I  on  H(B).      For   condition  II  we  have 
n  v   ' 


|Ux-x  L    =  ||Ux-PUx  L   ^     3sup  uu(Ux,6)/6      by   lemma  3.3.3  with  m=l, 
P  P        (K6^h 


and  xeH(|3)    implies  UxeH(f3).      We   can  assert  only   that   the   right   side    is 
bounded  and  not   that    it   converges.      Under   Gabdulhaev's    conditions    (3.5.1)    and 
the   assumption     xeH(8),    0<P<6<tf ,    then  UxeH(6)    and   the   right   side    is  bounded 

C        Q 

by  3llu|L   M2[x;6]h        .      Hence   condition  II    is   formally  satisfied  with  a  term 
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similar    to    that    in   the    convergence    condition    (3.5.4).      That    is   Gabdulhaev 
has   effectively   assumed   a   "compactness"    in   condition  for   the   operator  U , 
but   as   noted    in   section   3.5.1,    this    is   not   valid   since   we    are   not    justified 
in  restricting  x   to    the    space   H(6). 

It    is    concluded   that  neither  Kantorovich's    theory,    nor  Anselone's 
theory,    may   be    applied    to    justify   Gabdulhaev's   method  of   direct   approximation 
of   the    singular    integral  equation. 

3.7      Concluding  Remarks 

It  has   been  shown  by   the    analysis   of   Gabdulhaev's    paper   that   an 
attempt   to   ensure   uniform  convergence   by   considering    the   metric   of   a   coarser 
topology,   H((3),    for    the    space   of    functions   H(cv)  ,    (0<a)  ,    fails   basically 
because    the    space  with    this   metric    is   no    longer   complete.      Equivalently   an 
attempt    to   seek   a   solution    in   a   "smoother"   space   of   functions    fails    for   the 
same    reason.      The    introduction   of    a   sufficiently   smooth   function  k    into    the 
kernel,    k(s ,  t)  /  (t-s)  ,    does   not   help   since    for  xeH((3),    Kx   still  belongs    to 
H(3)    (section  3-4.3).      This    is   not   unexpected   since    an   operator   of    this    form 
readily   reduces    to    the    sum  of   a  Fredholm  operator   and   a  Hilbert   operator. 

The    only   possibility   that   would  be   of    any   assistance    in   this 
analysis    is    to   ensure    that    the   kernel    is    absolutely    integrable.      But    this 
essentially    implies    that   k(s,t)    has    a   zero   on   the    diagonal   s=t,    and    thus 
reduces    to    the   case   of   a  Fredholm  operator  with   a  weakly   singular   kernel. 
This    case   was    considered    in   section   2.5. 

It   was   noted    in   Chapter   2    that    for   a  Fredholm  equation  with   given 
functions    in  C,    every  bounded,    integrable    solution   also    lies    in   C.      Thus   we 
I  may   seek   approximate    solutions   which   converge   uniformly    (in    the    sup   norm)    to 
the   exact   solution.      If    the    right   side   of    the   equation    is    continuously 
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.ies    in   the 


.   ,_-,  A    nWs>t)|dt  =   e'(s)eC,    then  the   equation   1 

differentiable   and  J  \      ^ ^|aL        &   ^    > 

p   r'        A  similar   argument   shows    that  every  bounded,    integrable 
smoother   space  C    .      A  simiwi        6 

•  1l7    11p-    in   the    required   space   C'.      Thus    it   is 
solution  of    this   equation  necessarily    Ires    m  req 

again   sufficient   to    require    convergence   of    the    approximate   solutions    in   the    sup 

norm,    providing   the    approximate   operators   converge    pointwise    in    the   C    norm. 

f^k(s,t)N  .r><0<1'i    is 
That    is,   by  analogy  with   le^a  2.3.1,    providing   the   set  K^T" V° 

i  •         .-r,    t-v.^   ra<;p   of   a  Fredholm  equation 
equicontinuous.      The   same   argument   applxes    rn  the   ease 

in   the    space  H<f»    (If    f|Ms,t)|dt  -  ,(.)«»».      »   follows   from  lemma  3.4.3 

that  approximate  solntions  converging  in  the  sop  norm  also  converge  in  the  H(P) 

norm  providing  the  approximate  operators  converge  in  H(p)  ,  pointwise  for  xeR 

(i.e.,  if  [k  ),  {kt}  are  both  equi-H(p))- 

Hclver  this  argument  cannot  be  applied  to  a  singular  operator 
since  the  principal  value  integral  is  not  absolutely  integrable.   It  was       | 
noted  in  Section  3.4.3  that  it  is  necessary  that  the  operand  also  satisfies 

■   i  <„»  t-hat  a  direct  approximation  requires  con- 
some  Holder  condition,  implying  that  a  direct  PP  . 

vergence  in  the  Holder  space. 

i  a»a    t-W  the  Hilbert  operator  is  not  sufficiently 
Thus  it  is  concluded  that  cne  tiiiucit  v 

"smoothing.,  to  permit  direct  application  of  the  uniform  approximation  theory. 

This  would  include  a  scheme  involving  perturbation  of  the  Kernels  (compare 
possibilities  in  Titchmarsh  [45])  since  the  required  properties  ultimately 
depend  on  the  compactness  of  the  original  operator  (compare  Anselone  [4]). 
Similarly  a  sequence  of  operator  differences  still  does  not  avoid  the  basic 
problem  that  the  Cauchy  principal  value  integral  is  only  conditionally  and 

not  absolutely  convergent. 

Hence  we  finally  conclude  that  the  only  practical  way  to  deal 
numerically  with  a  singular  integral  equation  is  to  consider  some  semi-indirect 
method  in  which  the  singular  equation  is  first  reduced  to  the  more  tractable 
case  of  an  equation  with  a  Fredholm  operator. 
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4.   UNIFORM  APPROXIMATION  OF  THE  SOLUTION  OF 
SINGULAR  INTEGRAL  EQUATIONS  IN  A  HOLDER  SPACE 


4.1  Introduction 

It  has   been   shown   in   the    preceding   chapters    that  a   singular    integral 
equation    is   not   amenable    to   a   direct   method  of   uniform  approximation.      The 
method   proposed    in   this   chapter   might  be   classed   as    a   semi-direct  method 
since    the   equation   is   first   reduced   analytically    to    the   case   of   a  Fredholm 
equation.      In   principle    the    latter   equation   can  be   dealt  with  directly   accord- 
ing  to    the    theory    in  Chapter   2.      However,    a   practical   approximation   scheme 
involving  only   the    given   functions    (section  4.3)    requires   some   additional 
analysis    to   show   that   the    approximate   solution   of    the   equivalent  Fredholm 
equation   does    in   fact  converge   uniformly    to    the   desired   solution. 

It    is    shown   that    it    is   necessary   to   strengthen   the   hypothesis   on 
the   density   function   k(s,t)    so   that    this    function    is    required   to   satisfy   a 
Holder   condition  with  exponent  v>u.      with   respect    to  both    variables,  (instead 
of   exponent  \j»    as    in  Muskhelishvili   and   Pogorzelski)    in  order    to   obtain   a 
consistent   equation    in   the   space   H(|j,).      The    forms   of    the    reduced  Fredholm 
equations   are   given   and   an   analysis   of    the   convergence  of   an  approximation 
to   the  weakly   singular   equation   obtained  by  Noether's   method    is    given. 

4.2  The    Singular   Integral  Equation 

To   avoid   the    complications  at    the   end-points   of   a  non-closed   arc, 

we   consider    the    range   of    integration  to   be    a   smooth   closed   Jordan  contour  L 

(compare   [33,    Chapter   6]  )  .      Consider  the    linear    singular    integral   equation 
of    the   normal  1.38]     type; 

Mp(s)         (aI+K)cp(s)         a(s)cp(s)    +  -j  J   k(s't)CP  ^dt   =   f(s)  (4.2.1) 

"l 
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where  a,k,f  are  given  functions  on  L  belonging  to  some  Holder  class  H(y,). 

Muskhelishvili  [33,    p.    114]     seeks    the    solution  of    this   equation 
which   also   satisfies    the    "H   condition";    that    is,    the   Holder   condition   for 
some   unspecified   exponent.      It    is    implied    indirectly  [33,    p.    132]     that    the 
solution  belongs    to  H(p,'),    (Ky, '<jj,    (compare   also  [17,    p.    218]    and   section 
3.5.1).      The   solution  of   the  Hilbert  problem  belongs    to   the   same  Holder  class 
as    the   given  functions  [38,    p.    466]    and   thus   the   solution  of    the   equivalent 
dominant   equation  [38]    belongs    to  H(|J,)    if   the   given   functions   belong   to 
H(U.),  H<1. 

For   the   complete    singular   integral  equation   (4.2.1)    the   solution 
space    is    less   clear.      Despite   the    fact   that    some    authors    (Dolph  [15,    p.    152]) 
apparently    infer    that    the    preservation  of    the   Holder   class,    for   a   density 
function  of   two   variables,    follows   as   an   immediate   consequence   of   the   Privalov  ! 
theorem,    the   generalized   Privalov   theorem  [38,    p.    443]     implies   a  weaker   result' 
(see    discussion   at   the   beginning   of    section  3.4.3).      Thus   K:  H(p,  )-*H(|jl  )    if 
keH(v}|i),  v>u, ,   but  not   if  keH(|i). 

The    theoretical   investigation  and  solution  of   the   complete    linear 
singular    integral   equation    is    obtained  by   reducing   the   equation   to   an 
equivalent  Fredholm   (or   regular)    equation.      Any  of   the    three   methods   of 
regularization    (or   reduction)    in  use   [18,    Ch.    3]    ultimately    involve    the 
product  of    the    singular   operator  with    itself.      The   expression   for    the    product 
operator,   using   the   Poincare-Bertrand   transformation    (section  3.1)    involves 
a  kernel   of    the    form 

pk(s,u)k(u,t)du   =  J_    f  r.   k(s,u)k(u,t)         _    p   k(s?u)k(u,t)du| 
J     (u-s)  (t-u)  t-s    U  u-s  ■-'  u-t  J 
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By   the    generalized   Privalov   theorem,    the    first    integral   on   the    right   belongs 
to  H(|i.,u,')    and   the    second    integral  belongs    to  H(v',u.')    where  [X '  ,   v '    are 
arbitrary   positive    constants    less    than  |i,  ,v    respectively.      The   resulting 
density   function   thus    satisfies    the  H(u.,p')    class,    choos  ing  p,  '<jj,<v  '<v  . 
Hence    the    corresponding   product   operator  maps   H(u.)~*H(p  ' )    (by   the    generalized 
Privalov   theorem). 

Hence    in   order   that   the    reduced   equation    is    also   an  equation    in 
the    space  H(p)    we   require    that    the   density   function,    k,    satisfy  keH(u)  , 
v>|i .      With   this   change    in  hypothesis    in   equation    (4.2.1)    we   have   an  equation 
in   the    space   H((jl)    and  may   seek   a   solution   in   this    space. 

The    following   lemma   concerning   the    associated  weakly   singular 
operator  L    is    also  noted. 
Lemma  4.2.1  [38,    lemma   2,    p.    495] 

If   k(s,t)eH(v)    then   the   weakly   singular   operator,    L,   with   kernel 
4(s,t)    =    (k(s ,t) -k(s ,s) )/ (t-s) ,    transforms   each   bounded,    integrable    function, 
cp,    into  the    space   H(p)    for  p,<v  . 

Actually   Pogorzelski ' s    result   states    that    the    transformed   function 
satisfies    the    condition  with   exponent  v/2.      But   an   analysis   of    this    proof 
and   the   corresponding   result   of  Muskhe  lishvili  [33,    p.    135]     reveals    that    in 
effect    it    is    shown    that   Lcpcll  (> )    where   X    =  min(v  -a  ,  l~a)  ,    where  oi    is   an   arbi- 
trary  positive   number    less    than  v.      Thus    since    (Kv^l,    setting  a    =  v-p.    gives 
the   desired   result. 

4.2.1     Reduction   of    the   Singular   Integral  Equation 

The    theory   of    reducing   operators   may  be    found    in  Muskhelishvili 
1  33]     and   Pogorzelski  [38]  . 

The   operator   adjoint    to    the    dominant   operator,   M    ,    is   defined  by 
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M°V(s)    s  a(s)\|r(s)    -  ^r  J  b(^^t)dt      ,   where   b(t)    =  k(t,t).      (4.2.2) 

Li 

i 
Noether's  method  uses   the   operator  M        as   the   reducing  operator. 

Equation   (4.2.1)    is   reduced   to   the  weakly  singular  Fredholm  equation 


:  »p(s)     :   (ats)-b(s))cp(s)   ■!   ~|"  [  I    ;y^):,i,;(duicP(t)dt  =  M"   f 


TT       Jj        j_i 


(4.2.3) 
The   kernel  of   the    integral   operator  may  be  written 

n,s    t)    =    P    b(u)k(u,t)du   m  _l_   ,  p    b(u)k(u,t)du    _    j.    b(u)k(u,t)du]    =   q(s,t)-q(t,t 

v  J     (u-s)(t-u)  t-s     J  u-s  J  u-t  t-s 

Lj  JLf  J_j 

where  q(s , t)eH(v ,v ')  by  the  generalized  Privalov  theorem. 

i 
Although  solutions  of  the  Fredholm  equation  are  normally  sought  in  ' 

the  class  of  continuous  functions,  there  is  no  necessity  to  impose  before-   • 

hand  the  additional  condition  that  the  solution  of  (4.2.3)  belongs  to  the 

space  H(p>)  since,  by  lemma  4.2.1,  each  bounded,  integrable  solution  of 

(4.2.3)  necessarily  belongs  to  H(|i). 

The  index,  K,  of  M  is  defined  by  the  integer  K  =  - — r[  log — p3  _ 

J  °  2tti        °a+b  L 

where   the   symbol  []       stands    for   the    increment,    suffered  by   the    function   in 

■Li 

braces,    on  one   circuit  of  L   in   the   positive   direction.      The   theory  of 
Noether's   method  shows    that   in   the   case    that   the    index,  K,   of  M   is   non- 
negative,    the    reduced   equation    is   equivalent   to   the   original  equation.      If 
K<0   then    it    is    impossible    to   find  a  reducing  operator  of   this    type. 

However,   Vekua's   theorem  of  equivalence  [38,    p.    512]    shows   how  to 
find  an  equivalent  Fredholm  equation   in   the   case    the    index   is   negative. 
Introduce    the   new  unknown   function,   f ,    into    the    singular   equation  by  the 
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.substitution  cp(s)    =M     t  (s)  .      Then,    for  K<0,    the   equation    (4.2.1)    is    reduced 
to   the  weakly   singular  Fredholm  equation 


M(M°V)(s)    =    (a(s)-b(s))t(s)    +  W    (u-s)(t-u)du^(t)dt   =   f(s)  (4-2-4) 

J-i       ij 

As   before    the   determination   of    the    solutions,  cp ,    of    the   original 
equation    (4.2.1)    is   equivalent   to    the    determination   of    the    solutions,   \|f  ,    of 
(4.2.4),   where  cp    is    obtained  by   the    additional    step   of   evaluating  cp   =  M     ty  . 

Hence,    after    the    determination   of    the    sign  of    the    index   of    the 
equation,    the    singular    integral   equation  may  be    approximated  by   considering 
the   appropriate   form   (4.2.3)    or    (4.2.4). 

The    singular    integral  equation  may  not   have    a  unique    solution    (the 
number   of    linearly    independent   solutions    is    related    to    the    index) .      The 
approximation    theory    (section    1.3.2)    assumes    that   the   given  equation  has    a 
unique    solution.      It    is   noted    that    the    problem  of   obtaining  a   particular  solu- 
tion of   a  Fredholm  equation  with   a   non-unique    solution  has   been  considered  by 
Atkinson  [8]  ,   who  obtains    a   selective    inverse.      A   similar   approach  may 
yield   particular   solutions    in    the    case    of    a   singular    integral   equation,    but 
this   problem  was   not    investigated. 


4.3      Approximation   and  Convergence 

The    preliminary   analysis    reducing   the    singular   equation   to   a 
Fredholm  equation  before    the   application   of   approximate   methods    seems    to   be 
an  unfortunate   but   necessary    condition.      On   the   other  hand    this    reduction 
does   allow   the   approximation   to   be    considered    in    the   space   C  which  has    a 
much   simpler   norm   than   the   Holder   norms.      Thus    the   complications   of    the 
determination  of   the    appropriate   Holder   space   and  calculations   with    the 
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relatively  unwieldy  Holder  norms  are  avoided.   Further  for  a  given  singular 
equation,  the  appropriate  values  may  be  substituted  directly  into  the 
approximate  equations  obtained  in  this  section  and  so  for  practical  purposes 
the  details  of  the  reduction  process  may  be  ignored.   The  only  preliminary 
calculation  that  is  necessary  is  the  determination  of  the  sign  of  the  index 
of  the  equation  in  order  to  decide  which  of  the  forms  ((4.2.3)  or  (4.2.4)) 
of  the  Fredholm  equation  it  is  necessary  to  consider.   It  is  noted  that 
Gakhov  [18,  p.  89]  gives  a  computational  method  of  determining  the  index. 
For  the  purposes  of  this  section  we  suppose  a  singular  integral 
equation  (4.2.1)  is  given,  that  it  has  a  unique  solution,  and  that  it  has  a 
non-negative  index,  so  that  its  approximate  equation  will  be  determined  from 
equation  (4.2.3).   Equation  (4.2.3)  may  be  written  (compare  section  4.2.1) 


Np  =  [  (a2-b2)I  +  ~Q]cp  =  g  (4.3.1) 

TT 


where   I    is    the    identity  operator, 


and 


»(•)  ■J1<'^(t't><Kt)dt, 


g(s)    =  M°'f(s)    =   a(s)f(s)    -  -It    f   b<t>f<t)dt 

v  x  TTl    -J  t-S 


Since    the   kernel  of    the   operator  Q    is  weakly  singular   an   immediate   approach 
might   utilize    the    approximate    product    integration   technique   of  Atkinson  [  16J  . 
However    in   addition   to   the   difficulty   of   determining   a   suitable    split    in   the 
integrand    it    is   noted    that   the    functions    q   are    themselves   expressed   as    princ 
pal   value    integrals.      Thus,    in   general,    we   will   also  have   q   replaced   by   some 
approximation   q 
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in  order  to  simplify  the  arguments,  we  will  assume  here  that  the 
operator  P*  represents  piecewise  linear  interpolation  with  respect  to  the 
variable  t,  but,  as  noted  in  chapter  2,  interpolation  by  higher  order  poly- 
nomials involves  only  a  change  in  the  constants  and  not  in  the  substance  of 

the  argument.   Define  an  approximate  operator  0  by 

xn  y 

„  ,  ,    -  p?Cq(a,t)-q(t,t)<p(t)] 

Qn^s>=JL   ~s dt  (4.3.2) 

and  an  approximate  operator  0n  by 

n  J 

^    }        \  t-7~" dt  (4.3.3) 

where  qn  is  some  approximation  to  q. 

Defining  the  approximate  operator  N  =    (a2-b2U  +  — nn   ,'t-  *. 

n  '  2    n  ' 

necessary   to   obtain   pointwise    convergence    in    the    space    C,    oAhe   operators 
Nn    (and   thus   of   (£)    for    the    application   of    the    abstract    theory.      Thus    it    is 
required    to   show   that  Q^Q    pointwise,    and,    since    this    result  depends    on   the 
Holder   properties    of   q,    q^q   uniformly    in   the   Holder   norm  of   H(H)  . 
Define    the    approximation   q    (s,t)    by 


qn(S'fc)    =  J   ~s du  (4.3.4) 

JL_i 


Lemma  4.3. 1 


Let   qn   be    defined   by    (4.3.4).      Then   q^q   as   rr~ ,    uniformly  with 
respect   to   s,t    in   the    Hflj,)    norm  where   k(u,t)€H(V),   v>Ji . 
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Proof: 

Following   the   analysis    in  Chapter   3   and   the   generalized   Privalov 
theorem  we   have 

M1[qn(s,t)-q(s,t)]     <:  n   M1[p^b(u)kt(u)-b(u)kt(u)] 

+  C'M2[(P^b(u)kt(u)-b(u)kt(u));MJ 
£  TT  M2[bkt;v]hV   +  2c'M2[bkt;v]hV"^ 

using   the   relations    (3.3.1)    and    (3.3.7)    for  m=l.      Note    that   these  Holder 
constants   exist   finitely   since   b(u)    =  k(u,u)eH(v)  ,    and   that  h  =  0(— )  . 

Let  i   be   the   finite    length  of   the   arc  L  and  v1    an  arbitrary  value 
satisfying  p-<v  '<v  .      Then  from  the   generalized  Privalov  theorem  again 


M2[(qn(s,t)-q(s,t));M.]    ^  C^C  (P^b(u)kt(u) -b(u)kfc  (u)  )  ;uJ 


+     sup  {^[(P^^k^-b^k^u))^']  |trt2|V 


1  t1"t2l<  ' 


I  v'-ij, 


+  C3M2[(P^b(u)kt(u)-b(u)kt(u));v']|t1-t2| 

•log|t1-t2| 


where   C9,C^    are    constants.      Then   again   applying   the   estimate 


M2(P  x-x;|4.)   ^  2M2(x;v)hV_Mj    for  xeH(v)  ,    (relation    (3.3.7)) 


we   obtain    the    inequality 
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i 


M2E  (qn(s,t)-q(s,t));|a,]  <:  C^t  bkt  ;v]  hV_V  ,  C^   constant. 


Combining  these  results 


||qn(s,t)-q(s,t)||^  =  M1[qn(s,t)-q(s,t)]  +  KJ_  (qn(s  ,  t)  -q(s  ,  t)  )  ;pj 

£  C5  M2[bkt;v]hV"V  (4.3.5) 

Since   C,-    is    a  constant,    the  Holder   constants    are   uniform   in  u   and   t,    and 
h  =  0(— )  ,    the    relation    (4.3-5)    gives    the   required   convergence    as   rr100 . 

A   similar   result   of   course   holds    for    the    function   q(t,t). 
Lemma  4.3.2 

For   each  fixed  cpeH(p.)    the    functions   Q  cp(s)    obtained   from   (4.3.2) 
converge    to  Ccp(s)    in   the   norm  of   the    space   C. 
Proof: 

Let   the    sup  norm  be    denoted   by    the    subscript   C.      Let  |i  '    be    a 
positive   constant   satisfying  0<\x  '<ji  .      Then  by   an  argument    similar   to    that   used 
in  obtaining    (4.3.5)    in    the    previous    lemma  we   obtain 


(Qn-Q)<pHc   -  M1[(Qn-Q)cp] 

=S    B  Wp^(q(s,t)-q(t,t))cp(t)-(q(s,t)-q(t,t))cp(t)] 


+  M2[(P*(q(s,t)-q(t,t))cp(t)-(q(s,t)-q(t,t))cp(t));p/] 


where  B    is   a  constant.      Since   qeH(v')c  H(|i)    and  cpeH(M-)  ,    the  usual  estimates 
on   these    terms   give 
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(Qn-Q)<P!IG  *  B1.M2[(q(s,t)-q(t,t))cp(t)  ;n]  tf  "^  '        (4.3.6) 


and  the  required  convergence  is  obtained  as  rr*=°   since  h  =  0(— )  ,  B   is  a 
constant,  and  the  Holder  constant  is  uniform  in  s  and  t. 


The  results  of  these  two  lemmas  are  now  combined  to  obtain  the 


required  pointwise  convergence  for  the  approximate  operators  Q  . 
Lemma  4.3.3 

For  each  f ixed  cpeH(|i)  the  functions  Qcp(s)  obtained  from  (4.3.3) 
converge  to  Qcp(s)  in  the  norm  of  the  space  C. 
Proof: 

From  the  identity  Q  cp -Qcp  =  Q  cp-Q  cp  +  Q  cp -Qcp  we  obtain 


iQ>-wllcsll(Q;-QnHlc  +  l!Qn9-*llc 


The  second  term  converges  by  lemma  4.3.2.   Arguing  as  in  the  previous  lemma 


Wn-V^C^Wn-V^ 


tr 


£   B||p^[((qn(s,t)-qn(t,t))-(q(s,t)-q(t,t)))cp(t)]l|^ 


£    B.Hp'II    llcpll    (llq    (s,t)-q(s,t)||    +11  q    (t ,  t) -q(t ,  t)l|    }  (4.3.7) 

ii    n"|j.l|Tll|j,   "nnv  V    l|nnv        '    nv        ' 'p.  v 


Ijcpjl      is   a  constant   for   a  fixed  cpeH(|i)    and   it    is   easily   shown    (compare 

Gabdulhaev  [42,    p.    215])    that  ||p  ||      ^   3.      Thus  we   have   convergence   as   n-00  by 

n  M< 

lemma   4.3.1. 
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Actually  in  view  of  (4.3.7)  and  (4.3.6)  the  convergence  in  the 
above  lemma  is  uniform  over  a  bounded  set  of  functions  cpeH(u.)  but  such  a 
result  is  not  surprising  in  view  of  the  fact  that  Q  is  a  weakly  singular 
Fredholm  operator   and   a  bounded   set    in  H(u.)    is   equicontinuous    in  C. 

By   lemma   4.2.1,    the    set  tQp  rll^p'J  r  ^    l)    is   Holder   continuous    (since 
the   kernel    is   weakly   singular)    and   thus    is    a  bounded,    equicontinuous    set    in 
the   space   C.      That    is   Q    is    a   compact   operator    in   the   space   C.      In   fact  we   can 
make    the    stronger   statement   that  Q    is   a   compact   operator    in   the    space  H(M-). 
Lemma  4.3.4 

The   weakly   singular   operator   Q  defined    in    (4.3.1)    is   compact    in 
the   space  H(u.),  p.<l. 
Proof: 

Let  B   be   a   bounded   set   of    functions    in   H(M-).      By  definition   the 

operator  Q    is   compact    if   the    set  QB    is    precompact   and   this    latter   condition 

holds    if,  from  any  sequence   {(£p   }    of    its    points,  we   can   choose    a   subsequence 

f Qcp      }    which   is   convergent    in   the  H(m-)    norm.      The    functions   of    the    sequence 

"k 
{cp  }  C    B   are   equicontinuous   and  uniformly  bounded   and  hence,    by  Arzela's 


theorem,   we   may    select   a   subsequence   {cp      }    which    is    uniformly   convergent    (in 

nk 
the   usual   sense   of    the   sup  norm).      It   follows    that    for  each  e>0,    there 


exists   N      such    that 
e 


sup|cp      (t)-cp      (t)|    <  e    holds   for    j,k  >  N 
fc  .1  k 

Thus  M.tQp      -Qcp     ]     -  max  |f    q  (S  '  t>"^t ' t}  (cpn    (t)-cp      (t))dt| 
j        nk  s       l  c"s  nj  nk 

*  l^n.'^n  Ifc'^V^  •Tnax.[lt-s|V      ldtl 


£  C1-M2[q    ;v]  .e    for    j,k>N    ,    Cj^   constant.  (4.3.8) 
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and  M2L  (0^n  -Qcpn  )  ;p] 
J  k 


=      sup 

cx|  s1-s2|<je 


+ 


x  rn(q(s1,t)-q(t,t))-(q(s2,t)-q(t,t)) 


t-s. 


(CPn    (t)_Cpn    (t))d 
i  k 


'     (q(.2,t)-q(t,t))fi. „    (t)Kp        W)^    "    ^dt} 

L  j  k  1  2 


By   the   generalized  Privalov  Theorem  and  an  analysis   similar   to   that  used   in 
the   proof  of    lemma  4.3.1,    the   first   term  is   bounded  by 


C2'K  "Vllc-^  VV]  •U1-s2r^(H-log|s1-s2|) 
j       k 


<    C3.M2Lqt;v]  .e 


(4.3.9) 


for   some   constant   C-, ,    and   j,k  >  N    ,   where  v>y>  . 

By  the   same   argument  as   used    in  the   proof  of   the   Privalov  theorem  the   second 
term   is   bounded  by 


VK.^nJc'M2[Vv]'lsrs2 
J        k 

^  C,-.M0[q.  ;v]  .6     for  j,k  >  N 
5  2   t  e 


V-|J, 


for  v<l 


(4.3.10) 


Hence  for  [i<l   we  conclude  from  (4.3.8),  (4.3.9),  (4.3.10)  that  the  sequence 

{ Qcp   }  satisfies  the  Cauchy  criterion  for  convergence  in  the  H (fj, )  norm: 
nk 


101 


ftp      -Op  =  M,[ftp      -ftp     ]    +  M  J  (ftp      -ftp      );ll] 

J  k  J  k  J  k 


^  C,.M0[q    :v]  .e  for   j,k  >  N 

6     2      t  J  e 

Since    the    space  H(|J.)    is    complete    the   Cauchy  criterion    is    sufficient   for    the 

convergence  of   the   sequence   of   points   and  we   conclude   that  { ftp     }    is   con- 

nk 

vergent  in  the  H(p.)  norm. 

Note  that  the  result  may  also  be  obtained  using  corollary  3.3.2  since 

the  above  proof  also  shows  that  the  operator  Q  transforms  a  bounded  set  into  a 

bounded  set  in  H(v')  which  is  compact  in  H(jjl  )  .   The  result  may  also  be  obtained 

more  directly  by  expressing  the  kernel  in  the  form  -" — J — 1— li — i—L   =  q  ^  '  '      where 

t-s  I  .       I  l-oi 

a   is   an  arbitrary   positive   number    less    than  v'    and   q    (s,t)eH(v'-cy)   L  38 ,    lemma    1, 
p.    492]    and  using   lemma  4.2.1.      However    this    approach   is   not   amenable    to   ob- 
taining  a   similar   result   for    the   collective   compactness   of    the   operator   approxi- 
mations  as    is   suggested  by   theorem  2.3.3.      Note    that   the   proof  of   the   above 
lemma  gives: 

sup  U)(q    ,6)/8     -»  0  as  6-*0   implies  Q   is   compact   in  H(p,),    analogous   to   the 
'result  of   theorem  2.3.2,    and   thus   we   may  expect   the    same   condition   to   also 


;ive  {Q    :n^l}    is   collectively   compact. 
i,emma  4.3.  5 

The  approximate   operators   [q   :n^l}    defined  by   (4.3.3)    are   collectively 
ompact    in  the   space  Hfli),  \±<l 
roof; 

This    set    is    collectively   compact  by  definition    if   the    set 
^nB:n^l]    is   compact.      Proceeding   as    in    the    proof   of    the   preceding   lemma, 


{cp     }    be   a  uniformly  convergent    (in  C)    subsequence   obtained   from   the    set 
k 


by  Arzela's    theorem.      Thus 


i 
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M..[QncD      -Qncp      ] 


=  max 
s 


P  _[(q_(s,t)-q  n(t,t))(cp       (t)-cp      (t))]dt 
n       n  n  m  •  >\ 

L~~  t-s 


^  lip  |L.||cp      -cp     |L.M0[q    .;v]  .maxf'l  t-sl    "    Idt' 


^   C,.M0[q   „;v]  .e         for    j,k>M    ,    C,    constant, 
1     2   Mnt  J  e         1 


Note    that    IP        =1   for   linear    interpolation  but    the    argument  would   not  be 

n" 

■f-V» 

affected  if  P  was  piecewise  interpolation  by  m   degree  polynomials  since 

then  IP   =  L   where  L   is  the  Lagrange  operator.   Similarly,  following 

11    n"         "    m"  m  °        °        r  J  b 

the   proof  of   lemma  4.3.4, 


M0[  (Qncp      -Qncp      );jji]    £   C0M0[q      ;v]  .« 
2    vxnTm.      nTm.  ^" 

J  k 


'2112U4n. 


for  j,k  >  M 


Hence  the  result  will  be  established  if  we  can  bound  M„[ q  ^;v]  uniformly 

2  nnt 

with  respect  to  n. 


q  (Sl,t)-q  (s0,t)  =  f  P^^kfu.t)^ — ) 

Hnv  1  '    Mnv  2      '        JT  n  v  /    v  '  /vu-s,    u-s„ 


du 


=  n(Pnbkt(Sl)-Pnbkt(s2))  +  $  (srt)-$(s2,t)     (4.3.11) 


where 


P  bk  (u)-P  bk  (s) 

(s,t)  =J  du 

L 


u-s 


and  P  bk  (s)  represents  the  piecewise  linear  interpolation  of  the  function 
bk  ,  evaluated  at  the  point  s. 
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Suppose   xeH(|-0.      Then   P  x£H(l)  c   E([i)  .      Thus 


|p  x(Sl)-P  x(s0)|    ^  U)(P  xjs.-sj)   £  M0[p  x;|i]  -  j  s  n  -  s  0  |  ^ 
1    n        1        n        2    '  n  1      2''  2      n  '12' 


Further  M  J  P  x;pj    ^  M  J  x  ;pj    +  M0[  P  x-x;^u]    ^   3M0[x;^i] 
2      n  2  2      n  2 


Thus  |p  x(s1)-P  x(s0)|    s   3  Mix^.ls.-s.P. 

1    n    v    1        n    v    2/  '  2  '     1     2 

Applying  this  result  to  (4.3.11)  and  applying  techniques  similar  to  those  in 
the  proof  of  the  Privalov  theorem  in  the  estimation  of  the  integral  term,  $ , 
we   obtain    (compare    lemma   3.4.4) 

|qn(s1,t)-qn(s2,t)|    *   sup{  3ttM2[  bkt  ;v]  |  s  x"  s2  '  ^   +  3   c2M2^bkt'V-'  'lsl"s?l^ 

£   C3   M2[k(s,s)k(s,t)  ;v]  |  Sl-s2|V 

Thus  ,   since   keH(v)  , 


M 


2C  q   t'»vl    s   2C3M,[k]M2[k;v]       ,    independent   of   n. 


Hence 


~n  «n       ii  r    n  n       i  r  ,   n  n       .      -i 

Q  cp      -Q  cp  =  M,L  Q  cp      -Q  cp         +  M0L  (Q  9      -Q  9      )  ;m<J 

xnTm  .    xnTmJ'p,  1   xn^m  .    xnTra,  2    vxnTm.    xnTm, 


C.-M^k]  .M2[k;vJ  .e  for    j,k> 


M 


e 


md    it    follows    that    the    set   {q    :n^l}     is    collectively   compact    since 
Qncp:||cp||   ^l.n^l)     is    compact    in  H(\±)   by    the    previous    lemma. 


104 


Since  compactness  in  H(|-i)  implies  compactness  in  the  coarser 
topology  of  the  space  C,  the  lemmas  4.3.3,  4.3.4  and  4.3.5  show  that  the 
hypotheses  of  the  convergence  theorem  1.3.2  are  satisfied  and  thus  the 
solution,  cp  ,  of  the  approximate  equation,  N  cp   =  g,  converges  to  the  exact 
solution,  cp ,  of  equation  (4.3.1)  in  the  norm  of  the  space  C.   In  view  of 
lemma  4.2.1,  this  limit  of  this  sequence  of  approximate  solutions  satisfies 
the  required  Holder  condition  and  thus  is  a  solution  of  the  equivalent 
singular  integral  equation  (4.2.1). 

However  because  of  the  nature  of  the  right  side,  g,  it  would  be 
more  practical  to  consider  the  approximate  equation 


N  cp   =  g  (4.3.12) 

nTn   &n  v      ' 


where    g      is    some    approximation   to   the    right   side,    g   =  M     f,    of    (4-3.1).      Let 
the   approximation   g     be   obtained    in   the   usual  way   of   piecewise    linear    inter- 


polation of   the   numerator   of   the    integrand. 
Lemma  4.3.6 


The   approximation   g      converges   uniformly   to   g. 


Proof: 


|i  ,|  ii     ,        1      P   Pnbf (t)  .„  -    .1      (•  b(t)f(t)Jj| 

g   -g        =      af r       ■ dt    -    af   H — v—  ■< — ^— Mt 

11  &n   51lc        "  tii  .'L       t-s  TTl  JL       t-s  "c 

£  b||p  bf(t)-b(t)f(t)||  /2    (compare  lemma  4.3.3) 
£  3B  MJ  bf  ;ii]  h^ /      ->  0  as  n-*° ,    since  beH(v)  <=  H(p.)  and 
feH(u.)  implies  M?[bf;(i,]  is  a  constant. 

The  collectively  compact  operator  approximation  theory  implies 
that  if  N   exists  then  Nn~   exists  for  sufficiently  large  values  of  n  and 
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■1    ,  „-l 


conversely.   Thus  we  may  assume  that  both  N   and  N   exist  and  further  that 

n 

N   is  bounded  uniformly  in  n  for  sufficiently  large  n.   Hence  we  have  the 
convergence  theorem: 
Theorem  4.3.7 

The  solution,  cp  ,  of  the  approximate  equation  (4.3.12)  converges 
to  the  solution,  cp ,  of  the  singular  integral  equation  (4.2.1) 
Proof: 

From  (4.3.1)  and  (4.3.12) 


cp  -cp   =  N  g  -cp 
Tn  T     n  &n 


N_1(g  -N  cp) 
n  v  &n  rT ' 


-  N^1(gn-g+Mp-NnCp) 


Hence 

*  B{||gn-g|lc  +  ||  (Q-Q")o|!c)  (4.3.13) 

-'0  as  rr*00  by  lemmas  4.3.3  and  4.3.6,  where  B  is  the 

uniform  bound  for  N 

n 

Thus   [cp    }    converges    to   the    solution  of    (4.3.1),   has    the    required 
T  n 

Holder   property  by    lemma  4.2.1,    and   thus    converges   also    to   the    solution  of 
the   equivalent   singular    integral  equation    (4.2.1)    by   the   equivalence    theorem 
(section  4.2.1).      Since   q£H(v,v')    the   sets   [q   },{q   }    are   both  equi-H(M-)    and 
thus   convergence   of    the    approximate    solutions    in  H(M>)    follows   as    in    lemma   3.4.3 

Note   that   the   general   theory  gives   computable   error  bounds   for 
||cp   -<p||  .      These    are   normally   obtained   from  the   error   formula   for    the    particular 
quadrature    rule   but    it   has   been   shown   in   the    preceding    lemmas    that    the   error 
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: 
bounds  are  considerably  more  complex  in  this  case.   Bounds  for  (4.3.13) 

could  be  obtained  from  the  estimates  in  lemmas  4.3.6,  4.3.2  and  4.3.1  but  it 

would  not  seem  to  be  practical  to  try  to  estimate  these  constants.   At  best 

the  order  of  convergence  could  be  estimated  as  0(h  )  where  \    =  min(p,/2  ,v-v  ' ) 

thus  convergence  would  be  expected  to  be  slow  and  some  method  of  accelerating 

convergence  would  seem  advisable. 

As  usual  in  the  general  approximation  theory  discussed  previously, 

the  approximate  equation  (4.3.12)  is  not  solved  directly  but  rather  the 

equivalent  algebraic  linear  system 

0  0  1    r> 

N  cp    (t    .)=■(&   (t    ,)-b   (t    .))cp    (t    .)  +  .i-Q  cp    (t    .)   =  g   (t    .)  (4.3.14) 

nTnv   nj  v      v   ny  v   ny /Tnv   ny  2xnrn     nj  nv    nj  v 

j   =■   1,2,.  .  .  ,n. 

Recall  that  the  nodes  are  taken  on  the  arc  coordinate.   To  illustrate  the 
nature  of  this  system  consider  the  quadrature  approximation 


f  (s)  =  |  -S dt  = 

nv  '        v       t-s 


n-1 
E  i 
k=0 


-.  p  f(t)       ,.Lf(Vk^ 


t-s 


-dt 


where  cp,  (t)  is  the  piecewise  linear  function  defined  in  Section  3.5.   Suppose 

s=t . .   The  i    term  in  the  above  sum  vanishes  since  it  is 
J 


f(t.)       ; 

h       I  . 


fcj  t-t  /j+1     t       -t 

-E^iat  +  j        -Jii-dt 

J-l  J 


f  (t .)    f      ph  ph 

— ~-  |-    j      — ^du  +  j      — -duj     =   0,    by   the    substitution  u  =   t.-t    in  the 

0  0 
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first    integral   and    the    substitution  u  =   t-t.    in   the    second    integral.      It    is 
convenient   to   consider   the    remaining   terms    in  pairs,    symmetric   about   the    j 
term.      For   example    the    first   pair    is 


f(t      x)    rp    J-l   t-t  ?j   t-t  f(t  )      ,j+l      t-t 

h+2  t-t 

+  J  ^?t~dt 

t  J 

j+l 


Similar   substitutions    to   those   above    reduce    this    term   to 


{f(t.+  1)-f(t._1)}{(21og2-l)+l} 


Thus  we  obtain  a  scheme  of  "symmetric  pairing"  similar  to  that  obtained  by 
Bareiss  and  Neuman ,  mentioned  previously.   Note  that  in  the  present  case  the 
interval  of  integration  is  the  closed  contour  L  and  thus,  whether  in  this 
form  or  in  the  equivalent  parametric  form  considered  by  Gabdulhaev,  the 
functions  may  be  considered  as  periodic  over  the  range  of  integration.   Thus 
we  may  make  the  convention 


'(t.p  ■  f<Vj>  "  '<'a«-j> 


Thus  we   may  write 
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f  (t.)  = 


f  ^f(t)         r  i 


n 


-_dt  nf(t   )-f(t   )|  2  log  2 

L     J  J       J 


max(j,n-j) 
k=2 


Substituting  from  this  expression  into  the  system  (4.3.14),  where  the  term 

n 
involving  Q   is  evaluated  from  (4.3.3)  and  (4.3.4),  we  obtain  a  linear 


system  to  be  solved  for  cp  (t  .)  ,  j  =  l,...,n,  which  involves  only  function 
evaluations  of  the  given  functions. 

It  is  noted  that  this  system  is  considerably  more  complicated  than 
that  proposed  by  Gabdulhaev,  but  the  previous  analysis  shows  that  even  in 
this  more  complex  system  only  slow  convergence  is  guaranteed.   However  it 
is  assured  that  the  solution  does  converge  to  the  correct  result  and  this 
does  not  follow  for  Gabdulhaev's  system. 

4.4  Concluding  Remarks 

The  original  aim  of  this  investigation  was  simply  to  attempt  to 
extend  Anselone's  approximation  theory  to  the  case  of  linear  singular  inte- 
gral equations.   However  it  was  clear  that  the  occurrence  of  the  principal 
value  integrals  posed  special  problems  and  the  collectively  compact 
operator  approximation  theory  was  investigated  in  order  to  determine  whether 
this  might  be  extended  to  singular  integrands. 

The  Hilbert  transform  exists  in  the  space  of  Holder  continuous 
functions  but  approximations  of  such  operators  obtained  by  piecewise  inter- 
polation of  the  operand  do  not  necessarily  converge,  even  pointwise,  in  the 
Holder  space.   This  result  poses  a  serious  problem  for  a  direct  approximation 
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method   and   an  analysis   of   Gabdulhaev's   direct   method   shows    that  his    investi- 
gation of   the   convergence   of    the    solutions   obtained  by   polygonal   approxima- 
tions   cannot  be    considered   as    satisfactory.      This    argument    is    supported  by   an 
investigation  of   the   general    theory   of   approximation  methods   of  Kantorovich 
and  Akilov  which   again   shows    that   the   required  conditions   are    not   satisfied. 

Hence    it    is    concluded   that    it    is   necessary    to   follow  the   method  of 
theoretical    investigation   and   analytical   solution  of   this    singular    integral 
equation   and   reduce    it   to    the   more    tractable    case   of   a  weakly   singular 
Fredholm  equation.      The    theory   of   this    reduction    is  well  known   from   the 
standard   texts,    although   it  was    found   necessary   to    slightly   strengthen   the 
usual  hypothesis   on   the    density   function   in   order    to   obtain  a   consistent 
equation.      Finally    it    is    shown   that    the    approximate   solutions   obtained   from 
an   approximation  method   based   on   polygonal   approximations   of    the   numerator 
of    the    reduced  equation  do   converge    in   principle   but   more    precise   computa- 
tional  details    are    left    to   a   subsequent    investigation   of   some   particular 
example . 

Apart   from  computational  experiments   with    the    linear   system 
suggested    in   the    previous    section,    an   important   extension   to    this    investiga- 
tion would  be    to    the   case   where    the    interval   of    integration    is   a  non-closed 
are  or   finite    interval  on   the   real    line.      In   this   case   the   principal   value 
integral   preserves    the   Holder   condition   in    the    interior   of    the    interval   by 
the    Privalov   theorem  but    the    transformed    function  may  become   unbounded   at 
the   end-points.      Equations   of    this    type    are    solved   analytically  by 
Muskhelishvili    and   Pogorzelski   and   the   nature   of    the   unboundedness   of 
the   solution   at    the   end-points    can  be   bounded.      Thus    it  may  be    possible    to 
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introduce    a   special  norm  or   transformation   to    take   account   of    this   effect 
and   possibly  obtain   a   "uniform"   approximation  within   these    limitations. 

Other  generalizations  posing  additional  problems  include  systems 
of  linear  singular  integral  equations,  multi-dimensional  singular  integral 
equations    and  non-linear   equations. 


Ill 
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